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ABSTRACT 

Six  assumptions  concerning  sedimentation  in  the  east  equatorial 
Pacific  are  discussed  or  tested  by  computer  modeling  or  error 
analysis.  The  assumptions  considered  are:  (1)  depth  and  latitude  of 

deposition  are  the  prime  independent  variables  controlling  sedimenta- 
tion, (2)  rates  of  sedimentation  have  been  constant  since  the  end 
of  the  Eocene,  (3)  erosion  and  resedimentation  are  unimportant,  (4) 
fertility  and  dissolution  do  not  vary  with  longitude,  (5)  the  motion 
of  the  Pacific  plate  is  known  and  has  been  constant  since  the  Eocene, 
and  (6)  errors  in  measurement  of  sediment  properties  are  insignifi- 
cant. 

Two  numerical  models  are  used  in  the  test--  a forward-tracking 
model  (ESP)  and  a back-tracking  model  (ROPH).  Model  P.OPH  determines 
the  paleolatitude  and  paleodepth  of  a site  occupied  by  the  Deep  Sea 
Drilling  Project  (DSDP)  by  rotating  the  site  back  around  a pole  of 
rotation  and  by  using  a crustal  subsidence  curve  isostatically  cor- 
rected to  predict  depth.  The  paleodepth  and  paleolatitude  of  each 
site,  so  determined,  is  plotted  on  a grid  with  coordinates  of  depth 
and  latitude. 

The  sedimentation  rate  determined  from  DSDP  coring  records  and 
corrected  for  compaction  is  plotted  along  the  corresponding  track 
on  the  depth-latitude  grid.  The  carbonate  percentage  is  similarly 
plotted  along  its  track  on  a separate  grid.  Contouring  these  two 
grids  shows  the  extent  to  which  sedimentation  in  the  east  equatorial 
Pacific  can  be  explained  by  attention  only  to  depth  and  latitude. 
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The  contours  are  complexly  contorted.  Error  analysis  suggests  that 
the  contours  will  not  be  significantly  smoothed  by  assuming  different 
poles  or  rates  of  rotation.  Comparison  of  tracks  suggests  that  re- 
sedimentation or  erosion  in  preference  to  longitudinal  changes  in 
fertility  or  dissolution  can  explain  a portion  of  the  variation, 
which  caused  the  complexity  of  contours.  The  remaining  variation 
must  be  caused  by  temporal  dependency  or  measurement  errors. 

The  second  part  of  the  experiment  (ESP  or  forward-tracking) 
tests  for  temporal  dependency  by  working  the  back-tracking  process 
in  reverse,  i.e.  predicting  thicknesses  of  sediment  for  the  DSDP 
sites  by  working  the  process  in  reverse.  All  variation  is  assumed  to 
be  caused  by  temporal  dependency,  a correction  is  derived  for  this 
dependency,  and  ESP  is  executed.  Comparison  of  observed  thicknesses 
(P.0PH)  with  predicted  thicknesses  (ESP)  is  fairly  good,  which  shows 
that  depth,  latitude,  and  time  are  the  prime  variables  that  control 
sedimentation  in  the  east  equatorial  Pacific. 

This  study  suggests  that  only  long  period,  high  amplitude  changes 
in  sedimentation  rate  can  be  recognized  by  treatment  of  current  DSDP 
data  and  that  errors  in  measurement,  although  potentially  trouble- 
some, contributed  only  a small  amount  of  variation  to  the  prediction 
of  sediment  thicknesses. 

The  likely  magnitude  of  errors  of  measurement  as  they  influence 
track  elevation  and  latitude  is  given  quantative  expression,  as  is  the 
effect  of  changing  the  pole  of  rotation.  The  effect  of  errors  in 
estimate  of  basement  age  is  also  discussed. 
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INTRODUCTION 


Ewing  et  al_.  (1968),  using  normal-incidence,  seismic  reflection 
records,  demonstrated  that  sediment  is  thickest  at  about  4°  north 
latitude  in  the  east-central  Pacific  throughout  a belt  that  roughly 
paralles  the  equator.  These  sediments  are  composed  mainly  of  the 
calcareous  and  siliceous  remains  of  plankton.  Today  plankton  pro- 
duction is  greatest  near  the  equator,  but  only  a small  fraction 
of  the  plankton  shells  pass  into  the  sedimentary  record.  For  the 
calcareous  tests,  which  make  up  the  bulk  of  the  sediment,  preserva- 
tion is  related  to  depth  of  water--the  deeper  the  water,  the  more 
complete  the  dissolution  of  the  carbonate.  Destruction  takes  place 
both  in  the  water  column  and  on  the  sea  floor,  but  mainly  on  the 
sea  floor.  For  "normal"  crust,  the  depth  of  water  is  itself  a function 
of  the  age  of  the  crust,  water  depth  increasing  with  age  as  the  crust 
is  displaced  from  the  rise  crest  (Sclater,  Anderson,  and  Bell,  1971). 

As  a further  complication,  the  northward  motion  of  the  Pacific  plate 
causes  the  rise  crest  to  move  across  the  zone  of  high  productivity  at 
the  equator.  Thus,  sedimentation  can  be  described  as  the  product  of 
two  first-order  factors:  depth  of  water  and  proximity  to  the  zone 

of  high  productivity  (Winterer,  1973;  van  Andel  and  Heath,  1973). 

The  Deep  Sea  Drilling  Project  (DSDP)  completed  a number  of  holes 
in  the  thick  sedimentary  section  of  the  east-central  Pacific,  allowing 
rates  of  deposition  to  be  determined  for  the  past  40  million  years 
(Winterer,  et  aJL  , 1971;  Tracey,  et  aj . , 1971;  Hays,  et  aj. , 1972; 
van  Andel,  et  a_l.,  1973).  Employing  a process  he  termed  "backtracking", 
Berger  (1973)  traced  the  depth  of  selected  DSDP  sites  in  the  equatorial 
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Pacific  back  through  the  past  50  million  years,  using  subsidence 
data  of  Sclater,  Anderson,  and  Bell  (1971),  who  related  depth  of  the 
crust  to  the  age  of  the  crust.  To  facilitate  his  analysis  and  to 
test  for  temporal  changes  in  sedimentation  patterns,  Berger  divided 
the  50  million  years  in  9 parts,  each  part  represented  by  a grid 
that  had  as  its  abscissa  latitude  and  as  its  ordinate  depth  of  water. 
On  the  first  grid,  which  covered  the  Quarternary , he  recorded  the 
kind  of  sediment  that  drilling  records  showed  to  have  been  deposited 
during  that  time,  using  as  coordinates  the  present  depth  of  water 
and  the  present  latitude  of  the  drilling  site.  The  other  grids 
were  prepared  in  the  same  way,  using  the  projected  paleolatitude  and 
paleodepth  of  the  sites  for  the  appropriate  times. 

The  Quarternary  grid  showed  that  carbonate  sediments  tended  to  be 
symmetrically  arranged  around  the  equator,  the  CCD  (carbonate  com- 
pensation depth  or  depth  of  the  transition  from  calcareous  to  non- 
calcareous  sediment)  being  deepest  at  the  equator  and  rising  sharply 
to  the  north  and  south.  All  grids  delineating  older  sedimentation 
patterns  showed  the  same  symmetry,  but  had  maxima  that  were  shifted 
progressively  north,  which  is  to  be  expected  if  the  Pacific  plate  has 
been  displaced  northward  across  an  equatorial  zone  of  high  biologic 
productivity.  Using  a rate  of. 23  deg/my  (Winterer,  1973)  to  shift 
each  plot  southward,  Berger  (1973,  Fig.  10  and  11)  demonstrated  that 
in  a general  way  the  facies  patterns  younger  than  Eocene  were  con- 
gruent and  were  centered  on  the  equator.  Berger  further  suggested 
that  acomposite  grid  similar  to  his  Figure  11  might  be  used  to 
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predict  thicknesses  and  kind  of  sediment  that  overlay  any  portion 
of  the  crust  in  the  equatorial  central  Pacific  for  sediment  younger 
than  toceric  by  working  the  process  in  reverse  (forward-tracking), 
that  is,  by  tracing  the  point  back  across  the  lati tude-depth-of-water 
grid  to  determine  the  thickness  and  kind  of  sediment  that  might  be 
encountered,  were  a hole  to  drilled. 

A forward- tracki ng  model  based  on  composite  grids  has  6 assump- 
tions as  a foundation:  (1)  depth  and  latitude  of  deposition  are  the 

prime  independent  variables,  (2)  rates  of  sedimentation  have  been 
constant  since  the  end  of  the  Eocene  (37.5  mybp),  (3)  erosion  and 
resedimentation  are  unimportant,  (4)  fertility  and  dissolution  do  not 
vary  with  longitude,  (5)  the  motion  of  the  Pacific  plate  is  known  and 
has  been  constant  since  the  Eocene,  and  (6)  errors  in  measurement  of 
sediment  properties  are  insignificant. 

The  object  of  this  study  is  the  evaluation  of  these  assumptions 
by  performing  a numerical  experiment  using  DSDP  data  from  the  eastern 
equatorial  Pacific.  The  experiment  consists  of  four  parts.  In  the 
first  part  we  backtracked  23  DSDP  sites  on  2 grids  with  paleolatitude 
and  paleodepth  as  coordinates  (one  grid  for  rate  of  sedimentation,  one 
for  percentage  of  carbonate),  using  a pole  and  rate  of  rotation  to 
predict  latitude  and  a subsidence  curve  (Sclater  and  Detrick,  1973) 
to  predict  depth.  The  results  of  backtracking  were  contoured,  which 
allowed  a judgement  to  be  made  concerning  assumption  1 and  the  com- 
bined effects  of  assumptions  3 and  4. 

In  the  second  part,  the  same  23  sites  were  forward-tracked  through 


the  smoothed  fields  of  the  2 grids,  which,  in  effect,  built  23 


stratigraphic  columns.  By  comparing  the  differences  between  the  23 
constructed  columns  and  the  23  observed  DSDP  columns,  a measure  of 
the  validity  of  assumption  2 was  obtained. 

In  the  third  part,  the  importance  of  errors  in  the  estimate  of 
rate  and  pole  of  rotation  (assumptions)  and  errors  in  measurement 
(assumption  6)  are  given  quantative  expression  and  the  consequence 
of  these  errors,  as  they  influence  the  experiment,  is  discussed. 

In  the  fourth  part,  the  thickness  of  sediment  younger  than 
Eocene  is  generated  by  forward-tracking  for  the  east-central  Pacific 
and  compared  with  the  observed  thickness  as  determined  by  seismec 
investigations  (van  Andel , et  al_. , 1975,  adopted  from  Ewing,  et  al_. , 
1968).  This  exercise  provides  a display  of  the  over  all  effectiveness 
of  the  forward-tracking  method. 

COMPUTER  MODELS 

To  perform  the  experiments,  two  computer  models  were  developed. 
The  first,  ROPH,  backtracks  DSDP  sites  across  2 grids,  each  having 
latitude  and  depth  of  water  as  coordinates.  On  one  grid,  the  rate 
of  sedimentation  adjusted  for  compaction  is  plotted  in  1-million  year 
increments,  on  the  other,  the  percentages  of  carbonate.  Appendicies 
I and  III  document  the  program. 

The  second  model,  ESP,  forward-tracks  a point  across  the  grids 
generated  by  program  ROPH  to  predict  thickness  and  kind  of  sediment. 
Appendicies  II  and  III  document  program  ESP. 
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ROPH  corresponds  to  retracing  the  path  of  the  seafloor  through 
time  in  order  to  relate  the  kind  and  quantity  of  deposited  sediment, 
as  observed  in  DSDP  cores,  to  the  prevailing  latitude  and  depth  of 
the  sea  floor;  ESP  is  the  reverse.  Starting  with  a paleodepth  and 
paleolongitude,  the  kind  and  quantity  of  sediment  is  predicted  as  a 
function  of  paleolatitude  and  paleodepth. 

ROPH  Experiment 

The  latitude,  longitude,  and  water  depth  of  a DSDP  hole  provides 
the  starting  data  for  subsequent  incremental  rotation  of  the  point 
back  through  time  and  up  the  subsidence  slope.  At  each  step  the 
new  latitude  and  depth  of  water  determine  the  coordinates  of  the  point 
on  two  grids,  both  of  which  have  20°N  to  20°S  latitude  as  abscissa  and 
2300  m to  6300  m depth  of  water  as  ordinate.  The  rate  of  sedimentation, 
corrected  for  compaction,  that  prevailed  throughout  the  incremental 
rotation,  as  determined  from  the  data  in  the  DSDP  volumes,  is  plotted 
on  one  grid,  the  percentage  of  carbonate  on  the  other.  In  the  DSDP 
volumes  sediment  with  a carbonate  percentage  greater  than  75  percent 
usually  will  be  described  as  carbonate,  75  to  30  percent  as  alternating 
beds  of  calcareous  and  siliceous  ooze,  and  less  than  30  percent  as 
either  clay  or  siliceous  ooze;  thus,  the  percentages  have  a qualitative 
significance  with  respect  to  sedimentary  facies. 

Age-thickness  curves,  either  published  in  DSDP  volumes  or  con- 
structed from  the  data  of  DSDP  volumes,  provide  the  raw  material  for 
determination  of  sedimentation  rates.  Table  I lists  DSDP  holes  used 
in  this  model.  In  increments  of  1 million  years,  the  thickness  of 


Site 

(DSDP  No. ) 

Latitude 

(degrees) 

Longi tude 
(degrees) 

Depth 

(Corrected  meters) 

Sediment 

(m) 

Basement 
(m.y.b.p. ) 

65 

4.35N 

1 76 . 99W 

6130 

400 

130.0 

66 

2.39N 

166. 12W 

5293 

193 

95.0 

63 

16. 7 2 N 

T04.17W 

5467 

370 

100.0 

69 

6.00N 

152.86W 

4973 

360 

87.0 

70 

6.34N 

140.36W 

5059 

400 

60.0 

71 

4.47N 

140. 31W 

4419 

587 

60.0 

72 

0.44N 

138.87W 

4326 

363 

57.0 

73 

1.91S 

137.47W 

4387 

317 

54.0 

74 

6.23S 

136.08W 

4431 

103 

50.0 

75 

12.51S 

134.67W 

4181 

82 

37.5 

77 

0.48N 

133. 23W 

4291 

483 

41  .0 

78 

7.95N 

127.36W 

4363 

321 

35.0 

79 

2.55N 

121 . 57W 

4566 

413 

24.1 

80 

0.96S 

121 . 55W 

4399 

199 

22.5 

82 

2.59N 

106.94W 

3689 

223 

7.0 

159 

12.33N 

122. 29W 

4434 

108 

24.0 

160 

11.72N 

130.38W 

4940 

114 

34.8 

161 

10.24N 

139.95W 

4939 

245 

46.0 

162 

14.88N 

144. 04W 

4854 

153 

53.0 

163 

1 1 . 25N 

150. 29W 

5230 

284 

76.0 

164 

13.20N 

161 . 66W 

5485 

274 

115.0 

166 

3.76N 

175.08W 

4950 

310 

120.0 

168 

10.66N 

173. 55W 

5420 

246 

110.0 

Table  I.  Deep  Sea  Drilling  Project  Sites 

Latitudes,  longitudes,  depth  of  water,  sediment  thickness  overlying  the 
basement,  and  basement  age  of  drilling  sites  used  for  construction  of 
latitude-depth-of-water  grids. 


1 


sediment  that  is  younger  than  the  incremented  age  is  entered  in  the 
model  along  with  the  average  CaCO^  content  and  density  of  each  interval. 
Where  coring  is  not  continuous,  interpolation  is  used  to  determine 
missing  values.  Specifying  the  latitude,  longitude,  depth  of  water, 
age  of  crust,  thickness  of  sediment  cover,  rate  of  rotation,  and 
pole  of  rotation  readies  the  first  part  of  the  model  for  execution. 

Points  are  rotated  in  1-million  year  increments  counterclockwise 
around  the  pole,  using  standard  formulas  of  spherical  trigonometry 
to  obtain  new  positions.  The  new  latitude  provides  the  first  coordin- 
ate of  the  point  relative  to  the  latitude  -depth-of-water  grids. 

The  second  coordinate,  depth  of  water,  is  determined  by  reference  to 
the  Pacific  subsidence  curve  of  Sclater  and  Detrick  (1973).  At 
each  iteration  the  point  represents  a younger  stage  of  plate  evolu- 
tion. Accumulated  sediment  is  thinner"'  than  before,  and  the  crust  is 
younger.  The  subsidence  curve  predicts  the  depth  to  the  crust, 
stripped  of  sediment;  the  DSDP  hole  provides  an  estimate  of  depth 
of  crust,  overlain  by  sediment. 

The  difference  between  the  depth  predicted  by  the  Sclater-Detrick 
subsidence  curve  and  the  depth  observed  at  the  DSDP  site,  sediment 
cover  removed  and  the  crust  isostatical ly  adjusted,  is  held  constant. 

At  each  iteration,  the  new  depth,  according  to  the  subsidence  curve, 
plus  this  correction,  gives  the  depth  to  the  crust.  An  isostatic 
correction  using  the  reduced  thickness  of  sediment  gives  the  new  depth 
of  water.  Thickness  of  sediment  is  determined  by  subtracting  the 
amount  removed  for  all  previous  iterations  form  the  original  thickness 
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of  sediment.  The  density  of  sediment  for  each  interval  used  in  the 
calculation  is  the  value  recorded  in  the  DSOP  volume.  The  result 
is  a series  of  subsidence  curves,  each  corresponding  to  the  track  of  a 
drill  site  (Fig.  1). 

At  this  stage,  the  coordinates  of  the  grids  have  been  specified, 
and  the  rate  of  sedimentation  remains  to  be  calculated.  The  rate 
is  determined  by  calculating  the  thickness  that  was  removed  during 
the  1-million  year  rotation,  and  dividing  the  thickness  by  the  dura- 
tion of  the  interval,  in  this  model  1-million  years.  Bacause  compaction 
increases  with  depth  of  burial,  rates  determined  at  different  depths 
in  the  sedimentary  column  are  not  readily  comparable  and  must  be 
adjusted  to  some  standard  of  compaction.  The  standard  value  chosen 
is  arbitrary;  here,  all  thicknesses  are  adjusted  to  a porosity  of 
76  percent  (Hamilton,  1959;  Schlanger  et  aJL  , 1973),  using  bulk 

density  of  the  sediment  as  given  in  the  DSDP  volumes  and  adopting 
as  grain  density  the  value  2.60  gm/cm  . The  adjusted  thickness  is 
plotted  on  one  grid,  the  percentage  of  carbonate  on  the  other.  The 
pole  and  rate  of  rotation  (72°N,  83°W  at  0.81  deg/my)  is  approximately 
that  of  Jarrad  (1975).  For  isostatic  corrections,  the  crust  is  assumed 
have  a density  of  3.3.  Equations  and  a detailed  description  of  the 
method  is  contained  in  Appendix  I and  III. 


i 

V 

i « 

» * 


ESP  Experiment 

Model  ESP  predicts  the  stratigraphic  section  at  a point  on  the 
Pacific  plate  in  the  east-central  Pacific,  given  the  latitude,  longi- 
tude and  depth  of  water,  rate  and  pole  of  rotation  (same  as  ROPH), 


. JkL  k . -4 


72/'73 


162  160  161 


5000 


6000 


20°S 


Figure  1.  Latitude-depth  of-water  grid.  Traces  of  DSDP 
holes  are  shown  for  past  42  million  years,  with  depth  of 
water  and  latitude  as  coordinates.  Tick  marks  appear  at 
10  my  intervals;  OSDP  site  number  appears  at  lower  end  of 
each  trace;  pole  of  rotation  is  72°N,  83°W  at  0.8  deg/my. 
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and  two  grids  which  relate  sedimentation  rate  and  carbonate  percent- 
age to  latitude  and  depth  of  water. 

Because  ROPH  proceeds  by  counterclockwise  rotation,  one  might 
expect  ESP  to  use  clockwise  rotation,  but  this  is  not  the  case.  If, 
for  instance,  the  chosen  point  on  the  Pacific  plate  were  rotated  clock- 
wise, the  model  would  give  an  estimate  of  future  sedimentation  patterns 
on  the  corresponding  piece  of  sea  floor.  Accumulation  of  a historic 
record  could  be  simulated  by  rotating  the  point  counterclockwise  back 
to  its  position  of  origin  at  the  rise  crest,  and  then  clockwise, 
using  the  grids  to  preduct  the  thickness  and  kind  of  sediment;  how- 
ever, unnecessary  errors  will  result  from  this  procedure  if  the 
crust  is  older  than  about  42  my,  the  assumed  age  of  the  shift  in  pole 
position  for  the  Pacific  plate.  The  errors  occur  because  the  position 
of  the  pole  prior  to  42  mybp  is  imprecisely  known,  as  is  the  sedi- 
mentation pattern,  because  the  DSDP  holes  have  sampled  so  few  meters 
of  sediment  older  than  42  my  in  the  east-central  Pacific.  To 
illustrate  the  problem,  consider  a starting  point  where  the  crust  is 
80  my  old.  Rotate  the  point  counterclockwise  for  42  my  around  the 
Hawaiian  pole,  then  38  my  around  the  Emperor  pole  to  the  rise  crest. 
Because  the  position  of  the  Emperor  pole  is  imprecisely  known,  an 
error  in  latitude  results.  Now  rotate  the  point  clockwise  around  the 
pole  and  enter  the  grids  to  predict  the  thickness  of  sediment  laid 
down.  This  will  produce  an  error  in  the  estimate  of  depth  of  water 
because  sedimentation  pattern  is  not  known  for  the  first  38  million 
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years  of  rotation. 

The  solution  employed  here  uses  incremental  counterclockwise 
rotation  to  trace  the  path  of  the  point  back  through  time.  At  the 
end  of  each  increment  of  rotation,  the  grids  are  entered  at  the 
predicted  latitude  and  depth  of  water,  and  the  lithology  and  thickness 
for  that  increment  is  determined.  The  new  depth  of  water  is  calculated 
by  determining  the  isostatic  rebound  caused  by  the  removal  of  the  pre- 
dicted thickness  of  sediment.  Conceptually,  this  amounts  to  rotating 
the  point  counterclockwise,  stripping  off  a layer  of  sediment,  and 
determining  how  much  and  of  what  kind  has  been  removed. 

The  thicknesses  predicted  by  the  grid  has  a porosity  of  76  per- 
cent and  must,  as  a consequence,  be  corrected  to  reflect  compaction. 

The  correction  for  compaction  used  here  is  emperical,  being  derived 
from  the  curve  which  shows  the  relationship  between  the 
DSDP  thickness  and  the  DSDP  thickness  corrected  to  76  percent  porosity 
(Fig.  2.4).  The  derivation  of  the  relationship  is  discussed  at  Step 
1 5 of  Appendi x II. 

RESULTS  OF  ROPH 

The  backtracked  tracings  for  the  23  holes  drilled  in  the  east 
equatorial  Pacific  (Table  I)  are  illustrated  in  Figure  1,  and  the 
contoured  patterns  of  sedimentation  rate  and  carbonate  are  shown  in 
Figures  2 and  3.  Figures  4,  5 and  7 show  some  of  the  DSDP  site-tracks 
with  the  accompanying  rates  of  sedimentation  plotted  along  the  tracks. 
Figure  6 shows  some  of  the  carbonate  tracing. 


Figure  2.  Contours  of  rate  of  sedimentation.  Contours  are 
drawn  on  rates  of  sedimentation  for  past  30  my  plotted  along 
the  trace  of  each  site.  Figure  1.  Contours  are  in  meters/ 
million  years,  rates  corrected  to  porostiy  of  76  percent. 
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Figure  3.  Contour  of  percent  carbonate.  Contours  are  drawn 
on  percent  carbonate  for  past  30  my  plotted  along  the  trace 
of  each  site.  Figure  2. 
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Figure  4.  Traces  of  Sites  69,  70,  71,  72,  77,  79,  and  80 
showing  corrected  rate  of  deposition.  This  figure  is  an 
enlargement  of  a portion  of  Figure  1,  showing  rates  of  sedi- 
mentation adjusted  to  a porosity  of  76  percent.  Rates  greater 
than  5 m/my  are  rounded  to  the  nearest  meter.  Circled  and 
boxed  numbers  are  ages  in  millions  of  years  before  present. 
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In  an  ideal  world,  where  the  stated  assumptions  are  valid,  countour- 
ing of  the  backtracked  data  will  produce  isopleths  that  are  smooth, 
regular,  and  show  downward  deflection  at  the  equator,  where  rate  of 
sedimentation  is  greatest.  In  a crude  way,  the  contours  based  on  DSDP 
data  (Figures  2 and  3)  approach  the  ideal,  being  deflected  at  the 
equator,  but  complexities  abound.  Where  tracks  cross,  values  of 
sedimentation  rate  commonly  do  not  agree,  the  rates  of  some  neighbor- 
ing tracks  do  not  match,  and  some  isolated  tracks  can  be  accommodated 
only  by  bread  warpings  of  the  contours.  The  discrepancies  could  be 
caused  by  scaling  or  measuring  errors  of  original  data,  incorrect 
selection  of  pole  and  rate  of  rotation,  errors  in  basement  age,  errors 
in  the  subsidence  curve,  temporal  variations  of  sedimentation  rate, 
longitudinal  separation  of  points,  random  error,  or  local  variation 
i r.  rates  of  sedimentation.  The  results  presented  here  allow  some 
of  these  factors  to  be  evaluated. 

Sites  73,  74,  and  75 

In  the  portion  of  Figure  1 where  relationships  are  least  ambiguous. 
Sites  73,  74,  and  75  allow  smooth  contouring  and  fulfill  the  prophesy 
of  an  ideal  world,  which  leads  to  speculation  that  perhaps  south  of 
6°  latitude  and  above  4000  m,  sedimentation  on  the  Pacific  plate  may 
behave  with  probity. 

Sites  69  and  70 

At  the  crossing  of  their  traces,  Sites  70  (Figures  1 and  4)  shows 
sedimentation  rates  that  are  less  than  1 m/my  and  Site  69  shows  rates 


that  are  zero.  The  age  of  crossing  at  Site  69  is  2 mybp,  the  age  at 
Site  70  is  5 mybp.  The  crossing  must  be  considered  only  fair,  in 
spite  of  the  small  difference  in  rates,  because  the  time  of  crossing 
is  nearly  the  same  at  both  sites,  and  because  Site  70  clearly  received 
sediment,  while  69  did  not;  no  manipulation  of  reasonable  error,  as 
discussed  in  Appendix  IV,  for  any  of  the  modeled  variable  will  explain 
the  difference.  The  longitude  might  be  a factor,  because  approximately 
12  degrees  separates  the  two  sites,  but  local  variation  in  sedimentation 
rates  could  answer  just  as  well.  Both  holes  are  drilled  in  depres- 
sions on  the  sea  floor,  but  of  the  two,  the  depression  at  Site  70  is 
the  more  pronounced.  Within  2 miles  of  Site  69,  the  uppermost  layer 
is  eroded  (Tracey  and  Sutton,  et  , 1971,  p.  62).  At  Site  70  all 

layers  thicken  near  the  center  of  the  depression,  with  the  lowermost 
layers  showing  the  greatest  increase  throughout  its  history,  which 
suggests  that  Site  70  has  been  the  receptor  of  redeposited  sediment 
(Tracey  and  Sutton,  et^  aK , 1971,  p.  163).  A superabundance  of  sedi- 

ment at  Site  70  can,  and  probably  does,  explain  the  difference  between 
the  rates  of  the  two  sites. 

Sites  72,  77,  and  80 

The  tracing  of  Site  crosses  those  of  72  and  77  (Figures  1 and  4). 
Throughout  much  of  its  length,  rates  for  Site  80  are  one-half  to  one- 
fifth  those  of  72  and  77.  The  latitude  and  longitude  of  Sites  72  and 
77  differ  by  less  than  5 degrees  (Table  I),  and  so  both  will  trace 
similar  paths  across  the  grid,  regardless  of  the  poles  or  rate  used. 

Any  pole  or  rate  of  rotation  that  shifts  the  tracing  of  80  relative 
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Figure  5.  Trace  of  Sites  78,  82,  159,  160,  161,  163  showing 
corrected  rates  of  deposition.  This  figure  is  an  enlargement 
of  a portion  of  Figure  1,  showing  rates  of  sedimentation  ad- 
justed to  a porosity  of  76  percent.  Rates  greater  than  5 m/my 
are  rounded  to  nearest  meter.  Circled  numbers  are  in  millions 
of  years  before  present.  Solid  lines  mark  times  of  no  sedi- 
mentation. 
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to  77  and  72,  but  still  leaves  them  crossed  will  produce  an  unsatis- 
factory intersection  because  the  disparity  between  80  and  77  and  80 
and  72  is  so  great;  however,  could  a pole  or  rate  be  specified  that 
would  throw  80  to  the  right  of  72,  the  fit  would  be  satisfactary . 

Error  analysis  (Appendix  IV)  shows  that  such  an  alignment  can  be 
caused  only  by  assuming  rates  or  poles  of  rotation  that  differ  greatly 
from  any  of  those  previously  suggested  for  the  Pacific  plate  (Morgan, 
1972;  Clague  and  Jarrard,  1973;  Winterer,  1973;  Minster  et  al . , 

1974;  van  Andel , et  aJL  , 1975). 

Large  gaps  occur  between  cored  intervals  of  sediment  for  Site 
80.  According  to  Hayes  et  ah  , (1972,  p.  410),  large-scale  hiatuses 
may  be  present,  which  could  explain  the  low  sedimentation  rates,  but, 
if  present,  would  create  another  problem:  hiatuses  would  exist  where, 

according  to  the  assumptions  of  the  model,  none  should  be  expected. 
Thus,  Site  80  does  not  fit  the  assumptions  and  one  must  appeal  to 
longitudinal  differences  or  variation  in  sedimentation  rates  caused  by 
local  agents  to  explain  the  differences. 

Site  82 

Site  82,  (Figures  1 and  5)  records  a rate  of  103  m/my  at  6 m.y.b.p 
which  is  anomalously  high.  Paleontological  control  is  inadequate  for 
this  hole  and  this  can  explain  the  high  value;  however,  the  average 
percent  of  carbonate  is  only  73  at  2 m.y.b.p.,  77  at  4 m.y.b.p.,  and 
70  at6  m.y.b.p.,  which  is  low  compared  to  the  percentages  for  the 
nearby  tracing  of  Site  78  and  161  at  ages  younger  than  30  m.y.b.p. 
(Figure  6).  For  this  reason.  Site  82  must  be  judged  a poor  fit  within 
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the  grid. 

Sites  159  and  162 

At  the  crossing  of  Sites  159  and  162  (Figure  5).  Site  159  is 
7 m.y.  old  and  shows  a sedimentation  rate  less  than  1 m/my.  Site  162 
is  about  16  my  old  and  shows  a hiatus  from  the  present  to  31  m.y.b.p. 

The  crossing  is  unacceptable  in  an  ideal  framework,  causing  a long 
deflection  of  contours  in  Figure  2.  Here  no  change  in  pole,  rate  of 
rotation,  subsidence  curve,  or  age  of  basement  can  explain  the  dis- 
crepancy because  the  magnitude  of  error  and  its  persistence  through  time 
is  too  great.  Longitude  between  the  two  sites  differs  by  some  22°, 
which  might  explain  the  difference,  but  so  might  local  variations  in 
sedimentation  rate. 

Both  Site  159  and  162  were  drilled  through  a thin  sedimentary  layer 
that  overlies  a depression  in  the  basement. 

Sites  160  and  162 

These  two  sites  plot  as  adjacent,  nearly  parallel  traces  throughout 
much  of  their  length  (Figure  5).  Site  160  received  sediment  at  rates 
less  than  1 m/my  from  0 to  15  m.y.b.p.  and  then  increasing  to  10  m/my 
from  15  to  29  m.y.b.p.  Site  162  received  no  sediment  throughout  the 
last  31  my.  Thus,  as  with  Sites  162  and  159,  only  longitude  or  local 
variation  can  explain  the  difference,  which  is  unarguably  a significant 
one. 

Sites  78  and  J60 

Site  78  plots  close  and  nearly  parallel  to  Site  160  (Figure  5). 

Site  160  received  sediment  throughout  the  past  36  million  years;  Site 


Figure  6.  Traces  of  Sites  78,  82,  159,  160,  161  and  162  showing 
percent  of  carbonate.  This  figure  is  an  enlargement  of  a portion 
of  Figure  1,  showing  the  percent  of  carbonate.  Straight  line  seg- 
ments represent  hiatuses,  the  letter  "e"  a time  interval  for  which 
no  estimate  of  carbonate  is  available.  Circled  numbers  are  ages  in 
millions  of  years  before  present. 
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78  did  so  only  at  times  older  than  11  m.y.b.p.  Site  78  has  an 
anomalously  high  sedimentation  rate  for  ages  older  than  30  m.y.b.p. 

At  these  ages,  all  other  holes  show  a decrease  in  rate  of  sedimentation. 
Site  160  is  located  on  a down-faulted  block  in  an  area  where  the  relief 
of  the  sea  floor  changes  100  m or  so  in  a relatively  short  distance. 
Perhaps  the  lower  elevation  of  the  fault  block  allowed  redeposition 
of  sediment  eroded  from  the  surrounding  higher  areas.  Site  78  is 
located  in  an  area  where  the  sea  floor  is  smooth  (Hays  et  a 1 . , 1972 
p.  211).  Hays  et  al . (1972,  p.  217)  believe  that  nondeposition  rather 
than  erosion  is  responsible  for  the  lack  of  sediment  younger  than 
Middle  Miocene  at  this  site. 

DISCUSSION  OF  ROPH  RESULTS 

Many  of  the  discrepancies  noted  above  occur  where  sites  showing 
no  sedimentation  plot  close  to  sites  where  sedimentation  commonly  is 
less  than  1 m/my;  however,  some  of  the  sites  show  rates  that  exceed 
1 m/my,  approaching  10  m/my,  for  a period  of  8 million  years  (compare 
160  with  78,  162  and  161).  Although  the  quantitative  difference  is 
small  throughout  much  of  the  span,  the  difference  is  important  be- 
cause it  is  unambiguous  --  no  error  in  any  of  the  factors  entered  in 
the  model  can  explain  the  differences.  This  demonstrates  that  where 
sedimentation  rates  are  low,  errors  on  the  order  of  1 m/my  will  be 
common,  and  a question  arises  concerning  sites  where  sedimentation 
rate  is  high.  Should  discrepancies  be  expected  to  stay  near  1 m/my, 
or  might  not  the  error  be  proportional  to  the  rate  of  sedimentation  — 
the  higher  the  rate,  the  more  sediment  deposited,  and,  therefore. 
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the  greater  the  potential  for  error? 

Site  80  compared  with  77  and  72  suggests  that  greater  error 
might  be  expected  to  accompany  increased  rate  of  sedimentation.  Of 
all  the  sites  plotted,  only  73,  74,  and  75  behave  as  predicted  by  the 

' 

generalized  model.  Where  analysis  of  error  can  be  made  only  longitude 
or  local  variation  in  sedimentation  rates  can  conceivably  explain  the 
[ ' differences.  In  no  place  can  the  effect  of  longitude  be  discounted; 

however  in  many  places  the  data  indicates  that  local  variation  in 
sedimentation  rate,  probably  the  result  of  resedimentation,  i s important. 

Most  recently,  Berger  and  Johnson  (1975)  have  documented  slump- 
ing and  perhaps  creep  on  slopes  dipping  at  low  angles.  Seismic  re- 
flection profiles  commonly  show  that  sediment  is  thick  in  depressions 
and  thin  on  elevations  (Shor,  1956).  Most  of  the  DSDP  sites  considered 
here  are  drilled  in  depressions;  thus,  it  is  not  surprising  that  this 
generalized  model  of  sedimentation,  which  ignores  local  variation  in 
sedimentation  rates,  is  imperfect. 

ESP  RESULTS  AND  DISCUSSION 

After  contouring  of  the  sedimentation  rates  (Fig.  8)  and  the 

] percentage  of  carbonate  (Fig.  9)  had  been  completed  for  the  back- 

tracking data,  the  results  were  translated  to  two  40  x 40  arrays, 

which  were  used  by  the  computer  for  forward  tracking  (ESP),  taking 

I 

I as  starting  points  the  latitude,  longitude,  and  depth  of  water  of 

r * 

* the  23  DSDP  sites  (Table  I).  Several  executions  of  ESP  were  made; 

I i 

' • » after  each  execution,  the  arrays  were  altered  to  reduce  the  differences 

| * * 

r * 

between  the  observed  and  predicted  values,  always  maintaining  symmetry 


Figure  8.  Smoothed  contours  of  rate  of  sedimentation.  This  figure 
was  prepared  by  smoothing  the  contours  of  Figure  2.  Rates  are  in 
m/my,  depths  in  meters. 
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Figure  10.  Carbonate  variation  with  time.  The  ordinate  is  the 
quantity  (Pc  - 0C)/PC  where  Pc  is  the  predicted  percent  of  carbonate 
(ESP)  for  a DSDP  site  and  0C  is  the  observed  percent  (ROPH).  The 
abscissa  is  in  millions  of  years  before  present.  The  curve  is  con- 
structed by  connecting  the  mean  of  each  interval. 
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Figure  11.  Rate  of  sedimentation  variation  with  time.  The  ordinate 
is  the  quantity  (Pr  - 0r)/Pr  where  Pr  is  the  predicted  rate  of  sedi- 
mentation (ESP)  for  a DSDP  site  and  Or  is  the  observed  rate  (ROPH). 
Both  rates  are  corrected  to  76%  porosity.  The  abscissa  is  in  millions 
of  years  before  present.  The  curve  is  constructed  by  connecting 
the  means  of  each  interval. 
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of  the  contours. 

The  results  of  the  last  run  at  this  stage  of  the  investigation 
are  presented  in  Figures  10  and  11  as  percentage  differences  between 
observed  and  predicted  values. 

In  Figure  11,  age  of  sediment  in  millions  of  years  is  recorded 
as  abscissa,  the  difference  between  observed  and  predicted  sedimenta- 
tion rate  (corrected  for  compaction),  divided  by  the  predicted  rate, 
appears  as  the  ordinate.  Only  rates  of  sedimentation  greater  than 
5 m/my  are  plotted  in  the  figure  because  even  small  errors  in  the 
predicted  or  observed  rates  that  are  less  than  5 m/my  have  the  potential 
for  producing  large  errors  in  the  observed-minus-predicted-divided- 
by-predicted  ratio.  Such  errors  only  serve  to  obscure,  and,  thus, 
are  omitted. 

Figure  10  was  produced  in  a similar  way  using  the  observed  and  pre- 
dicted carbonate  percentages  and  restricting  the  analysis  to  percentages 
whose  sedimentation  rates  exceeded  5 m/my. 

In  Figures  10  and  11  the  mean  ratio  for  each  1 million  year  incre- 
ment has  been  plotted  and  the  successive  values  connected  by  line 
segments  to  form  a curve.  If  the  predicted  values  had  equaled  the 
observed  values  in  either  figure,  the  curves  would  have  plotted  as 
a straight  line  parallel  to  the  abscissa  and  intersecting  the 
ordinate  at  zero.  The  distance  from  zero  graphically  shows  the 
failure  of  ESP  to  reproduce  the  data  observed  from  the  DSDP  cores. 

These  results  clearly  show  that  depth  of  water  and  distance  from 
the  equator  cannot,  of  themselves,  explain  the  sedimentation  rates 
and  carbonate  percentages  observed  at  the  23  DSDP  sites,  which  is 
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Figure  12.  Mean  and  standard  deviation  of  the  differences  between 
predicted  thickness  (ESP  corrected  for  compaction)  and  observed 
thickness  (ROPH  comparted)  for  all  sites  at  6 million  year  intervals 
The  heavy  line  connects  the  mean  value  of  each  interval.  The  verti- 
cal bars  mark  off  1 standard  deviation  above  and  below  the  mean. 
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not  surprising,  because  Heath  (1969),  Berger  (1973),  and  van  Andel 
and  Moore  (1974)  have  suggested  that  the  carbonate  compensation  depth 
in  the  Pacific  has  temporal  dependency. 

In  the  succeeding  experiment,  all  of  the  variation  shown  in 
Figures  10  and  11  was  assumed  to  be  caused  by  temporal  changes  and  the 
predictions  of  thickness  for  the  23  i)SDP  sites  were  made  again  after 
applying  a correction:  The  ratios  traced  out  by  the  curves  (Figs. 

10  and  11)  are  applied  directly  to  the  predicted  rates  and  carbonate 
percentages.  For  instance  Figure  11  shows  that  at  23  mybp  the  sedimen- 
tation rate  was  54  percent  higher  than  the  mean  rate;  thus,  rates 
for  all  sites  of  this  age  will  be  increased  54  percent. 

The  result  of  this  experiment  is  summarized  in  Figure  12  and 
Table  II,  where  comparisons  of  total  accumulated  thicknesses  (DSDP- 
observed  minus  ESP-predicted)  is  made  at  6 million  year  intervals. 

Figure  12  records  the  mean  and  the  standard  deviation  of  the 
differences.  The  mean  difference  is  generally  less  than  6 m and 
the  standard  deviation  less  than  29  m.  Had  all  the  data  been  in- 
cluded in  the  construction  of  Figures  10  and  11,  the  subsequent 
corrections  to  the  rates  would  have  produced  a mean  difference  close 
to  zero  with  an  increased  standard  deviation. 

A fuller  appreciation  of  the  effectiveness  of  the  model  may  be 
realized  by  comparing  the  thickness  of  sediment  in  the  east-central 
Pacific  determined  from  DSDP  drilling  records  and  seismic  reflection 
profiles  with  the  thickness  predicted  by  the  model. 

Because  the  DSDP  results  suggest  that  the  rate  of  sedimentation 
decreased  drastically  for  times  older  than  Oligocene  but  younger  than  U. 


4. 


Figure  13.  Contoured  thickness  of  sediment  in  the  east  central 
Pacific  younger  than  36  my  determined  from  Deep  Sea  Drilling  Program 
records.  Contours  are  in  meters;  data  from  23  holes,  Table  1. 
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Figure  14.  Contoured  thickness  of  sediment  in  east  central  Pacific 
younger  than  36  my  predicted  by  model  ESP.  Contours  are  in  meters; 
predictions  made  at  5°  intervals  of  longitude,  2.5°  latitude. 
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Eocene,  36  my  was  chosen  as  an  arbitrary  but  reasonable  age  for  com- 
parison. In  Figure  13  the  observed  DSDP  thicknesses  of  sediment  younger 
than  36  mybp  (Listing,  Appendix  I)  has  been  plotted  and  contoured. 

The  predicted  thicknesses  were  generated  on  a grid  every  5°  of  longi- 
tude from  150°  W to  115°  W and  every  2.5°  of  latitude  from  15°  N to 
5°  S.  The  resulting  thicknesses  were  plotted  and  contoured  (Fig.  14). 

Comparison  of  the  two  figures  shows  that  almost  any  point  will 
show  a predicted  thickness  (ESP)  equal  to  or  slighly  greater  than  the 
observed  thickness  (DSDP)  and  that  the  center  of  maximum  predicted 
thickness  is  about  a degree  farther  north  than  the  center  of  observed 
maximum  thickness. 

Considering  all  of  the  sources  of  error  that  go  into  the  model, 
the  agreement  between  Figures  13  and  14  is  remarkably  good.  If  one 
accepts  the  agreement  as  satisfactory,  then  the  conclusion  follows  that 
depth,  latitude,  and  time  are, to  the  first  order,  the  variables  that 
describe  sedimentation  in  the  east  central  Pacific,  and,  therefore,  the 
curves  of  Figures  10  and  11  give  quantative  expression  to  the  variation 
of  sedimentation  with  time. 


I 


1 


CONCLUSIONS 


These  experiments  show  that  to  a first  order,  depth  of  water, 
distance  from  the  equator,  and  temporal  dependency  of  rates  and  carbonate 
percentage  are  adequate  for  the  prediction  of  sedimentation  patterns 
in  the  east  central  Pacific  Ocean.  The  agreement  between  the  predicted 
and  observed  thicknesses  is,  however,  not  perfect.  A portion  of  the 
difference  must  come  from  errors  of  measurement  of  sediment  properties 
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and  a portion  from  incorrect  assumptions.  Analysis  of  errors 
(Appendix  IV)  indicates  that  the  effect  of  errors  in  the  rate  or 
position  of  the  pole  of  rotation  will  cause  a relatively  small  error 
in  predicted  sediment  thicknesses. 

Longitudinal  effects  compared  with  the  effects  of  resedimentation 
and  erosion  are  insignificant  --  for  good  reason.  The  equator  and  the 
direction  of  spreading  are  subparallel  in  the  east  central  Pacific  so 
that  to  the  west,  away  from  the  rise  crest,  the  water  deepens  and 
dissolution  increases.  As  a result,  increased  dissolution  tends  to 
mask  changes  in  the  rate  of  sediment  production  in  the  surface  waters. 

As  an  example  of  the  effect  of  masking  by  dissolution,  consider 
a site  located  some  distance  west  of  the  rise  crest  where  the  water 
is  5000  m deep.  Assume  that  only  5 percent  of  the  sediment  produced 
in  the  surface  waters  passes  into  the  sedimentary  record.  If  the 
production  in  the  surface  waters  proceeds  at  the  rate  of  50  m/my, 
then  only  2.5  m/my  will  be  deposited  on  the  sea  floor.  If  the  rate 
increases  to  70  m/my,  only  3.5  m/my  will  be  deposited.  The  precision 
of  the  method  used  here  is  simply  not  good  enough  to  delineate  such 
small  differences,  and  thus,  the  effects  of  resedimentation  and  erosion 
can  be  expected  to  dominate  over  the  effects  of  longitude. 

The  sedimentation  rate  curve  (Fig.  11)  and  the  carbonate  curve 
(Fig.  10)  provide  a quantification  of  mean  sedimentation  rate  for  the 
east  central  Pacific.  From  2 to  16-18  mybp,  the  processes  of  sedi- 
mentation were  less  than  average  with  slight  increase  at  about  14  mybp. 
From  16-18  mybp  to  22-23  mybp  sedimentation  rates  and  carbonate  per- 
centages increased  to  greater  than  average  values  with  a climax  at 
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22-23  mybp.  From  22-23  mybp  to  about  40  mybp,  sedimentation  rate  and 
carbonate  percentage  started  a slow  decrease. 

These  rate  and  carbonate  curves  suggest  that  the  data  of  the  Deep 
Sea  Drilling  Project  are  useful  for  demarking  long  period  changes 
in  sedimentation  pattern  in  the  east  central  Pacific,  but  that  small 
scale  fluctuations  of  only  several  million  years  duration  (van  Andel 
et  al_. , 1975)  are  probably  beyond  the  resolution  of  the  data. 


APPENDIX  I 


Program  ROPH 
by 

Jerry  L.  Matthews 
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Software  and  Hardware 
Fortran  IV,  IBM  1800  computer. 

Purpose 

Given  the  latitude  and  longitude  of  a DSDP  hole  and  pertinent  data  on 
thickness,  age,  density,  and  carbonate  percent  of  the  sediment,  program 
ROPH  rotates  the  point  back  through  time  by  a sequence  of  iterations. 

At  each  iteration,  the  sedimentation  rate  or  the  carbonate  percent, 
depending  on  the  option  selected,  plots  on  a grid  with  latitude  as 
abscissa,  depth  of  water  as  ordinate. 

Variables 

Many  of  the  variables  appearing  under  the  In  column  are  printed  or  plotted. 
Only  variables  generated  by  the  computer  and  printed  or  plotted  appear  under 
the  Out  column. 

la 

MAPE  = 

I RITE  = 

IRATE  = 

JZAP  = 


Option  for  suppressing  the  printed  output.  MAPE  = 1 specifies 
output;  MAPE  = 0 suppresses  output. 

Option  for  suppressing  the  listing  of  input  data  and  depth 
correction  for  each  hole.  IRITE  = 1 specifies  output; 

IRITE  = 0 suppresses  output. 

Option  for  plotting  either  the  rate  of  deposition  or  lithology 
(expressed  as  percent  carbonate)  for  hole.  IRATE  = 0 specifies 
plotting  of  lithology,  IRATE  = 1 specifies  corrected  rate. 
Option  for  plotting  hole  number  in  conjunction  with  rate  of 
deposition  or  lithology.  JZAP  = -1,  the  hole  number  plots  with 
first  point  only;  JZAP  = 1,  the  hole  number  plots  with  all 
points;  JZAP  = 0,  suppresses  hole  number. 
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IG(M1)  = 

A = 

JJ  = 
BLAD1  = 


BLAD2  = 
BL0D1  = 


BLOD2  = 
TEST1  = 


TEST2  = 
TEST3  - 
RATE!  = 

RATE2  = 


RATE 3 = 


A 12  variable  array  used  for  conversion  of  floating  point 
to  alphameric  format.  Enter  at  0123456789.* 

Alphameric  data  printed  on  first  iteration  for  each  hole. 
Enter  as  A=NA  (abbreviation  for  not  applicable). 

Number  of  holes  to  be  read. 

Latitude  of  first  pole  of  rotation  entered  as  the  decimal 
fraction  of  degrees.  Use  north  projection  only  with  positive 
sign. 

Latitude  of  second  pole  of  rotation.  Conventions  same  as 
BLAD1 . 

Longitude  of  first  pole  of  rotation  entered  as  decimal 
fraction  of  degrees.  Use  positive  values  for  east 
longitude,  negative  for  west  longitude. 

Longitude  for  second  pole.  Conventions  same  as  BL0D1 . 

Age  (my)  for  change  from  first  pole  to  second  pole  of 
rotation  and  for  change  in  rate  of  rotation  from  RATE2  to 
RATE1 . 

Age  (my)  for  change  from  RATE3  to  RATE2. 

Age  (my)  for  change  from  RATE4  to  RATE3. 

Rate  of  rotation  (deg./m.y.)  for  first  pole;  determines 
rate  between  ages  0 and  age  TEST  3. 

Rate  of  rotation  (deg./m.y.)  for  first  pole;  determines 
rate  between  ages  TEST3  and  TEST2. 

Rate  of  rotation  (deg./m.y.)  for  first  pole;  determines 
rate  between  ages  TEST2  and  TEST! . 


BBI  I ' * 


RATE 4 = 
C,C2 ,C3, 
B( IN)  = 

ZLIM1  = 


MM(J)  = 

HOLE  = 
CLAD1  = 
CLOD1  = 
DEPTH  = 
AGE(B)  = 
AGE(M)  = 

DIST(M) 

FLITH(M) 

DENS(M) 

Out 


Rate  (Deg./m.y.)  for  second  pole;  determines  rate 
for  all  ages  greater  than  TEST1 . 

C4,C5,  and  C6  = six  variables  of  4 characters  each  that  incorporate 
heading  used  during  plotting. 

Seven  variables  that  incorporate  the  scales  for  the 
ordinate  and  abscissa  of  the  latitude  and  depth  of  water 
grid;  enter  as  -20.  0.020.030004005000. 

Entered  in  millions  of  years;  together  with  variable 
ZLIM2  provides  option  for  setting  limits  on  portion  of  hole 
to  be  plotted  and  printed;  points  whose  age  is  less  than 
ZLIM1  or  greater  than  ZLIM2  will  be  skipped. 

Number  of  points  for  each  hole  to  be  rotated  during 
execution . 

Alphameric  number  of  DSDP  hole;  used  by  plotter  and  printer 
Latitude  of  DSDP  hole  (same  convention  as  BLAD1). 

Longitude  of  DSDP  hole  (same  convention  as  BL0D1 ) . 

Depth  of  water  (m)  at  DSDP  site. 

Age  (mybp)  of  basement  at  DSDP  Site. 

Age  (mybp)  of  sediment  at  top  of  each  coded  interval  of 
DSDP  core. 

= Thickness  (m)  of  sediment  whose  age  is  less  than  AGE(M) 

= Lithology  (percent  carbonate)  of  interval  in  DSDP  core  of 
age  AGE(M)  to  AGE(M-l). 

= Density  of  interval  in  DSDP  core  of  age  (AGE(M)  to  AGE(M-l) 


Depth  of  crust  (m)  determined  from  subsidence  curve. 


PDEPH  = 
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DEPHC  = 

ADENS  = 
RATE  = 
HO  = 
URATE  = 
CERAT  = 
CLAD  = 
CLOD  = 


Depth  correction  (m)  carried  through  all  calculation  of 
water  depth;  measures  difference  between  calculated  and 
observed  elevation  of  crust. 

Weighted  mean  density  of  sediment  column,  computed  from 
densities  entered  from  DSDP  holes. 

• 1 

Sedimentation  rate  (m/my)  of  interval  uncorrected  for 
compaction.  »j 

Sedimentation  thickness  (m)  of  an  interval  corrected  for 
compaction. 

Sedimentation  rate  (m/my)  corrected  for  compaction  (calculated 
from  thickness  HO). 

Weighted  mean  rate  (m/m.y.)  of  sedimentation  for  sediment 
younger  than  AGE(M). 

Latitude  (deg.)  of  data  point  after  rotation;  negative 
values  are  south  latitude. 

Longitude  (deg.)  of  data  point  after  rotation;  negative 
values  are  west  longitude. 


Subroutines 

ALPHA,  ROCCW,  CRUST,  COUNT,  CPRIM,  XT AN , FT AN . 


Features 

Will  rotate  points  around  a first  pole  (Hawaiian)  and  a second  pole 
(Emperor).  Rotation  is  counterclockwise.  For  the  first  pole  three 
intervals  of  time  may  be  defined,  during  which  a different  rate  of  rotation 
applies.  Although  written  for  rotation  of  the  Pacific  plate,  the  program  is 
general  and  will  work  for  any  plate  requiring  counterclockwise  rotation,  as 
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long  as  north  projection  is  used  for  the  pole.  The  plotting  feature  may 
be  surpressed. 


J 


« 
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Input  of  data 

Images  numbered  1 through  11  require  1 card.  Variables  on  image  number  12 
will  normally  require  more  than  1 card. 


Image 


Number 

Read  Statement 

Format 

1 

MAPE 

(12) 

2 

I RITE 

(12) 

3 

IRATE 

(12) 

4 

JZAP 

(12) 

5 

IG(M1),  Ml  = 1,  12 

( 1 2A1 ) 

6 

A,  JJ,  BLAD1 , BL0D1 , BLAD2, 

BL0D2,  TEST1 , TEST2 , TEST3  (A2.13.7F10.2) 

7 

RATE4,  RATE3,  RATE2,  RATE1 

(4F10.2) 

8 

C,  C2 , C3,  C4 , C5,  C6 

(6A4) 

9 

B (IN),  IN  = 1,  7 

(7A4) 

10 

ZLIM1 , ZLIM2 

(2F10.2) 

11  MM(J),  HOLE,  CLAD1 , CL0D1  , DEPH,  AGEB  (16,  3X,  A3,  4F10.2) 

12  AGE(M) , DIST(M) , FLITH(M)  , DENS(M) , M = 1,  MMM  (2F10.1,  6X,  A4,  F10.0) 

Me  thod 

1.  Read  in  images  1 through  10. 

Date  that  is  pertinent  to  all  holes  is  read  in,  including  coded  values  that 
specify  plotting  and  printing  options,  headings  used  by  the  plotter,  and  poles 
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Figure  1.1.  Variables  TEST1 , TEST2,  and  TEST3  set  off  intervals  of 
time  during  which  the  rates  of  rotation  RATE1 , RATE2,  and  RATE 3 
govern  rotation  of  points  around  the  first  pole  (Hawaiian);  RATE4 
applies  to  second  pole  (Emperor),  and,  thus,  age  TEST3  marks  the 
change  of  pole  position. 


Figure  1.2.  Typical  sedimentation  curve  used  to  determine  thickness 
as  a function  of  age. 


and  rates  of  rotation.  For  the  first  pole  of  rotation  (BLAD1 , BL0D1 ) 
three  intervals  of  time  (TEST1,  TEST2  and  TEST3)  must  be  defined,  during 
which  three  rates  of  rotation  (RATE2,  RATE3,  and  RATE4)  apply.  RATE4 
applies  for  the  second  pole  (BLAD2,  BL0D2).  Figure  1.1  shows  the 
relationship. 

2.  Plot  and  annotate  grid. 

The  grid,  approximately  20  inches  by  17  inches,  has  latitude  as 
abscissa,  being  divided  in  40  parts  of  1 degree  each,  reaching  from  20°N 
to  20°S  and  depth  of  water  as  ordinate,  being  divided  in  28  parts  of  100  m 
each,  reaching  from  2700  to  5500  m;  however,  a wide  enough  margin  remains 
so  that  the  paper  will  accommodate  the  plotting  of  depths  as  shallow  as  2300  m 
and  as  deep  as  6300  m.  The  depth  is  annotated  at  3000,  4000,  and  5000  meters, 
the  latitude  at  20.0°,  0.0°,  and  -20.0°  (the  negative  value  indicates  south 
latitude).  The  heading  normally  specifies  the  span  of  time  (variables  C,  C2, 
C3,  C4,  C5,  C6)  for  which  the  plot  is  applicable.  Either  rate  of  sedimentation 
or  percent  of  carbonate  will  plot,  depending  on  the  option  selected  (variable 
I RITE ) . 

3.  Read  image  11 

This  is  a heading  card  that  specifies  the  number  of  points  that  will 
be  read  in  by  image  12,  and  also  includes  the  DSDP  hole  number,  the  latitude 
and  longitude  of  the  hole,  and  the  age  of  the  basement. 

4.  Read  image  12 

These  are  the  summary  cards  that  describe  the  data  of  the  DSDP  cores  and 
set  the  increment  used  during  the  rotation.  Variable  AGE(M),  Figure  1.2, 
is  determined  by  selecting  a set  of  ages  from  the  curve  (the  interval  between 
ages  need  not  be  constant).  Noting  the  thickness  of  sediment  for  each  aqe 
provides  the  values  for  variable  D1ST(M).  Variable  FLITH(M)  is  the  mean 


carbonate  percentage  and  DENS(M)  the  mean  density  of  the  interval  marked 
off  by  AGE(M)  and  AGE(M-l). 

5.  Determine  the  weighted  mean  density 

To  make  the  isostatic  correction  for  the  elevation  of  the  crust  at  the 
drilling  site  after  removal  of  the  sediment,  the  density  of  the  sediments 
must  be  estimated.  The  weighted  mean  density  provides  this  estimate: 

2 

1.5.1  5 = Z Di  T./ZT1 


where 

D = weighted  mean  density  (ADENS) 

D.  = interval  density  (DENS(M)) 

T.  = thickness  of  sediment  for  various  intervals 

6.  Determine  elevation  of  crust  from  subsidence  curve 

The  elevation  of  the  crust  is  determined  by  the  subsidence  curve  of 
Sclater  and  Detrick  (1973),  as  a function  of  age,  AGE(M).  The  descriptive 
parameters  are  contained  in  subroutine  CURVE.  Figure  1.3  shows  a schematic 
illustration  of  the  curve. 

7.  Determine  the  crustal  correction 

This  is  the  correction  that  adjusts  the  elevation  of  the  crust  at  the 
DSDP  site  to  compatibility  with  the  elevation  of  the  crust  predicted  by  the 
subsidence  curve.  Figure  1.3.  The  steps  involved  are:  a)  determine  the 
depth  of  the  crust  predicted  by  the  subsidence  curve  (Step  6),  b)  determine 
the  estimated  depth  of  the  crust  at  the  DSDP  hole  with  the  sediment  cover 
removed  and  isostatic  balance  restored,  and  c)  subtract  the  two  values 
(steps  a and  b)  to  determine  the  correction. 
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- Figure  1.3.  A schematic  representation  of  the  subsidence  curve 

f (Sclater,  Anderson,  and  Bell,  1971),  showing  the  effect  of  adding 

^ the  crustal  correction  (Eq.  1.7.2).  Curve  abc  is  the  original  sub- 

sidence curve;  curve  def  is  the  corrected  curve;  distance  eb  is  the 
correction  (Eq.  1.7.2). 


The  correction  is  expressed  by: 

1.7.1  W = [(Db  - D)]/(Db  - D$)]  Tm 

1.7.2  C = W + H - P 
where 

W = thickness  of  sediment  displaced  by  water  after  the  sediment  cover  has 
been  removed  and  the  crust  has  attained  isostatic  balance 
Db  = density  of  basaltic  basement 

D = weighted  mean  density  of  sediment  column  (Step  5) 

D$  = density  of  seawater 

Tm  = thickness  of  sediment  overlying  the  crust 

C = correction  for  difference  in  observed  and  predicted  elevation  of  crust 
H = depth  of  water  at  drilling  site 

P = depth  to  crust  predicted  by  subsidence  curve  (Step  6) 

During  subsequent  calculations  of  depth  of  water,  correction  C (Eq.  1.7.2) 
will  be  added  to  the  crustal  elevation  determined  by  the  subsidence  curve 
(Figure  1.3.).  The  effect  of  this  addition  is  to  produce  a second  curve 
parallel  to  the  subsidence  curve  and  separated  from  it  by  the  distance  C 
(Eq.  1.7.2).  This  second  curve  will  be  the  one  to  which  isostatic  corrections 
are  applied. 

8.  Rotate  the  point 

This  rotation  produced  the  new  latitude  of  the  point  which  is  plotted  as 
a coordinate  on  one  of  the  grids  (Fig.  1).  The  rotation  of  the  point  and  its 


subsequent  latitude  is  determined  by  solution  of  spherical  triangles  ABC, 
then  ABC1,  Figure  1.4,  given  2 sides  of  the  triangle  and  the  included 
angle.  Referring  to  Figure  1.4,  the  point  A is  the  spin  axis  of  the  earth, 
B the  pole  of  rotation,  and  C the  point  to  be  rotated  to  its  new  position 
C' . 

For  triangle  ABC  the  included  angle  at  CAB  is  the  absolute  difference 
between  the  longitude  of  the  point  at  C and  the  pole  at  B.  Angle  b is  90° 
minus  the  latitude  of  the  starting  point  C,  angle  c is  90°  minus  the 
latitude  of  the  pole  B.  Subroutine  ROCCW  makes  these  determinations  and 
controls  rotation  of  the  point.  The  two  sides,  b and  c,  and  the  included 
angle  ABC,  determined  by  ROCCW,  provide  the  substance  of  the  solution  for 
the  side  a and  the  angle  ABC  of  triangle  ABC,  Figure  1.4.  Because  rotation 
around  the  pole  at  B requires  the  point  C to  trace  out  a line  of  latitude 
around  the  pole,  side  a'  must  equal  side  a,  which  leads  to  the  solution  of 
triangel  ABC1,  using  sides  c and  a',  and  included  angle  ABC'.  Side  c has 
the  same  value  as  before  by  identity,  and  side  a'  is  equal  a.  The  included 
angle  ABC',  as  illustrated  in  Figure  1.4,  is  equal  to  angle  ABC  plus  angle 
CBC' , which  is  determined  in  subroutine  COUNT  by  multiplying  the  rate  of 
rotation  times  the  length  of  time,  AGE(M).  After  solution  of  triangle  C'AB 
by  subroutine  CPRIM,  the  latitude  of  point  C'  is  determined  from  angle  b', 
being  90°  minus  b'.  The  longitude  of  the  point  is  angle  CAB  minus  angle 
C'AB,  the  quantity  added  to  the  longitude  of  point  C.  These  conversions 
of  spherical  angles  to  latitude  and  longitude  are  made  by  subroutine  ROCCW. 
Solution  of  the  spherical  triangles  is  accomplished  by  subroutine  CPRIM 
using  the  following  equations: 
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Figure  1.4.  Great  circles  on  a sphere  that  form  spherical  triangles 
ABC  and  ABC'  used  for  rotation  of  points  on  a sphere.  Point  A Is 
the  spin  axis,  point  B the  pole  of  rotation,  and  point  C the  point 
to  be  rotated  to  position  C'. 
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tan  1/2  (C  + B)  = [cot  1/2  A cos  1/2  (c-b)]/cos  1/2  (c  + b) 

tan  1/2  (C  - B)  = [cot  1/2  A sin  1/2  (c  - b)]/sin  1/2  (c  + b) 

tan  1/2  a = [tan  1/2  (c  - b)  sin  1/2  (C  + B)]/sin  1/2  (C  - B) 

where  angles  b and  c are  the  sides  of  the  triangle  and  angle  A is  the 
included  angle. 

The  determination  of  angle  CBC ’ by  program  ROPH  is  not  so  straightfor- 
ward as  presented  above,  because  the  rate  of  rotation  can  be  made  to  vary. 
Figure  1.1  illustrates  this  point.  If  the  age,  AGE(M) , is  less  than  age 

TEST3,  the  angular  distance  is  RATE4  times  AGE(M) . If  the  age,  AGE(M) , is 

greater  than  age  TEST3,  but  less  than  age  TEST2 , the  angular  distance  is 
RATE4  times  TEST3,  plus  RATE3  times  the  quantity  AGE(M)  minus  TEST3.  Solutions 
to  the  two  remaining  possibilities,  greater  than  TEST2,  but  less  than  TEST3, 
and  greater  than  TEST3,  are  determined  by  analogous  procedures. 

For  the  arrangement  of  point  C and  pole  B illustrated  in  Figure  1.4, 
angle  ABC'  is  determined  by  addition  of  angle  ABC  and  CBD ' , but  if  point  C 
were  to  lie  counterclockwise  of  pole  B,  subtraction  is  appropriate.  Comple- 
cations  arise  if  the  longitude  of  point  C equals  the  longitude  of  point  B 
and  where  a rotation  swings  point  C from  a position  clockwise  to  B to  one 
counterclockwise  of  B.  Subroutine  COUNT  tests  for  these  possibilities  and 
uses  the  appropriate  algorithm  to  determine  angle  ABC'. 

To  summarize,  program  ROPH  determines  the  angular  distance  of  rotation 
(angle  CBC'),  subroutine  ROCCW  directs  the  rotation  by  calling  subroutine 
CPRIM,  which  solves  for  the  unknown  side  a and  angle  ABC  of  the  spherical 


triangle,  given  side  a and  c and  included  angle  CAB,  subroutine  COUNT 
calculates  the  angle  CBC'  from  which  angle  ABC'  is  determined,  and  sub- 
routine CPRIM  solves  for  sides  b and  angle  C'AB,  using  side  a1,  side  c, 
and  included  angle  ABC'.  Subroutine  ROCCW  converts  the  angular  measures  to 
latitude  and  longitude. 

9.  Determine  depth  of  water 

During  the  execution  of  the  program,  layers  of  sediment,  in  concept, 
are  shipped  from  the  pile  of  sediments  that  overlies  the  crust.  The  effects 
of  the  removal  of  the  sediment,  the  consequent  isostatic  rebound  of  the  crust, 
and  the  elevation  of  the  crust  predicted  by  its  decrease  in  age  (the  sub- 
sidence curve  is  a function  of  age)  must  be  calculated.  Depth  of  water  is 
determined  by  finding  the  crustal  elevation  from  the  subsidence  curve  (sub- 
routine CRUST),  adding  the  crustal  correction  (Step  7),  and  making  an  isostatic 
correcton  for  the  thickness  of  sediment  that  remains. 

This  may  be  expressed  as: 

1.9.1  P_  = P - W + C 
n s n 

where 

Pn  = new  depth  of  water 

P$  = depth  of  water  determined  from  subsidence  curve 

Wn  = thickness  of  sediment  displaced  by  water  (calculated  from 

equation  1.7.1,  where  Tm  is  equal  to  the  reduced  thickness  of 
sediment  and  D is  the  weighted  mean  density  of  this  sediment) 

C = correction  for  difference  in  observed  and  predicted  elevation  of 
the  crust  (Equation  1.7.2). 
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10.  Convert  latitude  and  longitude  to  positive  and  negative  values. 

Through  a sequence  of  tests  latitude  is  converted  to  a negative  value  if 

southern,  and  longitude  is  converted  to  a negative  value  less  than  180° 
if  western. 

11.  Plot,  or  plot  and  print  either  rate  of  sedimentation  or  percentage  of 
carbonate. 

Depending  on  the  options  that  are  set,  values  will  be  plotted,  or  plotted 
and  printed.  For  plotting  the  carbonate  percentages,  the  values  are  simply 
retrieved  from  storage.  If  the  rate  of  sedimentation  is  specified,  the 
uncorrected  rate  of  sedimentation  and  the  adjusted  rate  of  sedimentation 
(adjusted  to  rate  calculated  from  thickness  adjusted  to  76%  porosity)  are 


printed,  but  only  the  adjusted  rate  may  be  plotted. 

The  steps  in  computing  the  adjusted  sedimentation  rate  are  : a)  determine 

the  thickness  of  sediment  for  the  interval;  this  is  the  difference  between  two 


successive  thicknesses  entered  as  data  (Step  4,  Card  12),  b)  determine  the 
porosity  from  the  density  of  the  sediment  according  to  the  equation  of 
Hamilton  (1959), 


1.11.1 


D - D 
s c 


D - D 
w g 


where 


D = density  of  the  interval 

Dg  = mean  density  of  the  grains 

D = density  of  seawater 
w J 


N = porosity  of  sediment 


and  c)  transform  the  thickness  to  the  value  it  would  have  were  its  porosity  76% 


M.  v.  . 


by  the  equation  of  Schlanger  et  al . , 1973, 


1.11.2  Ha  = (1.  - N)/(l.  - .76))  Hq 

where 

N = porosity  of  sediment 
Hq  = unadjusted  thickness 

H = adjusted  thickness 

a 

and,  d)  determine  the  unadjusted  and  adjusted  sedimentation  rate  from  the 
appropriate  thickness. 

12.  Return  to  Step  8 
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ROl'll  LISTING 
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DIMENSION  DI  ST  ( 60)  , AGE  ( 60  ) , FI  IT  H(  60 ) , MM  ( 60 ) ■ R ( 2 0 1 . (1ENS  ( 60 ) , VE 
1L(60) , 1GI  12  ) 

DEFINITION  III-  VARIABLES - 

MAPfc option  fok  printing  uofut 

1 IF  OUTPUT  DESIRED 

0  IF  OUTPUT  NU I DESIRED 

I KITE  ... .OPT  ION  FOK  PRINTING  INPUT 

1  IF  PRINT 1NG  DESIRED 

0  IE  PRINTING  NOT  DESIRED 

IRATE .. ..OPT  ION  EUr  PLOTTING  RATE  OF  SEDIMENTATION  OR  LITHOLOGY 

1  IE  RATE  UE  SEUIMEN1  ATIUN 

0 IF  LITHOLOGY (PERCENT  CARBONATE) 

JZAP ..  . . .OPTIUN  FUR  PI.01TING  HOLE  NIIMREk  IN  CONJUNCTION  WITH  KATE 
OF  SEDIEEnTaT ION  OH  LITHOLOGY 
-1 IE  HOLE  NUMBER  WITH  FIRST  POINT  ONLY 

0  IE  NO  HOLE  NUMBER 

1  IE  HOLt  NUMBER  WITH  EACH  POINT 

1GIM1 ) . . . ALPHAMt RIL  ARRAY  USED  WITH  SUBROUTINE  ALPHA.  ENTER  AS  0 1 23A BG7 H9 . * 

A ALPHAMERIC  ABREV 1 AT  I ON  FOR  'NOT  APPLICABLE1.  ENTER  AS  NA 

JJ NUMBER  UE  HULLS,  CUNTRUI.S  MAIN  LOOP. 

BLAU1... .LATITUDE  Ur  RIkST  POLE  IN  DECIMAL  FRACTION  OF  DEGREES.  USE 
NORTH  RKUJEC110N 

BLAI)?. .. .L ATI TUUt  OF  SECOND  POLE  IN  DEGREES.  USE  NORTH  PROJECTION 
BLOUl ....LONGITUDE  OP  FIRST  POLE,  IN  DECIMAL  FRACTION  OF  DEGREES. 

+ . . . . I F EAST 
-....IE  WEST 

6L0U2. .. .LONG  I Tut  OF  SECOND  PULE.  SAME  CONVENTION  AS  PLOD  1 
TEST] ....AGE (MIL  LIU  NS  UE  YEARS)  FOR  SHIFT  FROM  FIRST  T()  SECOND  POLE,  FROM 
RATE  2 TO  KATE  I 

TEST?.. ..AGE (MILLIONS  UP  YEARS)  FOR  SHIFT  FROM  RATE3  TO  RATC2. 

TEST3. . . . AGE (MI LL IONS  UP  YEARS)  FOR  SHIFT  FROM  KATEA  TO  RATE3. 

RATE  1 ...  .RATE  OF  KU A T A i 1 ON ( OEG/M Y ) FOR  FIRST  POLE,  FROM  ZERO  YEARS  TO 
t ES  T 3 YPAKS. 

RATE?. .. .RATE  OF  KU T AT  ION ( DFG/MY ) FROM  TFST3  TO  TEST?  YEARS. 

RATE3... .RATE  (IE  RUT AT  I ON ( DEG/MY ) FROM  T F S T?  TO  TEST  1 YEARS. 

RATE*  ... .RATE  OF  Ru  T A T ION ( OF G/M Y ) FOR  AGFS  OLDER  THAN  TEST  1 YEARS.  GOVERNS 
ROTATION  A KUO NO  SECOND  PULE ( PL AD2 , BLOD?). 

C2-C6 SIX  ALPHAMERIC  VARIABLES  USED  AS  HEADINGS  BY  PLOTTER  FOR 

GRID.  E NT  f K AS  -20.  0 . 020. 030DA 0005 000 
71 1M1  ....  ITERATIONS  YUlINGEk  THAN  THIS  VALUE  NOT  PLOTTED  DR  PRINTED 
7LIM?.... ITERATION  OLDER  THAN  THIS  V/UUE  NOT  PLOTTED  OR  PRIN1E0 
MM I J I . . . .NOMBE K OR  PUlNlS  PER  HOLE. 

HOLE NUMBER  OB  HO l ( ( A L PHAMER  1 C I . 

CLAD1  ....LATITUDE  Ot  KILE.  (SAME  CONVENTION  AS  BLA01) 

CL0D1 ....LONGITUDE  OF  HOI E ( SAME  CONVENTION  AS  ftLUDl ) . 

DEPH. . . .DEPT H OF  WATER (METERS) . 

AGEB AGE  (IE  B A SEE  E Ml  I My  , . 

AGE (Ml... nGE  OE  SEDIMENT  AT  TOP  OE  INTERVALS  IN  CORE. 

DIST  IM ). .THICKNESS  OF  SEDIMENT  W) USE  AGE  IS  LESS  THAN  AGE ( M I . 

EL  I I H (M  | , I.  I T Hj  LO  GY  ( PERCE  NT  CARBONATE)  1IE  INTERVAL. 

DtNS (M) .. .DENSITY  UR  INTERVAL! IF  IMTFkVAL  HAS  HIATUS,  ENTER  DENSITY  OF  ZERO). 

9 EDRMiA  T II?) 

10  FORMAT ( 1 2 A 1 I 


- . * 


11  E0RMATIA2,  IS,  7F10.2) 

12  FORMAT  I4F  10.2  » 

13  FORMAT ( A2 , 13,  11F10.2) 

14  FORMAT  ( 6 A 4) 
lb  FORMAT ( 7A4 ) 

16  F0RMA1  ( 2F  10.2  I 

17  FORMAT  1 30X , 2F1S.2) 

18  F0RMA1 (2A4  ) 

19  FORMAT  I A4) 

20  FORMAT  116,  3X , A3,  4F10.2) 

21  FORMAT!//,  IX,  'Hue  NUMBER  =',A3I 

22  FORMAT (2F 10.1 ,6X , A4 , 2F 1 0 .0 > 

23  FORMAT ( IX, 'LATllUOt  IS,,F7.2,*  LONGITUDE  1S',F8.2,'  DEPTH  1S',F8. 
12,'  BASEMENT  AGE  IS',F7.21 

24  FORMAT ( IX,  'DEPTH  UF  CRUST(POEPH)  I S ' , F 8 . 2 , 2X , • DE PTH  CORREC T ION ( DEP 
1HC.)  I S ' , F H . 2 , ' AVERAGE  DENS  IT  Y I AOENS  ) IS',F6.2I 

25  FORMAT!/, 10X,' AGE  OF  POINT',  10X, 'DEPTH  IN  HOLE',  10X , ' L I THOLUGY  U 
1 F POINT',  BX  , 'DENSITY',  BX , 'VELOSITY'l 

26  FORMAT (6X,F10.2,  13X,  F10.2,  23X,  A4.18X,  F4.2,  19X,  F4.2) 

27  FORMAT!/,'  AGE  UF  DEPTH  UNCORC  COR  COR  COM  COR  LlTHOLO 
1GY  PRESL  PRES  WATER') 

28  FORMAT!'  PUINT  IN  HOLE  kATE  INTERVAL  RATE  RATE 

1 LAT  LUNG  DEPTH') 

29  FORMAT  ( F 7 . 1 ,F  F . 1 , A7 , 3A8  , A9  , F 1 0 . 2 , 2 F9  . 2 I 

30  FORMAT (F7. 1 , F8 . 1 , F 7 . 1 , 3F 8. 1 , A9 , F 1 0. 2 , 2F9 .2 ) 

31  FORMAT (SA1, A3  ) 

32  FORMAT  ISA1 ) 

33  FORMAT  <A4,  A3  ) 

34  FORMA T(A4) 

35  FORMAT (42X, 'VARIABLES  USED  THROUGHOUT',//) 

RF AD ( 1 , 9 ) MAPE 
RFADI  1,  9)  I RITE 
READ ( 1 ,9  ) IRATE 
RFADI 1,9)  JZAP 

RF AD! 1 , 10)  ( )G(M1  ) , Ml  = 1 , 121 

RE  AD  I 1 , 11)  A,  JJ,  BLAD1 , bLOOl,  BLAD2,  BL0D2 , TEST1,  TEST2, 

1TFST3 

READd,  12IRATE4,  RA1E3,  RATE2,  RAT  El 
SKIP  AU VANCES  PAPER  TU  TOP  OF  A NEW  PAGE 
CALL  SKIP 
WRITEI2,  3b ) 

WRITE ( 2 , 13)  A,  JJ,  bLADl,  BLOOl,  BLAD2,  BL002,  TEST  l,  TEST2, 
1TFST3,  RATt4,  RATES,  RAT E 2 , R ATE  1 
RF  AD  ( 1 , 14)  C,  C2,  (.3,  C4,  C5,  C6 

WRITFI2,  14)  C,  C2 « C3,  C4 , CS,  C6 
RF  AD ( 1,  lb ) (b<  IN) , IN  = 1,7) 

WR  ITE (2,  lb ) 1M(  IN) , IN  = 1,7) 

READd,  16)  ZL1M],  ZL  I M2 
WRITE (2,  17)  ZLIM1,  ZLIM2 
FRAME  CLEARS  IMAGES  FRUM  PLOTTING  SCOPE 
CALL  FRAME 

DRAW  GRID  AND  ANNO  TATE 
SENDS  PLOTTER  OFF  PAPER 

CALL  FPLUT  ( 1,  0.0,  29.0) 

SETS  NEW  SCALE 

CALL  SCALE !.b, .006,  -20. ,6)00.) 

SETS  CONSTANT  WHICH  INVERTS  LEGEND 


o o o 
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Y = 135.  * .0175533 
X = 180.  * .0175533 
C PLOTS  HE  AU  ING 

CALL  FCHAK! 12., 2 500 . , .50,  .50,  X) 
WRITE (8,15)  C,  L2,  C3,  C5,  C5,  C6 
C PUSIT1UNS  PIN  TO  DRAW  CKIO 

CALL  FPLOTtl,  -20.,  5500.) 

C Nb  XT  PU UK  STATEMENTS  OKAW  OK  I 0 

CALL  FGRIUI3,  -20.,  5500.,  100.,  28 ) 
CALL  bGKIO  (0,  -20.,  2700.,  1.,  50) 
CALL  FGK1U ( 1 , 20.,  2700.,  100.,  28) 
CALL  K;R1U(2,  20  . , 5500.,  1.,  50) 


C 

PLOTS  HEADING  '-20.0 

i 

CALL  FCHARI-17.9, 

2650. i 

, .3,  .35 

,X  1 

WRITE (8,  18)  8(1) 

, b(  7) 

c 

PLOTS  HEADING  >0  .()• 

CALL  FCHAK ( 1 .50 , 
WRITEIB,  19)  F ( 2 ) 

2650.  , 

.3,  .35, 

X) 

c 

PLOTS  HEADING  '20.0' 

CALL  FCHAR I 2 1 .8, 
WRITE (8,19)  8(3) 

2650., 

.3,  .35, 

X ) 

c 

PLOTS  HEADING  '3(00' 

CALL  FCHAK ( 2 3. 2 , 
WRITEIB,  19)  8(5) 

3025. , 

. 3 , . 35 , 

X) 

c 

PLOTS  HEADING  '5000' 

CALL  FCHARI23.2, 
WRITEIB, 19)  8(5) 

5025., 

.3,  .35, 

X ) 

c 

PLOTS  HEADING  '5000' 

CALL  FCHAK (23.2, 

5025.  , 

.3,  .35, 

X) 

WRITEIB,  19)  h l 5 ) 

SbT  LOOP  THAT  PROCESSES  EACH  HOLE 

00  3000  J = 1 , JJ 
TROSS  =0.0 
RTOT  = 0.0 
1 7 AP  = JZAP 
C READ  HEADING  EUR  EACH  HOLE 

KF  AO ( 1 , 20)  MM(J),  HOLE,  CLA01,  CL001 , OEPH,  AGEB 
WRITE (2,21)  HOLE 

WRITE (2,23)  CL  A 01 , C LUO  1 , OtPH,  AGEB 
MMM  = MM(J) 

RFAUd,  22)  (AGfc(M) ,01 ST(M) ,FL ITH(M) ,DENS(M> ,VEL(M>,  8=1, MMM) 

C FROM  TOTAI.  TU  GDENS  COMPUTES  THE  WEIGHTFD  AVERAGE  DENSITY  OF  THE  ENTIRE 
C SEDIMENT  PILE,  WHICH  IS  USED  TO  COMPUTE  THE  I N II  I AL  DEPTH  CORRECT  ION ( DFPHC  ) . 
TOTAL  = 0.0 
DO  1600  M = 2,  MM M 

TOT  = ( ( ( DEN  S ( M ) + OENSIM  - 1))*  .6)  A (l)IST(M)  - OIST(M  - 1))) 

TOTAL  = TOTAL  ♦ HIT 
1600  CONT INUE 

AOENS  = TOTAL  / l)lST(MMM) 

CALL  CKUSTIAGFH,  POEPH) 

W = ((3.3  - AD ENS  1/2.3)  * I>1ST(M8M) 

DOE  PH  = DEPH  W 
DFPHC  = UOLPH  - PDF  PH 

C IF  OF  PC  IS  NEGATIVE,  SEA  FLOOR  IS  TOO  SHAI LOW 
C IF  DFPHC  IS  POSITIVE,  SEA  FLOOR  IS  TOO  OFEP 
WRITEI2,  25)  PUEPH.DEPHC, AOENS 


5 . 


. -* 


-----  — ~ 
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IF  I I KIT  E)  lit),  110,109 

109  WRITE  12,  23) 

WR1TEI2,  2b)  (AGE(M),  0 1 ST ( M ) , F L IT H( M ) , DENS < M ) , Vtl(M),M=  1,  MMM ) 

110  IF  (NAPE  ) lit),  130, 165 
163  Wkl)tl2,  Zt  > 

WRITE 12,  26) 

IF(ZLI«1  - .0009)  120,  120,  130 

120  WRITE (2,29)  A06< 1 1 ,01 ST( 1 1 ,tt, A, A,A,FL ITH( 1 ),CLA01,CL0D1.0EPH 

C 

C KUTATF  POINTS 
C 

130  0(1  2000  H = 2,  MMM 

IF ( ( nGE ( M ) + .0009)  - 2LIM1)  2000,  140,  140 
1 AO  IF(AGEIM)  - t 2 L I M 2 + .00091)205,  205,  2000 
205  IF(AOktM)  - TEST3  ) 200,  210,  210 
200  BRPP1  = AGt(M)  * KATE4  s .0174333 

CALL  KUCCW  (bLAOl,  BLU01,  CLAU1,  CL001,  BBPP1,  CLAD,  CLOD) 

GO  10  600 

210  IFIAGfc(M)  - T E S f 2 ) 220,  230,  230 

220  HBPP1  = |TEST3*KATE4  + (AOttM)  - T t SI  3 ) *KATE3 1 * .0174533 
CALL  KUCCW ( BLAOl , BLOUl , CLAOl.CLOol,  BBPP1,  CLAD,  CLOD) 

GO  TO  600 

230  IF ( AGE ( M > - TEST1)  240,  250,  250 

240  B8PP1  = MEST3*kATE4  + ITES12  - TEST3 ) 4K AT  E3  + ( AGE  t M ) - TEST2 ) # 
IK AT 1 2 ) * .0174533 

CALL  KUCCw ( BLAOl , BLUDl , CLA01,  CLUD1 , BBPP1,  CLAU,  CLOD) 

GO  TU  600 

250  BBPP1 = ( TE  ST  3 #R A TE4  +ITEST2-  TEST3)»KATE3  + ( T EST  1 - TEST  2 )* 
1KATE2)  * .0174533 

CALL  ROCCW ( BLAU  1 , BLUDl,  CLAU1 , CL0D1,  BBPP1,  CLAU,  CLOD) 

CAO  = CLAO 
COO  = CLOU 

CALL  KUCCw  t BLAOl , BLOUl , BLA02,  BLUD2.  BBPP1,  BA2  , BU2 ) 

BBPP2  = IIAGE(M)  - TEST  1 ) *KATE  1 ) » .0174333 
CALL  KUCCw ( BA2 , BO 2 , CAD,  COO,  BHPP2,  CLAU,  CLOD) 

600  AGET  = AGEB  - AGE (M  ) 

XT  H 1C  = DIST(MMM)  - 0 1ST  t M ) 

C DETERMINE  DEPTH  UE  POINTS 

CALL  CRUST  t AGET , PUEPh) 

Ml  = M 
TOTA  = 0.0 

TTOT  = 0.0 

C DETERMINE  AVERAGE  DENSITY 

Of)  700  J1  = Ml  , MMM 

TTOT  = {(0LNS1J1)  +DENSIJ1  - 1 ) ) *. 5 ) * ( 0 1 ST ( J 1 ) - 01STIJ1-1)) 

TOTA  = TOTA  + TTUf 
700  CONTINUE 

BOEMS  = TOTA  /(OIST(MMM)  - UlSTtM  - D) 

W = ((3.3  - BOENS)  22.3)  * XT  NIC 
PDPPpi  = PUt PU  - w 
PDF  PH  = PUEPH  ♦ DEPHC 

KATE  = (OIST(M)  - OISKK  - 111  / (AGE(M)  - AGE  I M - 1)) 

RTUT  = kTUT  + kATL*(l)IST(M|  - D)ST(M-l)) 

ARATE  = KTUT/Ul ST (M ) 

604  IF  (CL'JO)  KOb,  605,  B06 

605  ACLUO  = CLUO 
GO  To  HOB 

606  IF (CLOU  - 160.)  R05 , b05,  607 
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807  ACIUU  = CLUO  - 360. 

C OfcTtRMIWt  CORRECTED  lNltkVAL  TnlCKNESStS 

808  lMl)tNS(M)  - .001)  812,  612,  6lA 
812  AVUEN  = (UtMS(MI  + DbN$(M-l)|  * .5 

GO  TO  819 

8 1 A IF(UENSIN-l)  - .001)  816,  6 1 6 , 812 
816  AVDEN  = OtNS(M) 

C DETERMINE  RATtS  OF  StU  I ME  NT  A T 1 UN 

819  PS1  = (AVUEN  - 2.6)  / (1.008  - 2.6) 

HO  »((l.  - PSI)/(I.  - .76))*  (f)JST(MI-OIST(M-lJI 
URATt  = HU  / ( AGt  ( M ) -AGE(M-l)) 

CKAU  = URaTE  * (UISflM)  -UlST(h-l)) 

TRUSS  = TbUSS  + CKAlfc 
CE  RAT  = T8USS/U1 ST (M ) 

IF (NAPE)  909,  909,  900 

900  WRITt (2,30) AGE ( M ) , 01 S T ( M ) , K AT e, hn , Uk A T E ,C EK6 T , F L I TH ( M ) , CL AO, AC  LUO, 
IPOtPH 

C PLOT  RA1ES  UR  CARBONATE  PERCENTAGES 
909  IF(IkATE)  916,  916,  911 

911  CALL  ALPHA(URATE, J1 , J2,J3, JA) 

CALL  ECHAklCLAO.PDtPH,  .06,  .07,  Y) 

IF(IZAP)  912,  913,  91a 

912  WRITE ( 8 , 31)  I G ( J 1 ) , 1 G ( J2 ) , I G ( J3 ) , I G ( 1 1 ) , I G ( JA ) , HUL E 
IZAP  = 0 

GO  TU  2000 

913  WRITt(8,  32)  I G ( J1 ) , IG< J2) , IGI Jr ) , I G ( 11  ) , lG(  JA) 

GO  TU  2000 

91 A WRITEI8,  31)  1GIJ1  ), I G I J 2 ) , 1 G ( J 2 ) , IG ( 11 ), IG(JA) , HOLE 
GO  TU  2000 

916  CALL  FCHARICLAO,  POEPh,  .06,  .07,  Y) 

IF(IZAP)  921,  922,  92A 

921  WRITE  ( 8 , 33)  FLITH(N),  HULL- 
IZAP  = 0 

GO  TU  2000 

922  WRITt (8,  3A  > FLITH(H) 

GO  TO  2000 

92 A WRITEIH,  33)  FLITHIM,  HOLt 
2000  CONT INUh 
CALL  SKIP 
3000  CONI lNUE 
CALL  fcXIT 
END 
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PROGRAM  ESP 

Software  and  Hardware 
Fortran  IV,  Burroughs  6700 

Preparation 

Each  of  the  grids  plotted  by  program  ROPH  must  be  contoured, 
then  encoded  into  arrays  for  use  by  Program  ESP.  The  data  is  entered 
in  two  (40  x 40  x 1)  arrays  as  variable  ALITH  (percent  carbonate)  and 
ATHIC  (adjusted  rate  of  sedimentation).  The  numbers  entered  in  each 
array  record  the  value  at  the  midpoint  of  the  latitude  and  of  the 
depth-of-water  intervals,  Figure  2.1. 

Purpose 

Program  ESP  predicts  thickness  and  kind  of  sediment  that  may  be 
expected  at  arbitrarily  specified  points  on  the  east-central  portion 
of  the  Pacific  plate.  The  prediction  is  accomplished  by  incremental 
counterclockwise  rotation  of  the  point  back  through  time,  using  the 
latitude  and  depth  of  water  as  the  independent  variables.  Rate  of  sedi- 
mentation and  carbonate  percent,  the  dependent  variables,  are  entered 
in  separate  (40  x 40  x 1 ) arrays,  with  latitudes  and  depth  of  water  as 
the  coordinates.  At  each  iteration,  these  grids  are  entered,  and  the 
extracted  values  used  to  predict  thickness  of  sediment  and  percent  of  car- 
bonate of  the  sedimentary  column. 

Variables 

R0PHT  = an  array  that  contains  interval  thicknesses  corrected  for 
compaction  determined  by  program  ROPH  for  each  DSDP  site; 
the  interval  (XINC)  for  ROPH  and  ESP  must  agree. 
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ROPHC 


CARRY (I ) 


CCARY(I) 


IMONY 


IAGE5 


TESTA 


TESTS 


an  array  that  contains  carbonate  percentages  entered  in  pro- 
gram ROPH  for  each  DSDP  site;  the  interval  (XINC)  for  ROPH 
and  ESP  must  agree. 

an  array  that  contains  the  digitized  curve  describing  varia- 
tion of  sedimentation  rate  with  time  (Fig.  11);  the  dimension 
is  17;  values  are  entered  at  2.5  million-year  intervals, 
starting  with  0 my. 

an  array  that  contains  the  digitized  curve  describing  the 
variation  of  carbonate  percentages  with  time  (Fig.  10);  the 
dimension  is  17;  values  are  entered  at  2.5  million-year  inter- 
vals, starting  with  0 my. 

An  option  for  supressing  the  printing  of  the  carbonate  array 
(ALITH)  and  the  rate  of  sedimentation  array  (ATHIC);  if  set 
to  1,  arrays  do  not  print;  if  set  to  0,  printing  occurs. 

An  option  for  supressing  printing  of  iterative  solutions; 

set  equal  1 or  greater,  solutions  will  be  printed  every 
IAGE5  years. 

Maximum  age  (mybp)  attained  at  the  base  of  the  modeled  sedi- 
mentary sequence.  If,  for  instance,  the  basement  age  is  110  my, 
but  solutions  for  the  past  37  my  are  desired,  TESTA  = 37.0 
will  stop  execution  at  this  time  for  all  sites. 

Duration  (my)  for  which  IB  = 1.  In  its  present  form,  the 
program  uses  1 set  of  2 latitude-depth-of-water  arrays.  With 
alteration  of  the  Dimension,  Read,  and  Write  statements  for 
arrays  ATHIC  and  ALITH,  a second  set  of  2 arrays  (IB  = 2) 
could  be  entered  to  reflect  a different  sedimentation  pattern. 


LL 

XINC 
BLAD1 

BLAD2 
BL0D1 

BLOD2 
TEST1 

TEST2 

TEST3 

RATE1 

RATE2 

RATE3 

RATE4 

ALITH( ID, 
IC, IB)  = 

ATHIC( ID, 
IC, IB)  = 

CLAD1 


The  number  of  points  for  which  a stratigraphic  section 
is  required. 

The  length  (my)  of  the  incremented  rotations. 

Latitude  of  first  pole  of  rotation  (Hawaiian)  entered  as 
decimal  fraction  of  degrees.  Use  north  projection  only  with 
positive  sign. 

Latitude  of  second  pole  (Emperor).  Same  convention  as  BLAD1. 
Longitude  of  first  pole,  entered  as  decimal  fraction  of  degrees. 
Use  positive  values  for  east  longitude,  negative  values  for 
west  longitude. 

Longitude  of  second  pole.  Same  convention  as  BL0D1. 

Age  (mybp)  for  change  from  Hawaiian  to  Emperor  pole.  For 
ages  greater  than  TEST1,  RATE1  applies. 

Age  (mybp)  for  change  from  RATE3  to  RATE2. 

Age  (mybp)  for  change  from  RATE4  to  RATE3. 

Rate  of  rotation  (my/deg)  for  Emperor  pole. 

Rate  of  rotation  (my/deg)  from  age  TEST1  to  TEST2  years. 

Rate  of  rotation  (my/deg)  from  age  TEST2  to  TEST3  years. 

Rate  of  rotation  (my/deg)  from  age  0 to  TEST3  years. 

Percent  of  carbonate  entered  in  latitude-depth-of-water  grid. 

Adjusted  sedimentation  rate  entered  in  latitude-depth-of- 
water  grid  . 

Latitude  in  decimal  fraction  of  degrees  for  point  selected 
for  prediction  of  stratigraphic  section;  positive  for  north 
latitude,  negative  for  south. 
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CL0D1  = Longitude  in  decimal  fraction  of  degrees  for  point  selected 
for  prediction  of  stratigraphic  section;  same  convention  as 
BL0D1 . 

AGEB  = Age  (mybp)  of  basement  at  latitude  CLAD1,  longitude  CL0D1. 

DEPH  = Depth  of  water  (m)  at  point  CLAD1,  CL0D1. 

NSITE  = Starting  points  may  be  given  designators  of  5 integer  length. 

Out 

CLAD  = Latitude  (convention  same  as  CLAD1)  of  point  after  incremental 
rotation. 

ACLOD  = Longitude  of  point  after  incremental  rotation  (convention 
same  as  CL0D1. 

AGET  = Age  (mybp)  of  points  after  incremental  rotation. 

PDEPH  = Depth  of  sea  floor  (m)  at  point  CLAD,  ACLOD. 

XLITH  = Percentage  of  carbonate  at  point  CLAD,  ACLOD. 

CAC03  = Cumulative  weighted  mean  carbonate  percent. 

XTHIC  = Uncompacted  thickness  (m)  predicted  for  each  increment. 

CINT  = Thickness  (m)  for  each  increment  after  correction  for 

compaction. 

COMPA  = Cumulative  adjusted  thickness  of  sediment. 

ID  = Subscript  for  depth  in  variable  ALITH  and  ATHIC. 

IC  = Subscript  for  latitude  in  variable  ALITH  and  ATHIC. 

IB  = Subscript  signifying  grid  number.  Ordinarily  IB  will  be  1, 

unless  program  is  altered. 

VEL  = Velocity  of  sediment,  determined  from  thickness  of  sediment. 

ROPHT  = Same  as  input. 


ROPHA 

TDIFF 

PDIFF 


ROPHC 

ROPHB 

CDIFF 

PCDIF 


= Three-point  average  of  ROPHT. 

= Difference  in  interval  thickness  (uncompacted  or  corrected) 
between  ESP  and  ROPH. 

= Percent  difference  in  thickness  between  ESP  and  ROPH. 

This  also  appears  in  output  table  as  "Error  in  uncompacted 
interval  thickness." 

= Same  as  input. 

= Three-point  average  of  ROPHC. 

= Difference  in  interval  carbonate  between  ESP  and  ROPH. 

= Percent  difference  in  interval  carbonate  between  ESP  and 
ROPH.  This  also  appears  in  output  table  as  "Error  in  inter- 
val carbonate." 


Subroutines 


ROCCW,  COUNT,  CRUST,  TABLE,  CPRIM,  XLOOK,  and  functions  FTAN  and  XTAN 
Features 

Two  poles  of  rotation  may  be  used.  Uses  one  set  of  2 arrays  for 
rate  of  sedimentation  and  percent  of  carbonate,  but  can  be  altered  to 
include  a second  set.  Includes  a routine  for  interpolating  values  from 
arrays  ALITH  and  ATHIC. 

Input  of  Data 
Image 


Number 

Variables 

Format 

1 

BLAD1 , 

BL0D1 , 

BLAD2 , BL0D2 

4 F10.2 

2 

TESTA, 

TESTB, 

XINC 

3 F10.2 

3 

TEST1 , 
RATE4 

TEST2 , 

TEST3,  RATE1 , RATE2,  RATE3, 

3 F10.2,  4 F10.1 

88 

Image 

Number 

Variables 

4 

LL 

5 

IM0NY,  IAGES 

6 

CARRY 

7 

CCARY 

8 

ALITH 

9 

ATHIC 

10 

CLAD1 , CL0D1 , AGEB,  DEPH,  NSITE 

11 

R0PHT 

12 

R0PHC 

Format 

110 

212 

20F4.2 
20F4.2 
1 0F8 . 1 

10(F3.1 , IX) 

3F10.2,  20X,  FI 0.2,  IX,  15 
20(F3. 1 , IX) 

26(F2.0,  IX) 


Method 

Step  1.  Read  and  print  images  1 through  9 of  the  Input-of-Data  list. 

The  cards  (images  1 through  9)  read  in  here  are  common  to  all  sites. 

In  addition  to  reading  and  printing,  this  step  includes  an  interpolati ve 
routine  that  reduces  the  effort  in  card  punching  required  to  translate 
sedimentation  rates  (determined  from  the  curves  of  Figure  8)  to  the  ele- 
ments of  array  ATHIC.  Rather  than  punch  a value  for  each  element  of  ATHIC, 
a 40  x 40  x 1 array,  the  user  need  only  enter  selected  elements  along  each 
row  of  the  array,  and  the  routine  will  interpolate  between  entries.  The 
sample  output  (Appendix  II)  shows  an  example.  The  values  punched  on  cards 
are  printed  under  the  heading  "Array  of  sedimentation  rates  before  interpo- 
lation." The  array  used  in  the  computations  after  interpolation  appears 
under  the  heading  "Array  of  sedimentation  rates."  Interpolation  is  linear 
and  proceeds  horizontally  across  the  rows  of  the  array.  As  an  example, 
consider  Row  1 of  the  sample  output.  A sedimentation  rate  of  7.0  m/my 
appears  in  Column  1,  zeros  in  Columns  2 through  7 and  a rate  of  10.0 
m/my  in  Column  8.  The  interpolati ve  routine  calculates  rates  for  Columns 
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ATHIC  (ID, 1C, IB) 

22  23  24 

-2°  -3°  -4° 


4000 


18- 


4100 


19- 


4200' 


20- 


4300- 

L_ 


15.0 

13.0 

10.0 

13.0 

12.2 

8.5 

11.0 

10.5 

7.5 

Subscript  IC 
Latitude  (deg.) 


Sedimentation  rate 
(m/m.y.) 


Depth  (m) 
Subscript  ID 


Figure  2.1.  A portion  of  array  ATHIC  (ID,  IC,  IB)  showing  position 
of  values  for  rate  of  sedimentation  relative  to  depth  and  latitude 
boundaries. 
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Figure  2.2.  An  enlarged  portion  of  Figure  2.1,  showing  interpolating 
procedure  --  interpolate  between  13.0  and  10.0  to  get  point  b;  inter- 
polate between  12.2  and  8.5  to  get  point  a;  interpolate  between  a and 
b to  get  rate  of  sedimentation,  point 


2 through  7 by  assuming  a linear  increase  in  sedimentation  rate  between 
7.0  and  10.0  m/my.  A zero  in  Column  1 suppresses  interpolation  for  all 
leading  zeros  along  that  row.  A zero  in  Column  40  suppresses  interpolation 
between  that  column  and  the  first  non-zero  entry  to  the  left  of  Column  40. 

Step  2.  Rotate  the  second  pole  (Emperor)  around  the  first  pole  (Hawaiian). 

Because  the  age  of  all  points  that  are  to  be  rotated  may  be  too  young, 
rotation  around  the  Emperor  pole  may  not  occur;  nevertheless,  the  possi- 
bility is  anticipated  here,  and  the  Emperor  pole  is  shifted  in  prepara- 
tion. Consider  a point  that  represents  a site  which  is  50  million  years 
old,  and  assume  that  the  Hawaiian,  pole  describes  rotation  of  the  Pacific 
plate  for  the  first  40  million  years  and  the  Emperor  pole  describes  it 
for  the  next  10  million  years.  The  Doint  must  first  be  rotated  around 
the  Hawaiian  and  then  around  the  Emperor  pole.  The  Emperor  pole,  as  its 
coordinates  are  generally  given,  describes  the  angular  distance  of  the  point 
from  the  Emperor  pole,  assuming  that  no  rotation  has  occurred  around 
the  Hawaiian  pole.  Once  rotation  of  the  point  has  occurred  around  the 
Hawaiian  pole,  the  Emperor  pole  no  longer  correctly  describes  the  rela- 
tionship, because  the  angular  distance  between  the  Emperor  pole  and  point 
are  changed,  unless  of  course,  the  Emperor  pole  has  been  rotated  the  same 
distance  around  the  Hawaiian  pole  as  the  point.  This  step  rotates  the 
Emperor  pole.  In  this  example,  the  angular  distance  will  be  age  (40  my) 
times  the  rate  of  rotation  (deg/my)  that  is  assumed  for  the  Hawaiian  pole. 

Step  3.  Read  and  print  image  10. 

This  card  specifies  the  latitude  and  longitude  of  the  starting  point, 
the  age  of  the  basement  of  the  point,  depth  of  water,  and  an  optional  5- 
digit  number  that  is  used  as  a designator  for  the  hole. 
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Step  4.  Determine  crustal  correction. 

The  subsidence  curve  predicts  the  depth  to  the  crust.  The  start- 
ing point  has  as  its  depth  the  depth  of  water  at  the  sediment-water 
interface.  The  difference  between  these  depths  comprises  the  crustal 
correction,  incorporating  two  components:  the  thickness  of  sediment 

overlying  the  crust  and  the  distance  the  crust  departs  from  its  ideal 
level  predicted  by  the  subsidence  curve.  This  correction  is  applied  at 
Step  9. 

Step  5.  Determine  number  of  iterations. 

The  number  of  iterations  (Kl)  is  determined  from  either  age  TESTA 
or  AGEB,  whichever  is  the  smaller,  and  dividing  by  the  length  of  time 
specified  as  the  increment  of  rotation  (XINC),  adding  .0009  to  prevent 
an  error  in  rounding  off,  and  converting  the  result  to  the  integer  Kl. 

Step  6.  Read  in  observed  thicknesses  (ROPHT)  and  carbonate  per- 
centages (ROPHC),  images  11  and  12. 

Step  7.  Determine  the  age  after  each  iteration. 

Age  is  determined  by  setting  a counter  AGET  equal  to  zero,  and  at 
each  iteration  adding  the  increment  of  rotation  (XINC)  to  the  previous 
resul t. 

Step  8.  Determine  angular  distance  of  rotation. 

Angular  rotation  is  determined  by  multiplying  the  rate  (deg/my) 
times  the  age  of  the  point  (my)  times  the  constant  0.0174533,  which  con- 
verts degrees  to  radians.  Rates  (RATE4,  RATE3,  and  RATE2)  are  used 
for  the  first  pole,  and  (RATE1)  for  the  second,  Figure  2.3.  TEST3,  TEST2, 
and  TEST1,  Figure  2.3,  determine  which  rates  will  apply.  If  the  age  of  a 
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Figure  2.3.  Variables  TEST3,  TEST2,  and  TEST1  set  off  interval  of 
time  during  which  the  rates  of  rotation  RATE4,  RATE3,  and  RATE2, 
determine  the  rotation  around  the  first  pole  (Hawaiian).  RATE1 
determines  rotation  around  the  second  pole  (Emperor). 


point  AGET  is  greater  than  TEST3,  but  less  than  TEST2,  then  RATE4  applies 
for  the  first  TEST3  years,  and  RATE3  applies  for  the  next  AGET  minus  TEST3 
years.  Solutions  for  points  older  than  TEST2  years  are  determined  by  simi- 
larly reasoned  processes. 

Step  9.  Calculate  initial  depth  of  water. 

At  this  point  in  the  calculations,  the  depth  of  water  is  that  which 
prevailed  at  the  end  of  the  preceding  increment  (previous  iteration),  but 
the  depth  required  is  that  at  the  end  of  the  current  increment.  If  the 
preceding  depth  is  used  to  enter  the  rate  of  sedimentation  array  (ATHIC), 
the  rate  and  consequently  the  thickness  determined  will  be  in  error, 
because  the  depth  of  water  is  too  shallow.  To  overcome  this  problem,  the 
depth  from  the  previous  iteration  is  adjusted  by  removing  the  thickness 
of  sediment  determined  from  the  previous  iteration  according  to  the  equa- 
tion: 

2.9.1  G = C + T + P 

where  G = depth  of  water  at  end  of  last  iteration 
C = crustal  correction  (Step  4) 

T = total  thickness  of  sediment  removed 
P = depth  of  crust  according  to  subsidence  curve 
(subroutine  CRUST) 

The  depth  determined  from  equation  2.9.1  is  used  to  enter  array 
ATHIC  once  again,  thereby  determining  a new  thickness  and  a new  depth. 

The  process  is  repeated  twice  more,  using  the  adjusted  depth  each  time 
to  determine  a new  depth  and  thickness.  The  results  at  the  end  of  the 
third  iteration  are  used  as  the  solution  for  depth  and  thickness  of  sedi- 
ment. The  three  iterations  begin  at  Step  10. 
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Step  10.  Determine  subscripts  for  variables  ATHIC  and  ALITH. 

Carbonate  percentages  are  stored  in  array  ALITH,  sedimentation 
rates  in  array  ATHIC.  Subscripts  for  these  3-dimensional  arrays  are 
determined  from  the  latitude  and  depth  of  water  by  subroutine  TABLE 
(Appendix  III). 

Step  11  . Determine  the  sedimentation  rate. 

The  sedimentation  rate  is  determined  by  entering  array  ATHIC  using 
subscripts  determined  by  subroutine  TABLE.  Subroutine  XL00K  (Appendix  III) 
interpolates  among  variables  of  array  ATHIC. 

Step  12.  Determine  the  percentage  of  carbonate. 

Carbonate  is  determined  as  for  sedimentation  rate,  Step  11  , using 
array  ALITH. 

Step  13.  Correct  thickness  and  carbonate  percent  for  variation  with  time. 

The  arrays  ALITH  and  ATHIC  describe  an  average  regime  of  sedi- 
mentation for  the  past  42  my.  Corrections  are  made  for  temporal  changes 
in  rate  of  sedimentation  and  percentage  of  carbonate  at  this  step  by  appli- 
cation of  the  1-dimensional  arrays  CARRY  (rate)  and  CCARY  (carbonate). 

Array  CARRY  is  constructed  by  digitizing  the  curve  that  relates  deviation 
from  the  average  sedimentation  rate  to  time  (Fig.  11).  The  values  entered 
in  array  CARRY  indicate  the  percentage  that  the  sedimentation  rate  (ATHIC) 
must  be  increased  or  decreased  for  any  particular  time. 

An  example  of  the  preparation  and  use  of  CARRY  follows:  from  Figure 

11  percentages  are  scaled  at  1 -0  my  intervals;  the  sign  of  the  percentage 
is  reversed  and  the  value  entered  in  array  CARRY;  the  sedimentation  rate 
is  determined  from  array  ATHIC;  array  CARRY  is  entered  using  the  appropriate 
age;  if  the  age  falls  between  digitized  elements  of  CARRY,  interpolation 
occurs;  the  percentage  determined  from  CARRY  is  multiplied  times  the  rate 
determined  from  ATHIC;  the  result  is  the  corrected  rate. 
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Array  CCARY,  which  determines  carbonate  percentages,  is  prepared 
from  Figure  10  and  used  in  conjunction  with  array  ALITH.  All  other  details 
of  CCARY  are  as  for  ALITH  and  CARRY. 


Step  14.  Determine  the  depth  of  water. 

This  step  contains  the  solution  for  depth  that  began  in  Step  10. 

The  depth  is  calculated  according  to 

2.14.1  P - P,  - H0  - (H0(l-n)(Dg-Dw))/Db) 

where  P = depth  of  water  after  removal  of  sediment  layer 

P.j=  depth  of  water  prior  to  removal  of  sediment  layer 
Hq=  thickness  of  sediment  removed  (determined  from  array 
ATHIC) 

n = porosity  (in  this  program  assumed  to  76  percent) 

Dg=  density  of  grains 

D, ,=  density  of  water 
w 

Db=  density  of  basement. 

If  the  number  of  iterations  for  depth  of  water  is  less  than  4,  con- 
trol returns  to  Step  10,  using  the  depth  of  water  P as  the  initial  depth 


of  water. 


Step  15.  Determine  the  cumulative  thickness  of  sediment  and 
apply  compaction  correction. 

The  correction  for  compaction  is  derived  from  the  relationship 
between  the  compacted  thickness  determined  from  DSDP  cores  and  the 
corrected  thicknesses  adjusted  to  76  percent  porosity,  Figure  2.4. 

The  ordinate  (Fig.  2.4)  is  the  ratio  HQ/Hd,  where  HQ  is  the 
observed  DSDP  thickness  and  H,  is  the  DSDP  thickness  adjusted  to 

a 

76  percent  porosity.  The  abscissa  (Fig.  2.4)  is  H . 

a 

The  points  plotted  in  Figure  2.4  were  determined  at  5 million 
year  intervals  from  the  data  of  all  23  DSDP  sites.  To  be  plotted  on 
the  figure,  the  corrected  thickness  adjusted  to  76  percent  porosity 
had  to  be  greater  than  20  m --  experience  showed  that  the  co- 
ordinates of  corrected  thicknesses  less  than  20  m had  extremely 
high  variation.  An  arbitrary  value  of  1.5  was  used  for  the  zero 
intercept  on  the  ordinate  (Fig.  2.4)  and  the  curve  was  fitted  by 
eye. 

The  curve  (Fig.  2.4)  was  digitized  at  25  m intervals  and 
entered  in  program  ESP  as  variable  CORTT.  Interpolation  occurs 
between  members  of  the  array. 

The  correction  for  compaction  is  made  using  the  relationship: 

2.15.1  Hc  = Ht  * H0/Ha 

where  Hc  = compacted  thickness  of  sediment 

= cumulative  thickness  of  sediment  predicted  from 


2 . w c 

+->  CD  • 

a>  a.  c 4-> 

rQ  QJ'r 

4->  00  U CO  • < 

OS-  T3 

CO  CD  Q-  CD  - 

•r-  03  0.0  4->  ■ 
OO  -M 

<D  C LO  Q O V 

4- >  •!-  r— 

03  ro  Cl 
C “O  4-  C\J  r 

•i-  a»  o a»  r 

T3  > i—  S-  « 

5-  S-  >>r-  03 

O CD  4->  ro 

CO  «r-  E 

a>  -O  to  E 

.c  o o o o 

h-  S-  S-  C\j 

to  O 4- 

to  Q.  “O  • 

• CD  to  CD  1 

C C rO  «—  T3 

O -X  <T3  CD 

•r-  U O > CD. 

4- >  *r-  4->  l-  CJ  ■ 

(J  -C  CD  X 

CD  4->  “O  4->  CD  *■ 

l-  CD  C 

5-  CD  4-1  -r-  « — « 

O -C  O O* 


C CO  S-  CD 

O *r-  O 10  l 

•r-  CJ  CO 

4- >  O C CD 

UI  (O  o c. 

03  CO  *r-  - 

Q_  CD  CD  < — U 

E S-  C r—  -r- 

o CD  *r-  -C  • 

o -c  o E +->  • 

s-  r cn-a  . 

O 034->  CD 

4—  m -»->  > 

CD  03  S- 

“D  CO  E <D  - 

CD  -r-  03  TD  to 

to  (O  CD  -Q  ■ 

3 03  CO- 

(OCD-r- 

CD  tO  _C  E <D  r 
> *r-  4->  i-  CO 

5-  CJ  CD  O 

3 CO  CO  4->  JC 
(J  JD  'r  CD  }6  - 
03  “O 

• CO  CO  • 

(D  i a;  -p 

• -C  s-  c 

C\J  D D 03  -r- 

C O C 
CD  *»  03  CO  CL  r 
S-  03  4-> 

31  CD  C >) 
CT>\  S-  *r-  •—  r 
O O O C r 
U.  I UQ-O 


■ 


98 


array  ATHIC  (Porosity  76  percent  and  uncorrected  for  com- 
paction; the  cumulative  thickness  Ht  is  the  abscissa  of 
Fig.  2.4) 

HQ/Ha  = the  ratio  determined  by  interpolation  in  array  CORTT, 
given  cumulative  thickness  Ht  (the  ratio  HQ/Ha  is  the 
ordinate  of  Figure  2.4). 

Step  16.  Determine  velocity  of  sediment. 

The  velocity  of  sediment  is  computed  from  a function  relating 
velocity  of  sediment  to  thickness  of  overburden.  Figure  2.5.  The 
function  is  determined  empirically  using  solutions  from  wide-angle 
reflection  records  and  by  computing  velocities  from  normal -incicence 
seismic  records,  where  acoustic  horizons  of  the  records  can  be 
correlated  with  lithologic  horizons  of  the  DSDP  cores  (basement 
or  top  of  Eocene  chert).  The  curve  is  fitted  by  eye. 

Step  17.  Optional  comparison  of  observed  and  predicted  interval 
carbonate  percentages  and  thicknesses. 

As  a test  of  the  effectiveness  of  program  ESP,  the  per- 
centages of  carbonate  determined  from  the  cores  of  the  various  DSDP 
sites  (values  used  by  program  ROPH)  can  be  compared  with  the  values 
predicted  by  ESP  (Fig.  10).  Similar  comparisons  can  be  made  be- 
tween the  calculated  and  predicted  thicknesses  (Fig.  11).  At  this 
step,  which  the  user  employs  at  his  option,  the  differences  between 
predicted  and  observed  values  of  carbonate  percentage  and  thickness 
are  calculated. 


The  interval  thickness  (variable  URATE,  program  ROPH)  is  entered 
in  ESP  as  variable  ROPHT , the  interval  carbonate  percentage  (variable 
FLITH,  program  ROPH),  as  variable  ROPHC.  A running  3-point  average 
(ROPHA)  is  made  of  ROPHT  and  a running  3-point  average  (ROPHB)  of 
ROPHC.  These  averages  are  compared  with  the  predicted  values  and 
the  results  printed: 

CDIFF  = XLITH  - ROPHB  (interval  carbonate  difference) 

TDIFF  = XTHIC  - ROPHA  (interval  thickness  difference) 
and 

PCDIF  = CDIFF/XLITH  * 100.  (percent  carbonate  difference) 

PDIFF  = TDIFF/XTHIC  * 100.  (percent  thickness  difference) 

If  the  results  of  this  step  are  of  no  interest  to  the  user, 
he  or  she  should  insert  5 blank  cards  for  variables  ROPHT  and 
ROPHC  (card  images  11  ar.d  12). 
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ESP  LISTING 


DIHENS ION  CO  BTT  ( 81 ) 

DIMENSION  VE LA  { 1 9) 

DIMENSION  CARRY(4l),R0PHT(45)  , NS  (25)  , PD  <25,  45) 

DIMENSION  ROPIIC  (45)  , PC  (25,  45)  ,CCAR  Y (41) 

DIMENSION  ATHIC  (40,41,1)  , ALITH  (40 , 4U , 1) 

INPUT  POLE  DATA,  AGE  DATA,  AND  NUMBER  OP  STARTING  POINTS 

READ (1,22) BLAD1 , BLOD1 , BL AD2 , BLOD2, T EST A.TESTB, XINC, TEST  1 , TEST2, 

1 TEST3.RATE1, RAT E2,R ATE3,RATE4,LL 
22  PORMAT(4F10.2/3F10.2/3F10.2,4F10.  1/110) 

WRITE (2,  4)  BL  AD  1,  BLODl,BLAD2, BLOD2, TESTA , TESTB , XI NC.TEST1 ,TEST2 , 

1 TEST3,RATE1 ,E ATE2 , RATE3 , EAT E4, LL 
4 FORMAT (8X,'E  S P I N P U T'// 

1 3X, ’LATITUDE  OF  HAWAIIAN  POLE  (BLAD1)  IS',P9.2/ 

1 3X, ‘LONGITUDE  OF  HAWAIIAN  POLE  (BLOD1)  IS‘,F9.2/ 

1 3X, 'LATITUDE  OF  EMPEROR  POLE  (BLAD2)  IS‘,F9.2/ 

1 3X, ‘LONGITUDE  OF  EMPEROR  POLE  (BLOD2)  IS‘,F9.2// 

1 3X, 'TESTA  IS' ,F6. 1 , • MILLION  YEARS'/ 

1 3X, 'GRID  IB= 1 APPLIES  FOR  ',F4.1,'  MILLION  YEARS'/ 

1 3X, ' AGE  INTERVAL  IS',F5.2,*  MILLION  YEARS'// 

1 3X, ' AGE  OF  TEST  1 IS',F6.2,'  MILLION  YEARS'/ 

1 3X, ' AGE  OF  TEST 2 IS',F6.2,'  MILLION  YEARS'/ 

1 3X, 'AGE  OF  TEST 3 IS',F6.2,'  MILLION  YEARS'// 

1 3 X, ' RATE  1 FOR  AGES  OLDER  THAH  TEST1  IS',F5.2,'  DEG/MY'/ 

1 3X, 1 R ATE2  FROM  AGES  TEST2  TO  TEST1  IS',P5.2,'  DEG/MY'/ 

1 3X, ' R ATE3  FROM  AGES  TEST3  TO  TEST2  IS‘,F5.2,'  DEG/MY',/ 

1 3X, • RATE4  FROM  AGES  0 TO  TEST3  IS',F5.2,'  DEG/MY'// 

1 3X, 'NUMBER  OF  STARTING  POINTS  IS', 14/) 

..PRINTING  OPTIONS 

IMONY  AND  IAGC5  ARE  FLAGS. 

IF  IHONY  = 1 OR  MORE,  THEN  CARBONATE  AND  SED  RATE  ARRAYS  WILL  NOT 
BE  PRINTED  OUT.  IF  I AGE  =1  OR  MORE,  THEN  FORWARD  TRACKS  WILL  BE 
PRINTED  OUT  EVERY  IAGE5  MILLION  YEARS. 

READ  ( 1 , 11)  IMONY,  I AGES 
11  FORMAT  (212) 

INPUT  AGE  DEPENDENCE,  ARRAYS  CARRY,  CCARY. 

READ(1,  3)  (CARRY  (I),  I =1,  41) 

HRITE  (2,  74)  (CARRY  (I)  , 1*1,  41) 

74  FORMAT  (3 X , 'AGE  DEPENDENCE' ,/,5X, 'SEDIMENTATION  R ATE' , /, ( 7X , 20 ( IX , F 

11.2))) 

READ  (1,3)  (CCARY  (I),  1=1,  41) 

3 FORMAT  (20F4.  2) 

WRITE  (2 , 75)  (CCARY(I),  1*1,  41) 

75  FORMAT (/,5X, 'CARBONATE  PERCENT' , /,  (7 X, 20 ( 1 X , F4. 2) ) ) 

C... INPUT  THICKNESS  CORRECTION  ARRAY (VALUES  AT  25  M.  INTERVALS). 

READ  (1,9)  (CORTT(I),  1*1,  81) 

9 FORMAT  (20F4. 3) 

WRIT E (2,  76)  (CORTT(I),  1=1,  81) 
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76  FORMAT (/,5X, 'COMPACTION  ABB A Y • ,/, (7X, 20 ( 1 X, F4. 3) ) ) 
BEAD  (1»  8)  (VELA  (I),  1*1,  19) 

WRITE(2,  12)  ( VELA  (I)  , 1*1,  19) 

12  FOB HAT (/,5X,' VELOCITY  ARRAY',/,  (7X,19(1X,P4. 2))) 

8 FORMAT  ( 19  ( IX  , F3.  2)  ) 

INPUT  CARBONATE  AND  SEDIMENTATION  RATE  ARRAYS 


READ(1,5,DATA=5001) 

1 {((ALITH  (ID,  IC,  IB)  ,ID*1,40)  ,IC*1,40)  ,IB*1,1) 

5 FORK  AT ( 1 0F8. 1 ) 

GO  TO  5006 

5001  WRITE (2 , 5002)  ID,  IC,  IB,  KEY 

5002  FORMAT  (1H0,  4110) 

STOP 

5006  KEY  =»  9 

BEAD  (1,7 ,DATA=50 01 ) 

1 ( ((ATHIC(ID,IC,IB) ,ID=1,40) ,IC  =1 , 40) , IB= 1 , 1 ) 

7 FORMAT  (10  (F3.1 ,1 X) ) 

IF(IMONY)  40  1,  402,  401 
402  IB=1 

DO  379  IC1=1, 21, 20 
WRITE  (2,358) 

358  FORMAT  (1H1, 3X, ’ARRAY  OF  SEDIMENTATION  BATES  BEFORE  INTERPOLATION*) 
WRITE  (2,351)  (IC,IC=IC1 ,IC1 *19) 

WRIT  E (2,  35  2)  (ID, ( ATHIC (ID, IC, IB) ,IC  = IC 1 ,IC 1 ♦ 19)  ,ID=1,40) 

379  CONTINUE 


FILL  IN  THE  BLANKS,  USING  LINEAR  INTERPOLATION 
INTERPOLATE  ATHIC  ARRAY  IN  IC  DIRECTION  WHERE  ZEROS  APPEAR. 
WILL  NOT  INTERPOLATE  FROM  ZERO  ON  LEFT  EDGE  OR  TO  ZERO  ON 
RIGHT  EDGE. 


401  IB*  1 

DO  7999  ID= 1 , 40 
DO  7999  IC=1  ,40 

IF  (ATHIC  (ID, IC, IB))  7899,7000,7899 
7000  IF  (IC-  1 ) 7899,7899,7050 
7050  IF  (ATHIC  (ID, IC~ 1 , IB) ) 7899,7899,7075 

7075  KOONT= 1 

7080  IF (ATHIC (I D, IC+ KCU NT, IB) ) 7200.7095,7200 

7095  IF  (IC+KOUNT-40)  7096,7200,7200 

7096  KOONT=KOUNT» 1 
GO  TO  7080 

7200  IF  (ATillC  (ID,  IC*KOUNT,  IB)  ) 7799,7799,7210 

7210  DD*  (ATHIC  (IT.IC  + KOUNT, IB)  - ATHIC  (ID , IC- 1 , IB) ) / (KOUNT* 1) 
DO  7299  KONT2=0, KOUNT-1 

7299  ATHIC(ID,IC+KONT2,IB)  * ATHIC  (ID, IC*KONT 2- 1 , IB)  ♦ DD 

7799  IC*IC+KOUNT 
7899  CONTINUE 
7999  CONTINUE 

IF  (IMONY)  403,404,  403 
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404  DO  377  IB*  1,  1 

DO  375  IC1  = 1 ,21,20 
HRITE(2,  353) 

353  FORMAT (1H1 , 3 X , • ARRAY  OF  CARBONATE  PERCENTAGES') 

WRITE  (2, 351)  (IC, IC*IC1 , IC1 ♦ 19) 

351  FORMAT  (3 X , ' IC * , 20  (3 X , 12 , 1 X)  , /•  ID') 

WRIT E (2,  352)  (ID,  (ALITH(ID,IC,IB) , IC=IC 1 , IC1 ♦ 1 9) ,ID*1,40| 

352  FORMAT(1X,I2,2X,20F6.  1) 

375  CONTINUE 

DO  376  IC1  = 1 ,21,20 
WRITE  (2,  354) 

354  F0RMAT(1H1,3X, 'ARRAY  OF  SEDIMENTATION  RATES') 

WRIT  E (2, 351)  (IC,IC=IC1 ,IC1 *19) 

WRITE  (2,  352)  (ID,  ( ATHIC (ID, IC, IB) ,IC=IC 1 ,IC 1 ♦ 1 9) ,ID=1,40) 

376  CONTINUE 

377  CONTINUE 


403  CONTINUE 

BBPP1  * (TEST3*R  ATE4  ♦ (TEST2-TEST3) *RATE3»  (TEST1-TEST2) * RATE2) * 

1 .0174533 

..ROTATE  EMPEROR  POLE  THE  DISTANCE  BBPP1  AROUND  HAWAIIAN  POLE. 

CALL  ROCCW  (BLAD1  , B LOD 1 , BL AD2 , BLOD2 , BBPP 1 , BA 2, BO 2) 

FILL  PERCENT  ARRAYS  WITH  LARGE  NUMBER  SO  THAT  ****  WILL  INDICATE 
NO  DATA  AS  WELL  AS  ESP  = 0 OR  PERCENT  VERY  LARGE 
DO  618  1=1,25 
DO  618  J=1,45 
PD  (I , J) =99999999. 

618  PC(I,J)=99999999. 

LOOP  THROUGH  ALL  HOLES 

DO  6000  L= 1, LL 

C. ..READ  LAT.,  LONG.,  AGE  OF  BASEMENT,  DEPTH  OF  WATER  OF  STARTING  POINT 
READ  ( 1 , 50)  CLAD1,CL0D1,AG£8  ,DEPH, NSITE 

50  FORMAT (3F10. 2,20 X,F1 0.2, IX, 15) 

WRIT E ( 2,  751) 

751  FORMAT  ( 1 H 1 ) 

WRIT  E (2,  608)  CLAD  1,  CLOD1,  AGEB 
WRITE  (2,  610)  DEPH, NSITE 

608  FORMAT  (3X, 'INITIAL  LAT  (CL A D 1 ) IS' , F6 . 2 , 3X , • I NITI A L LONGITUDE  (CLOD 
1)  IS* ,F8.2,3X, 'BASEMENT  AGE  (AGEB)  IS',F6.1) 

610  FORMAT (3X, 'DEPTH  OF  WATER  IS ' , F8 . 2, 70X , • SITE  NUMBER  IS  ',15) 

C., .DETERMINE  CRUSTAL  CORRECTION. 

CALL  CRUST  (AGEB,  PDEPH) 

DEPHC=  DEPH-PDEPH 
WRIT  E (2, 359) 

359  PORKAT  (/) 

WR ITE (2 , 708)  PDEPH,  DEPH,  DEPHC 

708  FORMAT (3X, 'PDEPH  = • , P8. 2, 1 0X , 'DEPH* • , F8. 2, 1 0X , • DEPHC  =',F8.2) 

WRITE  (2,51) 

51  FORMAT [///) 

WRITE  (2,  711) 

711  FORMAT  (29X , ' INT ' , 4X, 'DHCOMP' , 3X ,' COMPACTED ', 3X, 

1 'CUMUL' ,14X, 'SOUND') 

WRIT  E ( 2,  709) 

709  PORHAT (3 X , • L AT • , 4 X , • LONG • , 4X , • AGE' , 


IX, 'DEPT H', IX, 'CA CO 3' 


non 


1 2x, 'Thickness  thickness  thickness  id  ic  ib  veloc 

1 • ROPH  XV  DIFF  PCT  C*  XV  DIF  PCT',/) 

AGEC  * XGEB 

C...SET  VXRIOUS  XCCUHU LXTORS  TO  ZERO. 

TCXCO  =0.0 
TTHIC  = 0.0 
XGET  = 0.0 
COM  = 0. 

COHPXxO. 

C.,. DBTBRHINE  NUMBER  OF  ITERATIONS. 

XGEN  a TESTB 

IP  (TESTB  - XGEB)  712,712,714 
712  K1  * (TESTB/XINC  ♦ .0009) 

GO  TO  715 

714  K1  = (XGEB/  XINC  ♦ .0009) 

715  REXD  (1, 999)  (ROP HT (J) , J= 1. 45) 

999  FORMAT(20(F3.1,1X) } 

REXD  (1,  998)  (ROPHC(J)  ,0  = 1,45) 

998  FOBMXT  (26  (F2.0,1 X) ) 

LOOP  THROUGH  TIHE 

DO  2000  J=1,K1 
AGBT  * XGET  ♦ XINC 
BOPHX=0 
ROPHB=0 

C. . .DETERMINE  ANGULXR  DISTXNCE  OF  BOTXTION. 

IF  (AGET  - (TEST3  ♦ .0009))  200,  200,  210 
200  BBPP1  = AGET  * RATE4  * .0174533 
GO  TO  3005 

210  IF (XGET  - (TEST2  ♦ .0009) ) 220,220, 230 

220  BBPP1  * (TEST3  * RXTE4  ♦ (XGET  - TEST3)  * RATE3)  * .0174533 
GO  TO  3005 

230  IF  (XGET  - (TEST1  + .0009))  240,240,250 

240  BBPP1  = (TEST3*R XTE4  ♦ (TEST2  - TESTS) *RATE3*  (AGET-TEST2) *RATE2) * 

1 .0174533 
GO  TO  3005 

250  BBPPl  = (TEST3*H  ATE4  ♦ (TEST2-TEST3) *RATE 3* (TEST 1-TEST2)  *RATB2)  * 

1 .0174533 

CALL  ROCCH (BLXD1 ,BLOD1 ,CLXD1 ,CLOD1 , BBPPl, CLX D, CLOD) 

CXD  » CLAD 
COD  = CLOD 

BBPP2  = ((XGET  - TEST1)  *RATE  1)  *.0174533 
CALL  ROCCH  (B A2 , B02,  CXD,  COD,  BBPP2,CLAD,  CLOD) 

GO  TO  3006 

C.. .ROTATE  POINTS  YOUNGER  THAN  TEST1  YEARS  (HAWAIIAN  POLE  ONLY). 

3005  CALL  ROCCH (BL AD  1 , BLO D 1 ,CLAD1 .CLOD  1 , BBPPl, CLAD,  CLOD) 

3006  AGEC  = AGEC  - XINC 
CALL  CRUST (AGEC,  PDEPH) 

C.. .DETERMINE  DEPTH  OF  HATER. 

PDEPH  * PDEPH  ♦ DEPHC  ♦ COHPA 
D = PDEPH 
DO  1000  IT ER=  1 , 3 

C. . . DETERMINE  SUBSCRIPTS  FOR  ATHIC  AND  ALITH. 

CALL  TABLE  (AGEN,  PDEPH,  CLAD,  TESTA ,RFLAG , ID ,IC,IB) 

C. ..INTERPOLATE  TO  FIND  RATE  OF  SEDIMENTATION. 

CALL  XLOOK (ATHIC, ID, IC, IB, PDEPH, CLAD,  XTHIC) 

C.. .INTERPOLATE  TO  FIND  CARBONATE  PERCENT. 

CALL  XLOOK (ALITH, ID, IC, IB, PDEPH, CLAD, XLITH) 

C...END  LOOP  IF  HFLAG (SET  BY  TABLE)  IS  1,  BECAUSE  POINT  LIBS  OUTSIDE  GRI 
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C AND  ATRIC) . 

IF(MPLAG)  902,  902,  6000 

C... REDUCE  UNCOMPACTED  INTERVAL  THICKNESS  AND  PERCENT  CARBONATE  BT 
C... AMOUNTS  IN  ARRAYS  CARRY  AND  CCARY.RESP.  IF  ARRAYS  ARB  ZERO  THEN 
C... THERE  IS  NO  REDUCTION. 

902  IF  (AGET- 40.)  910,909,909 

909  XTHIC  = XTHIC  - XTHIC  * CARRY  (41) 

XLITH*XLITH-XLITH*CCA8 Y (41) 

GO  TO  601 

910  I = IFIX (AGET  ♦ .00001)  ♦ 1 
FAGET  = FLOAT  (I  - 1) 

XI  = AGET  - FAGET 

XII  = (CARRY  (I)  - CARRY  (I  ♦ 1))*  XI 

XIII  = CARRY  (I)  - XII 

XIIC= (CCAR Y (I) -CCARY  (I*1))*XI 
XIII C=CC ARY (I) -XIIC 
XLITH*-XLITH*XIIIC*XLITH 
XTHIC  = XTHIC  * (-XIII)  ♦ XTHIC 
801  FTHIC  = XTHIC 

INT  = (IFIX  ( FTHIC  + ,00001>)/25.  ♦ 1 
FINT  = FLOAT  (INT  - 1)  * 25. 

FI  = (FTHIC  - FINT)/  25. 

FII  = (CORTT  (I  NT ) - CORTT (INT  ♦ 1))  * FI 
Fill  = CORTT  (IHT)  - FII 
FTHIC  = FTHIC  * Fill 

1000  PDEPH  * D ♦ FTHIC  - . 24* 1 . 6* FTHIC/2 . 3 
917  TTHIC  * TTHIC  * XTHIC 
C. ..COMPACT  THICKNESS 

INT  = (IFI X ( TTHIC  ♦ .000001))/  25.  ♦ 1 
FINT  * FLOAT  (INT  - 1)  * 25. 

FI  = (TTHIC  - FINT) /25. 

FII  = (CORTT  (INT)  - CORTT  (INT  *1))  * FI 
Fill*  CORTT ( INT)  - FII 
COMPA  * TTHIC  * Fill 

C... DETERMINE  ADJUSTED  INTERVAL  THICKNESS. 

CINT  = COMPA  - COM 
COB  * COBPA 

C.. .DETERMINE  VELOCITY. 

INT  * IFIX( (COMPA  ♦ .00001) /50.)  ♦ 1 
FINT  = FLOAT  (INT  - 1)  * 50. 

FI  * (COBPA  - FINT) /50. 

FII  = (VELA  (INT)  - VELA (INT  *1))  * FI 

VEL  = VELA  (INT)  - FII 
IF  (J.EQ.  1 .OR.  J .GE.K1)  GO  TO  956 
ROPH A= (ROPHT  (J-1) *ROPHT  (J+1) -2*ROPHT (J) )/3 
ROPH  B~  (P.OPHC  ( J-1)  ♦ROPHC(J*1)  «2*ROPHC  (J)  )/3 
956  ROPHA=ROPHA*ROPHT (J) 

ROPHB  = ROPHB*  ROPilC  ( J) 

CDIPF=XLITH-ROPHB 
TDIFP=XTHIC-KOPH A 
IP  (XTHIC)  960,960,950 
950  PDIFF=TDIFF/XTHIC*100. 

GO  TO  970 

960  PDIFF=99999999. 

970  PD(L,J) *PDIPP 

NS  (L)*NSITE 
IF(XLITH)  961,962,961 

961  PCDIF=CDIFF/XLITH* 100 
GO  TO  963 

962  PCDIF=999999999. 
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963  PC  (L,  J)  *PCDIF 

IP  (CLOD)  90S,  905,  906 

905  ACLOD  * CLOD 
GO  TO  90  1 

906  IP  (CLOD  - 180.)  905,  905,  907 

907  ACLOD  * CLOD  - 360. 

901  IP  (IAGE5)  921,921,924 
924  IAGET=IFIX  (AGET) 

IF (HCD (I AGET, I AG E5) ) 2000,921,2000 
C.. .OUTPUT  PORWARD  TRACKS 

921  WRITE(2, 900)  CLAD,  ACLOD,AGZT,PDEPH,XLITH, XTHIC ,CI NT,  COHPA 
1 ,ID,IC,IB,VEL,BOPHT  (J) , ROPHA ,TDIFP, PDIFF, ROPHC(J) , ROPHB, CDIFF, 
1PCDIF 

900  FORMAT (1X,P6.2,1X,F7.2,1X,F5.1,1X,I4,2X,F4.1,3P10.1,  16, 

1 2 (IX, 12) ,2X,F5. 3.  IX, 2F5.1.F6. 1,14,213,214) 

2000  CONTINUE 
6000  CONTINUE 

END  OF  ALL  LOOPS 

...OUTPUT  PERCENT  ERROR  TABLES 
DO  930  J1=1 ,22,21 
WRITE(2,972)  (J,  J* J 1,  J 1*20) 

WRITE  (2,973)  (NS  (L) , (PD(L,J)  , J=J1, J1*20) ,L=1,LL) 

930  CONTINUB 

DO  933  J1=1 ,22,2  1 

WRIT E (2, 97 4)  (J, JsJ 1, J 1*20) 

WRITE  (2 , 97 3)  (NS  (L) , (PC(L,J)  , J 1 , J 1* 20) ,L*1,LL) 

933  CONTINUE 

972  FORMAT  ( 1 H 1 , • ERROR  IN  ONCOHPACTBD  INTERVAL  THICKNESS (PERCENT)  • , 
1 //, 7X, 2 115/) 

973  FORMAT  (IX, 13, 3X, 2115) 

974  FORMAT ( 1H 1, • EKBOB  IN  INTER VAL  CARBONATE  (PERCENT) ',/// 7X ,2 1 15/) 
CALL  EXIT 
END 
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APPENDIX  III 
Subroutines 

by 

J.  L.  Matthews 
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Subroutine  ALPHA  (V,  I,  J,  K,  L) 


Purpose 

Subroutine  ALPHA  is  used  by  ROPH  to  convert  variables  whose  value 
is  less  than  1000.0  from  floating  point  format,  rounded  to  the  nearest 
tenth,  to  alphameric  format. 

Variables 

In 

V = variable  in  floating  point  format  to  be  changed  to  alphameric  format 
Out 

I = subscript  for  the  unit's  division 
J = subscript  for  the  ten's  division 
K = subscript  for  the  hundred's  division 
L = subscript  for  the  tenth's  division 

Method 

Subroutine  ALPHA  rounds  the  incoming  variable  (floating  point)  to 
the  nearest  tenth  of  a unit,  then  by  a series  of  tests  determines  the 
hundred's  division,  the  ten's  division,  the  unit's  division,  and  the 
tenth's  division  of  the  variable  and  assigns  each  division  a variable 
subscript.  The  subscript  refers  to  a one-dimensional  array  of  integers 
and  characters  entered  with  an  alphameric  format  in  the  main  program. 
Subroutine  ALPHA,  by  determining  the  appropriate  subscript  for  each  divi- 
sion, allows  the  main  program  to  transform  the  number  to  alphameric  for- 
mat by  reference  to  the  alphameric  array.  As  an  example,  assume  that  95.2 
is  to  be  converted  from  floating  point  to  alphameric  format.  Array  IG 
(JI)  already  has  been  read  into  the  main  program  with  the  following  corres 
pondence  of  subscript  and  value: 
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Subscript  of  IG  1 2 3 4 5 6 7 8 9 10  11  12 

Alphameric  value  012345678  9 . * 


Subroutine  ALPHA  will  determine  that  the  hundred's  division  is  0 
and  assign  the  subscript  1 to  variable  K,  the  ten's  division  is  9 and 
assigns  the  subscript  10  to  variable  J,  the  unit's  division  is  5 and 
assigns  the  subscript  6 to  variable  I,  and  the  tenth's  division  2 and 
assigns  the  subscript  3 to  variable  L.  The  main  program  will  write  5 
variables  of  one  character  each  --  IG(J1),  IG(J2),  IG(J3),  IG(ll),  and 
IG(J4),  where  J1  corresponds  with  K,  J2  with  J,  J3  with  I,  and  J4  with 
L.  Note  that  the  fourth  variable  always  has  the  subscript  11,  and  causes 
a decimal  to  be  printed.  If  a value,  sent  to  subroutine  ALPHA,  is  equal 
or  greater  than  1000.0,  all  subscript  variables  of  subroutine  ALPHA  are 
assigned  the  value  12,  which  causes  asterisks  to  print  for  all  4 variables. 

I 

Subroutine  COUNT  (BL0DF,  CL0DF,  BBPPF,  BBF,  IFLAG,  BBPF,  CLADF,  BLADF) 


Purpose 

Subroutine  COUNT  is  called  by  subroutine  R0CCW  and  determines  the 
included  angle  ABC'  of  the  spherical  triangle  ABC'  (Figure  1.4)  prepara- 
tory to  solution  of  the  triangle  by  subroutine  CPRIM  and  sets  flags  used 
by  R0CCW  to  determine  the  new  latitude  and  longitude  of  the  rotated  point. 


Variables 


BL0DF  = longitude  of  pole 

CL0DF  = longitude  of  point  to  be  rotated 

BBPPF  = angular  distance  of  rotation  (CBC',  Figure  1.4) 

BBF  = included  angle  of  spherical  triangle  (ABC,  Figure  1.4) 


*L. 
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CLADF  = latitude  of  point  to  be  rotated 
BLADF  - latitude  of  pole 

Out 

BBPF  = angle  (ABC1,  Figure  1.4)  to  be  used  in  solution  second  spherical 
triangle  by  subroutine  CPRIM 

I FLAG  = value  is  set  at  -1,  1,  or  2;  used  by  subroutine  ROCCW  to  deter- 
mine the  latitude  and  longitude  of  rotated  point. 

Method  ^ 

If  the  position  of  the  point  C relative  to  the  pole  B is  as  illus- 
trated in  Figure  1.4,  determination  of  angle  ABC'  is  made  by  addition  of 
angles  ABC  and  CBC';  however,  should  the  relative  positions  of  the  point 
and  pole  be  interchanged,  subtraction  is  appropriate.  Additional 
complications  will  arise  if  the  rotated  point  crosses  the  line  of  longitude 
at  360°,  or  if  it  crosses  the  line  of  longitude  through  the  pole  of 
rotation. 

Figure  3.1  illustrates  the  algorithms  used  by  COUNT  to  calculate 
angle  ABC' . 

Subroutine  CPRIM  (AB,  BS.  CS,  BB,  AS,  CB) 

Purpose 

This  subroutine  solves  a spherical  triangle,  given  two  sides  and 
an  included  angle.  Used  in  programs  ROPH  and  ESP  for  rotation  of  points. 
Calling  subroutine  is  ROCCW. 

Variables 

In 


j 
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Figure  3.1.  Stereographic  projection  of  spin  axis  A,  pole  of  rota- 
tion B,  and  Point  C to  be  rotated  to  point  C',  showing  various  deter- 
minations of  angle  ABC',  given  angle  ACB  and  angle  CBC'.  The  solu- 
tion for  angle  ABC1  depends  on  the  longitude  of  points  C and  C' 
relative  to  longitude  of  pole  B.  The  latitude  of  the  points  and 
pole  must  be  considered  where  points  C and  C'  lie  on  opposite  sides 
of  the  line  of  longitude  passing  through  point  B and  where  the  longi- 
tude of  C and  B are  equal. 


iniT'~iaii«i  fr  ii'  ~'  * 
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AB  = included  angle  of  spherical  triangle  (radians) 

BS  = side  adjacent  (radians)  angle  AB 
CS  = side  adjacent  (radians)  angle  AB 

Out 

BB  = angle  (radians)  opposite  side  BS 
AS  = side  (radians)  opposite  angle  AB 
CB  = angle  (radians)  opposite  side  CS 

Subroutines 

Called  by  subroutine  ROCCW,  requires  functions  XTAN  and  FT AN. 

Method 

This  subroutine  solves  the  following  equations: 

tan  1/2  (C  + B)  = [cot  1/2  A cos  1/2  (c  - b)]/cos  1/2  (c  + b) 

tan  1/2  (C  - B)  = [cot  1/2  A sin  1/2  (c  - b)]/sin  1/2  (c  + b) 

tan  1/2  a = [tan  1/2  (c  - b)  sin  1/2  (C  + B)]/sin  1/2  (C  - B) 

where: 

A,  B,  and  C = the  three  angles  of  the  spherical  triangle 
a,  b,  and  c = the  sides  of  the  spherical  triangle,  where  side  a 
is  opposite  A,  b opposite  B,  and  c opposite  C. 
Subroutine  CRUST  (GAGE,  WATER) 

Purpose 

This  subroutine  determines  the  elevation  of  the  crust  from  the  sub- 
sidence curve,  given  the  age  of  the  crust.  Subroutine  CRUST  is  called  by 
programs  ROPH  and  ESP. 


Variables 


In 

GAGE  = age  (my)  of  crust 
Out 

WATER  = elevation  (m)  of  crust 
Method 

Subroutine  CRUST  solves  one  of  three  polynomial  equations  to  deter- 
mine tne  elevation  of  the  crust,  using  the  subsidence  curve  of  Sclater 
and  Detrick  (1973).  Dividing  the  subsidence  curve  in  three  parts  --  0 
to  25  mybp,  25  to  50  mybp,  and  older  than  50  mybp  --  gives  a better 
fitting  of  the  data  than  a fitting  by  a single  polynomial. 

The  equations  are: 

0-25  mybp 

Q = 2716.3934  + 143.5062Y  - 9.21059Y2  + 0.3957787Y3  - 0.0063777Y4 
25  - 50  mybp 

Q = 3381.7032  + 35.057943Y  - 0.0251748Y2 
> 50  mybp 

Q = 2577.2461  + 80.83089Y  - 0.7438308Y2  + 0.0031010Y3  - 0.0000049Y 

where: 

Q = elevation  of  crust 


Y = age  of  crust 
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Subroutine  ROCCW  (BLADX,  BLODX,  CLADX,  CLODX,  BBPPX,  XCLAD,  XCLOD) 

Purpose 

Given  the  latitude  and  longitude  of  a point  on  a sphere,  the  pole 
of  rotation,  and  the  angular  distance  of  rotation,  this  subroutine  rotates 
the  point  counterclockwise  and  determines  the  new  latitude  and  longitude 
of  the  point. 

Variables 

In 

BLADX  = latitude  (deg.)  of  pole  of  rotation 
BLODX  = longitude  (deg.)  pole  of  rotation 
CLADX  = latitude  (deg.)  point  to  be  rotated 
CLODX  = longitude  (deg.)  point  to  be  rotated 

BBPPX  = angular  distance  (deg.)  around  pole  points  are  to  be  rotated 
Out 

XCLAD  = latitude  (deg.)  of  point  after  rotation 
SCLOD  = longitude  (deg.)  of  point  after  rotation 

Subroutines 

Called  by  program  ROPH  and  ESP.  Calls  subroutines  COUNT  and  CPRIM, 
which  requires  functions  XTAN  and  FTAN. 

Method 

Subroutine  ROCCW  rotates  a point  counterclockwise  and  determines 
its  subsequent  latitude  and  longitude  by  solution  of  2 spherical  triangles, 
given  2 sides  and  the  included  angle  of  each  triangle.  In  Figure  1.4, 
the  included  angle  CAB  is  the  absolute  difference  of  the  longitude  of 
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point  C and  pole  B.  Side  b is  90°  minus  the  latitude  of  point  C,  and 
c is  90°  minus  the  latitude  of  pole  B.  Subroutine  CPRIM,  called  by  ROCCW, 
uses  the  included  angle  CAB  and  the  2 sides,  b and  c,  to  solve  for  side  a 
and  angle  ABC  of  triangle  ABC.  For  triangle  ABC1,  side  a has  the  same 
value  as  before.  Because  rotation  around  pole  B requires  that  point  C 
trace  a line  of  latitude  around  the  pole,  side  a and  a1  must  be  equal. 
Angle  ABC  is  known  from  the  solution  of  triangle  ABC.  Angle  CBC'  is 
the  angular  distance  of  rotation,  determined  by  the  main  program,  and  is 
equal  to  the  rate  of  rotation  (deg/my)  times  the  duration  (my)  of  the 
rotation;  thus,  the  included  angle  ABC'  equals  angle  ABC  plus  angle  CBC', 
the  addition  being  done  by  subroutine  COUNT  (for  some  rotations,  angle 
ABC'  is  not  a simple  addition  of  angles).  Subroutine  CPRIM  solves  this 
second  triangle,  ABC',  using  side  a',  side  c,  and  included  angle  ABC' 
to  give  side  b'  and  angle  C'AB.  The  latitude  of  C'  is  90°  minus  side  b ' ; 
the  longitude  of  C'  is  the  longitude  of  pole  B plus  or  minus  angle  C'AB. 

Subroutine  TABLE  (AGE,  DEPTH,  CLAD,  TESTB , MFLAG,  ID,  1C,  IB) 

Purpose 

Using  values  determined  by  program  ESP  for  depth  of  water,  latitude 
and  age  of  the  point,  this  subroutine  transforms  the  values  to  subscripts 
ID,  IC,  and  IB  for  use  by  arrays  ALITH  and  ATHIC. 

Variables 

In 

AGE  = age  (mybp)  of  the  point  for  the  current  iteration.  This  variable 
is  compared  with  variable  TESTB  to  set  the  value  of  variable  IB. 
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DEPTH  = depth  of  water  (m)  for  point  after  rotation. 

CLAD  = latitude  (deg.)  of  point  after  rotation. 

TESTB  = age  (mybp)  at  which  shift  is  made  from  first  set  of  grids  to 
second  set. 

Out 

MFLAG  = a flag  that  is  set  if  a rotated  point  plots  off  the  paper. 

Used  by  main  program  to  terminate  loop. 

ID  = water  depth,  coded  in  100  m intervals  from  2300  to  6300  m. 

Value  of  2300  m is  1,  6300  m is  40. 

IC  = latitude,  coded  in  1°  intervals  from  20°N  to  20°S.  A latitude 
of  20°N  is  1 , 20°S  is  40. 

IB  = grid  number,  one  set  of  2 grids  is  used,  but  by  changing  the 

read  statement,  2 sets  could  be  used.  Variable  TESTA  specifies 
which  set  of  grids  is  to  be  used. 

Features 

Prints  warnings  if  any  of  the  points  fall  outside  the  limits  of  the 
latitude-depth-of-water  grid.  Sets  flag  (MFLAG)  that  is  used  by  main  pro- 
gram to  terminate  loop  if  point  falls  outside  the  grid. 

Method 

Permissible  values  for  the  depth  of  water  range  from  2300  m to  6300  m. 
The  depth  range  2300  to  2400  m is  assigned  a subscript  ID  = 1.  Successive 
intervals  of  100  m receive  a successively  higher  index  number,  progressing 
to  subscript  of  40,  which  defines  the  range  6200  to  6300  m.  The  subscript 
is  determined  by  multiplying  the  depth  of  water  times  .01,  truncating  the 
decimal  fraction,  and  converting  the  result  to  an  integer,  which  sets  the 
subscript  ID. 


_ ».  Mt  .<4  j. 
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Subscript  IC  marks  the  latitude,  19°N  to  20°N,  receiving  the  sub- 
script 1,  19°  to  20°S  the  subscript  40.  To  accomplish  this  transforma- 
tion, the  latitude  is  stripped  of  its  decimal  fraction.  If  the  latitude 
is  southern,  21  is  added  to  the  absolute  value  of  the  stripped  transfor- 
mation; if  northern,  the  transformed  value  is  subtracted  from  20. 

Subscript  IB  denotes  the  set  of  grids  to  be  used.  The  value  of 
IB  will  be  1 unless  the  age  of  the  iteration  point  exceeds  the  test  value 
TESTA,  which  is  entered  in  the  program  ESP. 

Subroutine  XL00K  (VAL,  ID,  IC,  IB,  PHEPH,  CLAD,  VVAL) 

Purpose 

The  subroutine  is  used  by  ESP  to  interpolate  values  on  the  latitude- 
depth-of-water  grids  ALITH  and  ATHIC. 

Variables 

In 

VAL  = value  extracted  from  either  array  ALITH  or  ATHIC 

ID  = subscript  for  depth  of  water  used  by  VAL 

IC  = subscript  for  latitude  used  by  VAL 

IB  = subscript  for  grid  number,  which  depends  on  age  used  by  VAL 

PDEPH  = depth  of  water  of  point 

CLAD  = latitude  of  point 

Out 

VVAL  = interpolated  solution  from  either  ALITH  or  ATHIC 
Method 

The  arrays  ATHIC  and  ALITH  are  constructed  so  that  the  values  for 
each  set  of  coordinates  marks  the  center  of  a grid.  Figure  2.1.  Thus, 
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for  example,  in  Figure  2.1,  12.2  (ID  = 19  and  IC  = 23)  indicates  the 
rate  of  deposition  at  4150.0  m depth  of  water  and  =2.5°  latitude.  Any 
point  not  at  the  center  of  this  or  any  other  grid  requires  interpola- 
tion. Figure  2.2,  an  enlarged  portion  of  Figure  2.1,  shows  the  direc- 
tion in  which  interpolation  must  move,  if  the  point  c falls  in  the  first 
quadrant  of  the  grid  ID  = 19,  IC  = 23.  As  an  example,  consider  the  depth 
coordinate  of  point  c to  be  4150.0  m and  the  latitude  to  be  -2.75°.  Inter- 
polation will  occur  between  ID  = 19,  IC  = 23,  and  ID  = 19,  IC  = 23  v/hich 
establishes  the  value  of  point  a.  The  value  of  point  b is  interpolated 
between  ID  = 18,  IC  = 23  and  ID  = 18,  IC  = 24.  Now  interpolation  between 
a and  b gives  the  answer,  the  value  at  point  c.  Points  that  plot  in  the 
other  quadrants  follow  the  same  scheme,  although  a different  set  of  4 co- 
ordinates must  be  used  for  each  quadrant: 


Quadrant 

1 

2 

3 

4 


Point 

Grids 

a 

ID 

and 

IC,  ID 

and 

IC+1 

b 

ID- 

-1  and  IC, 

ID-1 

and 

IC+1 

a 

ID 

and 

IC,  ID 

and 

IC-1 

b 

ID- 

-1  and  IC, 

ID-1 

and 

IC-1 

a 

ID 

and 

IC,  ID 

and 

IC-1 

b 

ID+1  and  IC, 

ID+1 

and 

IC-1 

a 

ID 

and 

IC,  ID 

and 

IC+1 

b 

ID+1  and  IC, 

ID+1 

and 

IC+1 

Function  FTAN(D) 

Purpose 

To  compute  tangent  of  an  angle  from  the  sin  and  cos  of  the  angle; 
used  by  machines  with  no  tangent  function. 
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In 

D = angle  in  radian 
Out 

FTAN  = tangent  of  angle  in  radians 
Function  XTAN(XT) 

Purpose 

To  convert  arc  tangent  with  negative  value  to  positive  angles  less 
than  180° 

Variables 

la 

XT  = angle  in  radians 
Out 

XTAN  = angle  ( radians)  that  is  positive  and  equal  to  or  less  than  90° 
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SUBROUTINE  ROCCW  (BLADX , BLODX  .CLADX,  CLGDI , BBPPX,  (CUD,  XCLOD) 
IF(CLODX)  400.400,410 
400  CLODX  * CLOD X ♦ 360. 

410  IP  (BLODX)  420,  420,  430 

420  BLODX  = BLODX  ♦ 360. 

430  AB  = ABS (BLODX  - CLODX) 

IF(AB)  422,421,422 

421  AB  = .00001 
GO  TO  490 

422  IF (AB  - ISO.)  460,470,480 
460  AB  = AB  * .0174533 

GO  TO  490 

470  AB  = 179.99999  * .0174533 
GO  TO  490 

480  AB  * (360.  - AB)  * .0174533 
490  BS  = (90.  - CLAD  X)  * .0174533 
CS  = (90.  - BLADX)  * .0174533 
CALL  CPRIB  (AB,BS,CS,BB, AS,CB) 

BBX  = BB 

CALL  COUNT (B LODX .CLODX , BBPPX , BBX, IFLAG, BBPX, CL ADX, BLADX) 

CSP  = CS 
ASP  * AS 
BBP  = BBPX 

CALL  CPRIH  (BBP, CSP, ASP.CBP.BSP, ABP) 

IF (IFLAG)  1400,1500,1600 
1400  XCLOD  = BLODX  - (ABP/  .0174533) 

GO  TO  1616 

1500  XCLOD  = BLODX  - 180. 

GO  TO  1616 

1600  IP (IFLAG  - 1)  1610,1610,1620 
1610  XCLOD  * BLODX  ♦ (ABP/  .0174533) 

IF (360.  - XCLOD)  1615,  1615,1616 

1615  XCLOD  * XCLOD  - 360. 

GO  TO  1616 

1620  XCLOD  * BLODX 

1616  XCLAD  * 90.  - (BSP  / .0174533) 

RETURN 

END 
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SUBROUTINE  COUNT (BLODF.CLODF, BBPPF, BBP, IFLAG , BBPF,CLADF,BLADF) 
503  DIFF  a CLODP  - BLODF 
ADIFF  * ABS  (DIPF) 

IF(DIPF)  500,  510,  520 
500  IF  (ADIFF  - 180.)  530,  560,  550 
530  BBPP  * BBK  ♦ BBPPP 

IF  (BBPF  - (180.  * .0174533))  601,605,610 
601  IFLAG  = -1 
GO  TO  600 

605  BBPF  = 179.999  * .0174533 
IFLAG  = 2 
GO  TO  600 

610  BBPF  = (360.  *.0174533)  -BBPF 
IFLAG  » 1 
GO  TO  600 

520  IF  (ADIFF  - 180.)  550,  560,  530 

550  IF  (BBPPF  - BBF)  551,  552,  553 

551  BBPF  * BBF  - BBPPF 
IFLAG  « 1 

GO  TO  600 

552  BBPF  * .00001 
IFLAG  a 0 

GO  TO  600 

553  BBPF  » BBPPF  - BBPF 
IFLAG  = -1 

GO  TO  600 

510  IF  (B LADF  - CLADF)  560,  570,  540 
540  BBPF  = (180.  * .0174533)  - BBPPF 
IFLAG  * 1 
GO  TO  600 
560  BBPF  = BBPPF 
IFLAG  = -1 
GO  TO  600 
570  BBPF  a o.O 
600  RETUBN 
END 


SUBROUTINE  CRUST  (GAGE,  WATER) 

IP  (GAGB  - 25.)  100,  100,  200 

100  HATER  « 2716.3934  ♦ 1 4 3 . 5062 *G AG E - 9. 21059*GAGE**2.  ♦ .3957787 
1*G*GE**3.  - . 0063777*GAGB**4, 

RETURN 

200  IF  (GAGE  - 50.)  300,300,400 

300  HATER  a 3381.7032  ♦ 35. 05794 3*G AGS  - . 0251748*GAGE**2. 

RETURN 

400  HATER  a 2577.2461  ♦ 80.830B9*GAGE  - . 743830B*GAGE**2.  ♦ .0031010* 
1 GAGE**3 . - .000004931 1*GAGE**4. 

RETURN 

END 


SUBROUTINE  TABLE  (AGE, DEPTH, CLAD, TBSTA , HFLAG , ID , IC,IB) 

HFLAG  = 0 

IF  (AGE  - (TBSTA  ♦ .00001))  20,20,10 
10  IB  *=  2 
GO  TO  30 
20  IB  * 1 

30  ID  * DEPTH  *.01 
ID  * ID  ~ 22 
IF  (ID)  820,  820,821 

820  WRITE(2,819) 

819  FORM  AT  (2  X , 'ID  IS  LESS  THAN  1,  DEPTH  IS  LESS  THAN  2300  HETERS*) 
HFLAG  * 1 
GO  TO  500 

821  IF (ID  - 40)  826,  826,  824 

824  WHITE (2,  825) 

825  FORMAT (2X, 'ID  IS  GREATER  THAN  40,  DEPTH  IS  GREATER  THAN  6300H.  •) 
HFLAG  = 1 

GO  TO  500 

826  IC  = CLAD 

IF  (CLAD)  100,101,101 

100  IC  = I ABS  (IC)  ♦ 21 
GO  TO  103 

101  IC  = 20  - IC 

103  IF  (IC)  810,810,812 
810  WRITE  (2,  813) 

,813  FORMAT (2X, 'IC  IS  NEGATIVE,  LAT.  OF  POINT  IS  GREATER  THAN  20  DEG*) 
HFLAG  = 1 
GO  TO  500 

812  IF  (IC  - 40)  500,500,815 
815  WRIT E (2,  814) 

814  FORMAT  (2X,*IC  IS  GREATER  THAN  40,  LAT.  OF  POINT  IS  LESS  THAN  -20  D 
1EG«) 

HFLAG  - 1 
500  RETURN 
END 
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SUBROUTINE  CPRIN  ( AB, BS, CS, BB, AS.CB) 
BPC  = (BS  ♦ CS)  / 2. 

BAD90  » 90.  * .0174533 


RAD  1 8 * 180.  * .0174533 
If  (BPC  - RAD90)  300,300,310 
300  SBPC  =>  SIN  (BPC) 

CBPC  = COS  (BPC) 

GO  TO  320 

310  SBPC  * SIN  (SADI  8 - BPC) 

CBPC  = -COS (RADI  8 - BPC) 

320  BBC  = ABS  (BS  - CS)  / 2. 

IF  (BBC  - RAD90)  330,330,340 
330  SBBC  = SIN  (BBC) 

CBBC  = COS (BBC) 

GO  TO  350 

340  SBI'C  = SIN  (RAD  18  - BBC) 

CBBC  = - COS  (RAD  1 8 - BBC) 

350  HAB  = AB  / 2. 

COTAB  = 1./FTAH  (HAB) 

XTAN1  =*  ATAN  ( (COTAB*SBHC)  /SBPC) 
XTAN2  = ATAN  ((COTAB  * CBBC)  / CBPC) 
IF  (CS  -BS)  355,355,356 

355  BB  *=  XTA N ( XTAN2)  + XTAN(XTANI) 

CB  = ABS (XTAN (XTAN2)  -XTAN ( XTAN 1 ) ) 

GO  TO  357 

356  CB  = XTAN  (XTAN 2)  ♦ XTAN (XTAN 1) 

BB  = ABS  (XTAN  (XTAN2)  - XTAN  (XTAN1) ) 

357  BBBC  = ABS  (BB  - CB)  / 2. 

BBPC  * (BB  ♦ CB)  / 2. 

IF  (RAD90  - BBKC)  360,  370,370 


COTAB 
XTAN  1 


360  SBBBC  = SIN  (RADI  8 - BBBC) 

GO  TO  380 

370  SBBHC  = SIN(BBtiC) 

380  IF  ( RAD90  - BBPC)  390,  400,400 
390  SBBPC  = SIN  (RADI  8 - BBPC) 

GO  TO  410 

400  SBBPC  = SI N (BBPC ) 

410  TBHC  = FT AN ( BBC) 

440  TAB  - (SBBPC  * TBHC)  / SBBHC 
IF  (TAB)  450,  460,  460 
450  AS  = 2.  * ATAN  (RADIS  ♦ TAB) 


RETURN 

460  AS  - 2.  * ATAN  (TAB) 


RETURN 

END 


O o o 
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FUNCTION  XTAN(XT) 

IF  (XT)  900,910,910 
900  XTAN  = (180.  * .0174533)  ♦ XT 
RETURN 

910  XTAN  = XT 
RETURN 
END 


FUNCTION  FTAN(D) 

IF  (D  - (90.  * .0174533))  600,  600,  610 
600  FTAN  = SIN (0)  / COS(D) 

RETUR  N 

610  FD  = (180.  * .0174533)  - D 
FTAN  * - (SIN  (FD ) / COS  (FD) ) 

RETURN 

END 


SUBROUTINE  A LPH A ( V, I , J , K , L ) 


| 


ROUND  OFF  TO  NEAREST  TENTH 
V = V + .05555 

C IF  V IS  GREATER  THAN  UR  EQUAL  TO  100  , **.*WlLl  PRINT 
IF (V  - 1000.  ) 40,55,55 
55  I = 12 
J = 12 
K = 12 
L = 12 


RETURN 

40 

I = 1 

45 

< 

II 

< 

- 100.00001 

IF  ( V ) 

50,  46 , 46 

46 

I = I 

+ 1 

GO  TO 

45 

50 

J = 1 

V = V 

+ 100.00001 

100 

V = V- 

- 10.00001 

1F  (V  ) 

300 , 200 , 200 

200 

J = J 

+ 1 

go  ru 

100 

300 

K = 1 
V = V 

+ 10.0000  1 

350 

V = V 

- 1.00001 

IF  (V) 

500,  400,  400 

400 

K = K 

+ 1 

GO  TO 

350 

500 

1.  = 1 

V = V 

+ 1.00001 

450 

V = V 

- .10001 

IF  (V) 

700,  600,  600 

600 

L = L 

+ 1 

GO  TO 

450 

700 

RETURN 

END 

\ 


' 
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SUBROUTINE  XLOOK  (VAL,  ID, IC, IB, PDEPH , CLAD, VVAL) 
DIMENSION  VAL(40,40,  1) 

IF  (ID  - 1)  5,5,10 
5 V V AL  = VAL  (ID, IC, IB) 

RBTURN 

10  IF  (ID  -40)  15,  5,  5 
15  IF  (IC  -1)  5,5,20 
20  IF (IC  -40)  30, 5, 5 
30  DEPH  = PDEPH  * .01 
IDEP  = IFIX(DEPH) 

FDEP  = FLOAT (IDEP) 

XFDEP  = Di’PH  - FDEP 
ICLAD  = IF1X  (CLAD) 

FCLAD  * FLOAT (ICLAD) 

XCDIF  = ABS  (CLAD  - FCLAD) 

IF  (XFDEP  - . 5)  40,  40,  140 
40  IF  (CLAD)  50,  60,  60 
50  IF  (XCDIF  - . 5)  70,  70,  110 
60  IF  (XCDIF  - .5)  110,  110,  70 

70  V A 1 = V AL ( ID,  IC,  IB)  - VAL(ID,  IC  - 1,  IB) 

VA 2 = VAL  (ID  - 1 , IC,  IB)  - V AL  (ID- 1,  IC-1,  IB) 
80  IF  (XCDIF  - . 5)  90,  90,  100 
90  FLAG  = 1.0 
FFLAG  = 1.0 
GO  TO  130 
100  FLAG  = -1.0 
FFLAG  = 1.0 
GO  TO  130 

110  VA  1 = VAL  (I D,  IC,  IB)  - VAL(ID,  IC»1,  IB) 

VA2  = VAL( ID-1, IC, IB)  - VAL (ID- 1 , ICO,  IB) 

120  IF  (XCDIF  - .5)  90,  90,  100 

130  VI  * VAL (ID,  IC,  IB)  ♦ (XCDIF  - .5)*VA1*FLAG 
V2  * V AL  (I  D- 1 , IC,  IB)  ♦ (XCDIF  - .5)*VA2*FLAG 
V V AL  = VI  ♦ (VI  - V2)  * (X  FDEP  - . 5)  *FFLAG 
RETURN 


140 

IF 

(CLAD)  150 

9 

150,  160 

150 

IF 

(XCDIF  - . 

5) 

170, 

170,  210 

160 

IF 

(XCDIF  - . 

5) 

210, 

210,  170 

170 

VA 

1 = VA  L (I D 

,IC,  IB) 

- V AL ( I D,  IC-1, 

, IB) 

VA2  =*  VAL  (ID 

♦ 1 

9 IC  0 

IB)  - V AL  (ID*  1, 

IC-1, 

180 

IF 

(XCDIF  - . 

5) 

190, 

190,  200 

190  FLAG  = 1. 

FFLAG  = -1.0 
GO  TO  230 
200  FLAG  = -1.0 
FFLAG  * - 1.0 
GO  TO  230 

210  VA  1 = VA  L (ID , IC,  IB)  - VAL(ID,  ICO,  IB) 

VA2  = VAL  (IDO,  IC.  IB)  - VAL(IDO,  IC*  1 , IB) 

220  IF  (XCDIF  - .5)  190,  190,  200 

230  VI  = V AL  (ID, IC,  IB)  ♦ (XCDIF  - .5)*VA1*FLAG 
V 2 = VA  L (I  DO , IC,  IB)  ♦ (XCDIF  - .5)*VA2*FLAG 
V V AL  = VI  ♦ (VI  - V2) * (XFDEP  - . 5) *FFLAG 
RETURN 
END 


fc..  . 
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APPENDIX  IV 

The  effect  of  error  in  basement  ages  and  pole  of  rotation 
on  back-tracking  and  forward-tracking  models 
in  the  central  equatorial  Pacific  (ROPH,  ESP) 

Richard  F.  Johnson 
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A 


C constant  aeptn  correction  in  subsidence  equation 

E elevation  above  basin  depth  (z  ) of  ridge  crest  when  C = 0 

h elevation  above  basin  depth  of  a point  on  a track 

hQ  initial  elevation  above  basin  depths  of  a point  on  a track 

hC(,D  elevation  above  basin  depth  of  carbonate  compensation  depth 
k time  needed  after  creation  for  oceanic  crust  to  subside  to 

an  elevation  above  basin  depth  of  e ^ ■ E 
P path  of  integration  of  ESP 

S(^,z)  scalar  field  of  sedimentation  rates,  such  as  Fig.  8 

Sq  average  sedimentation  rate  above  carbonate  compensation  depth 

S sedimentation  rate  when  tracks  cross 

x 

t time  before  present 

t,  basement  age 

D 

t time  before  present  at  which  a track  crosses  carbonate  compensation 


depth 

time  before  present  when  tracks  cross 
thickness  of  sediment  predicted  by  ESP 
depth  of  water  to  the  seafloor 

depth  at  which  rate  of  supply  of  calcium  carbonate  equals  rate 
of  dissolution 


basin  depth  (when  C = 0) ; i.e.,  depth  of  very  old  oceanic  crust 
error  in  basement  age 


error  in  predicted  thickness 

error  in  depth  of  a point  on  a track  at  a given  age 
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X* 

V 

X 

o 

X* 


4> 


A <r 


u* 


error  in  longitude  of  pole  of  rotation 

error  in  latitude  of  a point  on  a track  at  a given  depth 

error  in  latitude  of  pole  of  rotation 

error  in  rate  of  rotation  of  pole 

initial  longitude  of  a point  on  a track 

longitude  of  pole  of  rotation 

initial  longitude  of  a point  on  a track  in  coordinates  of  pole 
of  rotation 

latitude  of  a point  on  a track 
initial  latitude  of  a point  on  a track 
latitude  of  pole  of  rotation 

colatitude  of  a point  on  a track  in  coordinates  of  pole  of 
rotation 

relative  error  in  latitude  of  two  points  on  two  different 
tracks  due  to  error  in  pole  of  rotation;  gauged  along  a 
given  contour  of  sedimentation  rate 
rate  of  rotation  of  pole 
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Introduction 

There  are  four  main  sources  of  error  in  models  ROPH  and  ESP: 

1)  error  in  basement  ages;  2)  error  in  selection  of  a pole  of  rotation* ; 

3)  error  in  the  equation  we  use  to  describe  the  subsidence  of  the  seafloor 
as  it  moves  outward  from  the  spreading  center;  and  4)  error  in  sediment  data 
(viz.,  sedimentation  rates,  carbonate  percentages,  and  sediment  densities). 
In  addition,  ESP  is  also  subject  to  a contouring  error. 

Error  in  our  equation  of  subsidence,  error  in  sediment  data,  and  error 
in  contouring  are  important  but  difficult  to  quantify.  In  sediment  data, 
errors  in  assigning  absolute  ages  to  biostratigraphic  boundaries,  and  errors 
in  defining  these  boundaries  in  the  sediment  can  have  drastic  effects,  espe- 
cially on  sedimentation  rates  (see  van  Andel  and  others,  1975). 

Because  of  the  limited  information  available  on  the  other  sources  of 
error,  we  will  discuss  in  this  report  only  the  effects  of  errors  in  basement 
age  and  pole  selection. 

Model  ROPH 

Effect  of  error  in  basement  ages 

The  effect  on  ROPH  of  error  in  basement  ages  is  relatively  easy  to 

estimate.  Let  £ be  the  error  in  a given  basement  age.  This  error  will 
b 

cause  an  error  in  depth,  £ ^ = £z(t),  at  each  time  step,  t,  along  the  tracks 
in  ROPH.  We  now  approximate  the  equation  of  subsidence  of  the  seafloor  by  a 
simple  exponential  curve: 

z = z(t)  = z - E’exp(-(t  -t)/k)  + C,  (]) 

00  D 

*Here  and  hereafter,  "pole"  means  a vector  pole  of  rotation, 
characterized  by  a latitude,^*,  and  a longitude,  X*,  and  a rate,U>*. 
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where 

z = depth  of  the  seafloor, 
t = time  before  present, 
z^3  basin  depth,  a constant, 

E = ridge  elevation,  a constant, 
t^  = basement  age  at  t = 0, 
k = time  constant,  and 

C = zQ  - Zqq  + E-expf-t^/k),  a constant  depth  correction 
depending  on  the  local  depth,  zQ  = z(0). 

By  simply  guessing  a fit  to  Fig.  2a  of  Sclater  and  others  (1971,  p.  7894), 
we  find  that  approximate  values  for  the  constants  of  eq  (1)  are 
z00=  6000  m, 

E = 3300  m,  and 
k * 40  m.y. 


J 


From  eq  (1),  the  error  in  depth,  £^,  due  to  an  error  in  basement  age. 


E , must  be 
b 


VS. 


= ^E-exp(-(tb-t)/k)  - E-exp(-tb/k)J  • £l  - exp(- Et  /k)J 


= (z  - z)-  £ /k,  for  E <<  k. 

Tb  Zh 

The  fractional  elevation  error  has  the  particularly  simple  expression. 


(2) 


e„/(z  - z>  = e /k. 

z o t. 

b 


(3) 


*1 


For  example,  at  a depth  1000m  below  the  ridge  crest,  a 4 m.y.  error  in 
basement  age  will  result  in  an  error  in  backtracked  depth  of  100m. 


...<*•  _ 


Van  Andel  and  Bukry  (1973)  give  estimates  of  the  confidence  limits 
for  the  basement  ages  of  12  of  the  DSDP  sites  used  in  this  study  (74,  75, 
77,  78,  79,  80,  82,  159,  160,  161,  162,  and  163).  These  confidence  limits 
range  from  about  il  m.y.  for  hole  82  to  17  m. y.  for  holes  74  and  163. 

From  eq  (3)  then,  the  elevation  error  theoretically  ranges  from  2% 
for  hole  82  to  nearly  20%  for  holes  74  and  163.  For  hole  163,  the 
error  could  be  as  high  as  ±400  m at  the  ridge  crest;  this  probably 
represents  a maximum  error  among  the  23  DSDP  sites  of  this  study. 

Berger  (]973)  obtained  a similar  value  graphically. 

Effect  of  an  error  in  pole  of  rotation 

The  effect  of  an  error  in  choice  of  the  pole  of  rotation  is 
more  complicated  than  the  effect  of  an  error  in  basement  age.  An  error 
in  the  pole  displaces  tracks  on  the  latitude-depth-of-vater  grid  (Fig.  1) 
right  or  left.  The  amount  of  the  displacement  varies  along  a track's 
length  and  among  tracks.  For  a given  track,  the  displacement  is  zero  at 
the  initial  latitude  and  depth  and  a maximum  at  the  final  latitude  and 
depth. 

If  we  use  eo  (1)  to  relate  time  to  depth,  we  can  describe  a track 
by  its  latitude,  4 = 4(z;4*,  X*,  u>*,40,  X0,t(j),  as  a function  of  depth 
(the  independent  variable)  and  the  parameters,  4*,  X*,  ui*,  4a,  \0,  and 
t^,  where,  in  addition  to  the  symbols  already  defined,  4q  and  XQ  are 
the  latitude  and  longitude  at  t = 0.  From  spherical  trigonometry. 


we  can  deduce  this  function  to  be 
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4>  = arcsin 


A A 

sin  ^*>cos  i>  + cos  ^*.sin  i>  *cos(\  + OJ**t) 
P P P0 


where,  in  the  simplest  case, 
A 


A I 

= arccos  sin  ^*>sin  i>  + cos  ^*-cos  4>  • cos(A * - X ) 
Pi  o o a J 


X = arcsin  cos  <b  *sin(X*  - X )/sin  <t 


PJ 


t = k-log^Cz^  + C - z)/E 


+ V 


(4) 


Here  <t>  and  X are  the  colatitude  and  initial  longitude  of  a given 

P po 

track  in  the  coordinates  of  the  pole  of  rotation;  this  longitude 
is  measured  from  a prime  meridian  passing  through  the  positions  of  the 
poles  of  rotation  and  daily  spin. 

An  error  in  the  pole  of  rotation,  E^  or  E^,  causes  an  error 
in  track  latitude. 


(z ; , x*  ,uax,  «SQ , Xj . tb , ,e^A) 

i (z ; «**+£j* , X*+£x*  >w*+Eayc  ’ ’ Xo  ’ tb) 


- $(z  X*,oj*,<S  > X„>  t.)  • 

o o b 


(5) 


j 


► 

i t 


When  the  error  in  the  pole  is  small,  we  can  approximate  E^  by  the 
differential. 


A* = §f*d** + ff*dX* + t>dw*- 


The  differentials,  d^*,  dX*,  and  do )*,  are  now  the  errors  in  the  pole. 

Table  4.1  gives  some  typical  values  for  the  partial  derivatives 
in  eq  (6)  for  the  DSDP  sites  used  in  this  report.  The  size  of  these 
derivatives  indicates  the  sensitivity  of  the  tracks  to  small  changes 
in  the  choice  of  the  pole  of  rotation. 


(6) 


, * 
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For  tracks  no  longer  than  40  m.y.,  a shift  in  the  position  of 
the  pole  by  a degree  of  latitude  in  any  direction  will  displace  the 
latitude  of  the  tracks  by  less  than  a degree.  Thus,  the  tracks 
included  in  this  study  are  stable;  nothing  violent  occurs  by  shifting 
the  pole  slightly. 

Also,  note  from  Table  4.1  that  the  more  eastern  holes  (see 
Table  I)  are  insensitive  to  errors  in  pole  latitude;  that  increasing 
pole  longitude  decreases  track  latitude;  that  increasing  pole  latitude 
tends  to  increase  track  latitude;  and  that  errors  in  rotation  rate 
of  .l°/m.y.  have  the  same  magnitude  of  effect  as  shifting  the  pole 
by  a distance  of  1°  of  latitude. 

In  lieu  of  confidence  limits  for  our  tracks,  we  offer  Table  4.2. 
This  table  shows  the  range  of  track  latitudes  produced  from  five 
published  poles  of  rotation: 


Morgan  (1972) 

67°N, 

73°W, 

.85°/m.y. 

Clague  and  Jarrard 

(1973) 

72°N, 

83°W, 

•9°/m.y. 

Clague  and  Jarrard 

(1973)* 

69°N, 

68°W, 

•9°/m.y. 

Winterer  (1973)** 

67°N, 

45°W, 

l.l°/m.y. 

Minster  and  others 

(1974)*** 

67°N, 

59°W, 

•83°/m.y . 

(The  pole  used  for  ESP  and  ROPH  is  72°N,  83°W,  .81°/m.y.)  Unfortunately, 
only  Clague  and  Jarrard  (1973)  have  published  confidence  limits  for  their 
pole,  and  their  limits,  according  to  Jarrard  (personal  communication),  are 
not  calculated  from  reliable  statistics. 

*This  determination  omitted  the  "Line  Cross"  trend.  Also,Qthe 
rate  we  show  here  is  a time  average  of  what  they  published,  1.3  /m.y. 
for  0-20  m.y.  BP  and  .5°/m.y  for  20-40  m.y.  BP. 

**Winterer's  rotation  is  relative  to  the  athenosphere . 

***We  assume  Minster's  instantaneous  pole  holds  for  the  past  40  m.y. 


. <v, 
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For  the  holes  in  Table  1,  Winterer's  pole  consistently  yielded  the 
minimum  track  latitudes  of  Table  4.2;  Clague  and  Jarrard's  western 
pole  consistently  yielded  the  maximum  track  latitudes.  The  pole  we 
use  in  this  study  produces  more  northerly  tracks  than  the  majority 
of  the  published  poles. 

Note  also  from  Table  4.2  that  the  range  of  track  latitudes  varies 
by  a factor  of  2 or  3 from  hole  to  hole;  that  the  range  increases  more 
or  less  linearly  with  age;  and  that  the  range  is  as  much  as  10°  at 
40  m.y.  It  is  clear  from  this  last  fact  that  today  any  estimate  of  paleo- 
latitude  for  older  sediments  in  the  central  equatorial  Pacific  is  unreliable. 

Model  ESP 

To  evaluate  the  effect  of  error  in  basement  ages  and  pole  of  rotation 
for  model  ESP,  we  recall  first  that  ESP’s  main  function  is  to  calculate 
thickness,  T,  from  the  integral. 


= j S (^ , z)  dt 
*'P 


S(j5,z)  is  the  scalar  field  of  sedimentation  rates  found  by  smoothing 
data  along  tracks  calculated  by  ROPH;  P is  the  path  of  integration, 
which  in  fact  is  a track  given  by  eq  (4) . 

Errors  in  pole  selection  and  basement  age  affect  ESP's  calculation 
of  thickness  in  two  ways:  1)  Erroneous  tracks  calculated  by  ROPH  will 

give  erroneous  contours  of  S(<£,  z)  (contours  such  as  those  of  Fig.  8);  2)  erro- 


neous 
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However,  in  those  applications  where  ESP's  tracks  are  identical  to 
ROPH's,  such  as  in  the  present  study,  these  two  kinds  of  errors  tend 
to  cancel  each  other.  This  is  especially  true  in  the  case  of  an  error 
in  selection  of  a pole  of  rotation,  as  discussed  below. 


Effect  of  error  in  basement  ages 

Errors  in  basement  ages  affect  ESP  primarily  as  errors  in  the 

paths  of  integration.  The  effect  of  error  in  basement  ages  on  ROPH 

contours  is  presumably  lessened  by  smoothing. 

Consequently,  thickness  errors,  £ , caused  by  errors  in  basement 

age,  P , can  be  estimated  easily.  To  good  precision,  these  thickness 
b 

errors  are  given  by 


c = 91  a.  . 9_ 

T " 3tbdtb  ' atb 


U‘“- 


z)  dt 


dt. 


(8) 


Because  of  the  abruptness  of  the  carbonate  compensation  depth  (see  Fig.  8), 

we  can  approximate  S(^,z)  as  a step-function  in  depth:  Let  S(^,z)  be  a 

constant,  S for  all  depths  above  the  carbonate  compensation  depth,  z . 

° CCD 

Let  S(«S,z)  be  zero  below  the  carbonate  compensation  depth. 


With  this  approximation,  eq  (7)  becomes 


where,  from  eq  (1), 

CCCD  = Cb  + k*log((Zo  " ZCCD  + E-exp(-tb/k))/E), 
tCCD  being  the  time  at  which  a track  crosses  the  carbonate  compensation 


depth. 
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We  can  now  perform  the  differentiation  in  eq  (8)  : 


£T  = So-E-exp(-tb/k)-dtb/(Zoo  + C - ZcCD).  (9) 

Also,  the  fractional  error  in  thickness  must  be 

£X/T  = -E-exp(-tb/k)/((Zoo  + C - ZCCD) •l°g((zao  + C - Zca))/E)).dtb 

(10) 


If  we  let  h = h(t)  be  the  height  of  a point  above  the  basin  depth, 
Zqo  , these  equations  have  a simpler  appearance: 


£I  ■ 


and 


ex/T  = hQ/Ch 

= (ho/h 


CCD 

CCD 


•lo8(hCCD/E))-dtb 

) • (dtb/(tb-tCCD)) , 


where  hQ  = h(0)  and  h^  = Mt^)  . 

Table  4.3  shows  £ and  £ /T  for  various  values  of  z and  basement 

1 1 L>VjU 

age.  The  values  of  thickness  error  and  percent  thickness  error  are 
per  million  year  error  in  basement  age.  Note  that  the  percent  error 
has  a maximum  when  the  compensation  depth  is  both  shallow  and  deep, 
that  the  percent  error  decreases  with  basement  age,  but  that  if  the 
basement  age  is  a constant  10%  in  error,  the  thickness  may  be  in  error 
by  only  5-10%  for  4000  in  < — 5000  m. 

DSDP  hole  75  is  an  example  in  which  the  thickness  computed  by 
ESP  has  a large  uncertainty  because  of  basement  age  uncertainty. 
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To  find  the  uncertainty  in  thickness  from  Table  4.3,  we  first  note 

that,  according  to  Fig.  8,  ESP  assumed  a sedimentation  rate  for  hole 

75  of  about  10  m/m.y.  above  the  CCD.  From  Fig.  9,  we  note  that  ESP 

took  z to  be  about  3700  m.  Now,  van  Andel  and  Bukry  (1973)  give 

hole  75  confidence  limits  of  ±4  m.y.  Therefore,  with  a basement  age 

of  37.5  m.y.  (Table  I),  hole  75  has  a thickness  uncertainty,  according 

to  Table  4.3,  of  ±20%.  In  meters,  this  amounts  to  ±22  m:  We  must  halve 

the  meters  of  error  because  S for  hole  75  is  half  that  used  in 

o 

calculating  Table  4.3. 

Effect  of  error  in  pole  of  rotation 

In  model  ESP,  the  effect  of  error  in  the  pole  of  rotation  depends 
on  the  relative  error  in  the  resulting  tracks.  That  is,  the  effect 
depends  on  the  relative  displacement. 


A?S  £«i (Z  (1) 5 **  ’ X* ,W*  ’ (1)  ' Xo  (1)  ’ Cb  (1)  ’ £«(* ,£X*,£UJ*) 


VZ(2)  ;«Kx*,UJ'Srf0(2)  ’Xj(2)  ,tb(2)’e^*’£X*,ecu*) 


SKI)  " %(2)’ 

where  eq  (5)  defines  the  absolute  error  and  the  subscripts  (1)  and 
distinguish  one  track  from  another.  To  be  meaningful,  we  must  gauge 
the  relative  displacement  along  a given  contour  of  sedimentation  rate; 
depths  zq)  and  z(2)  are  t*ie  depths  that  the  contour  cuts  the  two 
tracks . 

Fig.  4.1  shows  an  example  of  relative  error. 


A 
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To  demonstrate  the  significance  of  relative  error  to  ESP,  we 
present  Figs.  4.2  (a)-(c).  In  these  examples,  track  e represents 
a track  along  which  ESP  predicts  thickness.  Two  other  tracks,  r^  and 
bracket  track  e and  determine  the  contours  of  sedimentation  rate. 
This  dichotomy  of  tracks  is  an  idealization  when  ESP  predicts  thickness 
along  the  same  ROPH  tracks  used  in  contouring.  However,  even  in  this 
case,  smoothing  tends  to  decouple  a track  from  the  contours  and 
Figs.  4.2  (a)-(c)  are  a good  analogy.  Conversely,  without  smoothing, 

ESP's  predictions  along  ROPH  tracks  are  virtually  error-free  (i.e.,  the 
input  is  the  prediction) . 

Fig.  4.2  (a)  indicates  that  small  relative  error  in  tracks 
results  in  small  error  in  ESP's  thickness.  In  general,  small  relative 
displacements  occur  when  the  sites  to  be  tracked  have  similar  initial 
positions  and  similar  basement  ages.  DSDP  holes  70-74  and  77  fall 
in  this  category.  In  addition,  any  tracks  that  appear  parallel 
will  probably  have  small  relative  displacements.  DSDP  holes  159 
and  160  and  161  and  78  are  examples.  As  a consequence  of  the  small 
relative  error  in  tracks,  ESP's  thickness  for  most  of  these  holes 
is  insensitive  to  error  in  the  pole  of  rotation.  However,  this 
statement  assumes  that  the  smoothed  contours  are  not  affected  by  those 
tracks  with  large  relative  errors.  This  is  probably  a valid  assumption 
for  holes  73-75,  and  78  and  161. 

Figs.  4.2  (b)  and  (c)  indicate  the  effect  of  large  relative  error 
in  tracks.  For  such  track  error,  ESP's  predictions  of  thickness  can 
be  in  error  by  a large  amount.  This  presents  us  with  the  possibility 
that  the  differences  in  Table  III  are  due  to  error  in  pole  of  rotation; 
or,  in  other  words,  that  Fig.  2 represents  Fig.  8 poorly  because  we  have 
chosen  the  wrong  pole. 


4~  A 
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Fig.  4.2  (a),  (b) , and  (c) . Effect  of  an  error  in  pole  of  rotation  on 
model  ESP.  Tracks  labelled  r^  and  are  ROPH  tracks  and  determine  contours 
of  sedimentation  rate,  S.  Track  labelled  e is  track  along  which  ESP  attempts 
to  predict  thickness;  in  these  hypothetical  examples,  it  is  merely  a track, 
without  sedimentation  rates  of  its  own.  The  contours  of  sedimentation  rate 
begin  at  a rate  S^  and  have  a contour  interval  of  As.  Solid  lines  indicate 
the  tracks  and  contours  found  with  the  true  pole  of  rotation.  Dashed  lines 
indicate  the  tracks  and  contours  found  with  an  erroneous  pole  of  rotation. 

(a)  In  this  diagram,  the  erroneous  pole  of  rotation  and  the  initial 

positions  of  the  sites  are  such  that  there  is  little  relative  displacement 

between  the  three  tracks.  There  is  little  error  in  ESP's  prediction  of 
thickness  along  track  e. 

(b)  In  this  diagram,  the  erroneous  pole  of  rotation  and  the  initial 

positions  of  the  sites  are  such  that  there  is  much  relative  displacement. 

There  is  a large  error  in  ESP's  prediction  of  thickness  along  track  e. 

(c)  In  this  diagram,  the  erroneous  pole  of  rotation  and  the  initial 
position  of  the  sites  are  such  that  the  tracks  r ^ and  r 2 cross  with  the 
erroneous  pole.  This  alters  the  contours  drastically.  Consequently,  there 
is  a large  error  in  ESP's  prediction  of  thickness  along  track  e. 


- ’*<1 
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contour  interval2 AS 
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Table  4.4  is  a test  of  this  possibility.  If  a better  pole  exists 
(within  the  limits  of  our  search) , all  four  intersecting  hole-pairs  in 
Table  4.4  should  show  a common  pole  in  the  "better  pole"  column.  They  don't 


Summary 

1)  Track  elevations,  z - z,  may  be  in  error  by  an  amount,  £ /k, 

° D 

where  6 is  the  error  in  basement  age,  zq  is  the  initial  depth  of 
b 

a given  track  and  k « 40  m.y. 

2)  Track  latitudes  may  be  in  error  by  an  amount  shown  in 
the  RANGE  rows  of  Table  4.2. 

3)  Predicted  thicknesses  are  uncertain  to  within  5-20%  because 
of  uncertainties  of  basement  age. 

4)  Predicted  thicknesses  for  certain  tracks  are  sensitive  to 
choice  of  the  pole  of  rotation,  but  no  change  in  the  pole  relieves 
the  contortions  near  the  equator  in  Fig.  2. 


*- 
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Table  4.1  Typical  values  of  the  partial  derivatives  of  track 
latitude,  <5.  This  table  indicates  the  sensitivity  of  tracks  to  changes 
in  the  choice  of  the  pole  of  rotation.  We  computed  these  partial 


derivatives  by  differentiating  eq  (4)  and  using  a pole  of  rotation 
of  72°N,  83°W,  .81°/m.y.  The  values  are  typical  for  poles  within 
10°  of  latitude  and  .2°/m.y.  of  the  one  we  used. 


The  derivative,  — , has  been  divided  by  cos  «S*  so  that  it  can 
dX* 

0 J, 

be  compared  to  on  a distance  to  distance  basis.  Thus,  the  units  of 

M ^ J 

and  are  degrees  of  latitude  per  degree  of  latitude.  The  units 
of  are  degrees  of  latitude  per  degree  of  rotation  per  m.y. 
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DERIVATIVES 
OP  HOLE  LATITUDE 

KITH 

RESPECT  TIME  (HI  B.P.) 


TO 

POLE 

5 

10 

15 

20 

25 

30 

35 

40 

HOLE 

LATITUDE 

0.08 

0.15 

0.22 

0.29 

0.35 

0.42 

0.48 

0.54 

65 

LONGITUDE 

-0.  19 

-0.20 

-0.20 

-0.21 

-0.23 

-0.25 

-0.28 

-0.31 

RATE 

-1.55 

-3.08 

-4.59 

-6.06 

-7.45 

-8.76 

-9.95- 

11.02 

LATITUDE 

0.05 

0.  11 

0.  18 

0.24 

0.30 

0.36 

0.42 

0.48 

66 

LONGITUDE 

-0.19 

-0.21 

-0.23 

-0.25 

-0.28 

-0.32 

-0.36 

-0.40 

RATE 

-1.52 

-2.99 

-4.40 

-5.72 

-6.93 

-8.02 

-8.97 

-9.75 

LATITUDE 

-0.01 

0.01 

0.03 

0.04 

0.05 

0.06 

0.06 

0^06 

68 

LONGITUDE 

-0.08 

-0.  15 

-0.22 

-0.29 

-0.36 

-0.43 

-0.50 

-0.57 

RATE 

-0.97 

-0.76 

-0.86 

-0.77 

-0.49 

-0.02 

0.64 

1. 48 

LATITUDE 

0.02 

0.09 

0.14 

0.20 

0.26 

0.31 

0.36 

0.40 

69 

LONGITUDE 

-0.15 

-0.18 

-0.21 

-0.25 

-0.30 

-0.34 

-0.39 

-0.45 

RATE 

-1.41 

-2.73 

-3.95 

-5.04 

-5.98 

-6.77 

-7.38 

-7.81 

LATITUDE 

0.00 

0.06 

0.11 

0.16 

0.20 

0.24 

0.28 

0.32 

70 

LONGITUDE 

-0.13 

-0.  17 

-0.22 

-0.27 

-0.33 

-0.38 

-0.44 

-0.50 

RATE 

-1.24 

-2.36 

-3.33 

-4.15 

-4.80 

-5.26 

-5.52 

-5.58 

LATITUDE 

-0.00 

0.05 

0.  10 

0.15 

0.  19 

0.24 

0.27 

0.31 

71 

LONGITUDE 

-0.  14 

-0.  18 

-0.23 

-0.28 

-0.34 

-0.39 

-0.45 

-0.51 

RATE 

-1.24 

-2.  36 

-3.  33 

-4.  14 

-4.78 

-5.23 

-5.48 

-5.53 

LATITUDE 

-0.02 

0.03 

0.08 

0.13 

0.  17 

0.21 

0.25 

0.28 

72 

LONGITUDE 

-0.  16 

-0.  20 

-0.  25 

-0.31 

-0.  36 

-0.42 

-0.48 

-0.54 

RATE 

-1.22 

-2.30 

-3.24 

-4.00 

-4.59 

-4.98 

-5.  17 

-5.  15 

LATITUDE 

-0.03 

0.02 

0.06 

0.  1 1 

0.15 

0.19 

0.22 

0.25 

73 

LONGITUDE 

-0.17 

-0.22 

-0.27 

-0.3  2 

-0.38 

-0.44 

-0.50 

-0.56 

RATE 

-1.  19 

-2.25 

-3.  15 

-3.88 

-4.42 

-4.76 

-4.89 

-4.81 

} 

! 

! • 
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WITH 


RESPECT 

TII1E  (MY  8 

i.  P. ) 

TO 

POLE 

5 10 

15 

20 

25 

30 

35 

40 

HOLE 

LATITUDE 

-0.10  -0.06 

-0.01 

0.03 

0.07 

0.  10 

0.13 

0.  15 

75 

LONGITUDE 

-0.23  -0.28 

-0.  34 

-0.39 

-0.45 

-0.52 

-0.58 

-0.65 

RATE 

-1.14  -2.13 

-2.95 

-3.57 

-4.00 

-4.21 

-4.  19 

-3.93 

LATITUDE 

-0.03  0.02 

0.06 

0.  10 

0e  14 

0.17 

0.21 

0.23 

77 

LONGITUDE 

-0.15  -0.20 

-0.25 

-0.31 

-0.37 

-0.43 

-0.50 

-0.56 

RATE 

-1.  12  -2.09 

-2.89 

-3.52 

-3.96 

-4.19 

-4.22 

-4.03 

LATITUDE 

o 

• 

o 

o 

• 

0 

1 

0.07 

0.  1 1 

0.  14 

0.17 

0.20 

0.22 

78 

LONGITUDE 

-0.11  -0.17 

-0.22 

-0.28 

-0.35 

-0.41 

-0.48 

-0.  54 

RATE 

-1.01  -1.86 

-2.53 

-3.03 

-3.34 

-3.45 

-3.36 

-3.05 

LATITUDE 

1 

o 

• 

0 
*= 

1 

o 

• 

o 

o 

0.03 

0.06 

0.09 

0.11 

0.13 

0.  14 

79 

LONGITUDE 

-0. 12  -0.  18 

-0.24 

-0.  31 

-0.37 

-0.44 

-0.51 

-0.58 

RATE 

-0.88  -1.59 

-2.11 

-2.44 

-2.56 

-2.48 

-2.  18 

-1.67 

LATITUDE 

-0.05  -0.02 

0.02 

0.05 

0.07 

0.09 

0.11 

0.12 

80 

LONGITUDE 

-0. 14  -0. 19 

-0.26 

-0.32 

-0.39 

-0.45 

-0.52 

-0.S9 

RATE 

-0.88  -1.58 

-2.  10 

-2.41 

-2.52 

-2.41 

-2.09 

-1.55 

LATITUDE 

-0.06  -0.03 

-0.02 

-0.00 

0.01 

0.01 

0.02 

0.01 

82 

LONGITUDE 

-0. 10  -0.  17 

-0.24 

-0.31 

-0.38 

-0.45 

-0.52 

-0.59 

RATE 

-0.S3  -0.87 

-1.00 

-0.93 

-0.65 

-0.  17 

0.53 

1.43 

LATITUDE 

-0.00  0.04 

0.07 

0.  10 

0.13 

0.16 

0.18 

0.20 

159 

LONGITUDE 

-0.10  -0. 16 

-0.22 

-0.28 

-0.34 

-0.4  1 

-0.48 

-0.54 

RATE 

-0.90  -1.64 

-2.20 

-2.58 

-2.77 

-2.77 

-2.56 

-2.15 

LATITUDE 

0.01  0.06 

0.10 

0.15 

0.  19 

0.22 

0.26 

0.29 

160 

LONGITUDE 

-0. 10  -0.  15 

-0.21 

-0.26 

-0.  32 

-0.38 

-0.44 

-0.51 

RATE 

-1.13  -2.13 

-2.97 

-3.64 

-4. 14 

-4.45 

-4.56 

-4.48 

j 


D 

i E R I 

VAT 

IVES 

( C 0 N T D ) 

WITH 

RESPECT 

TO 

POLE 

5 

10 

15 

TIRE  ((IV  B 
20  25 

.P.) 

30  35  40 

HOLE 

LATITUDE 

0.03 

0.09 

0.  14 

0.19  0.24 

0.29  0.33  0.37 

162 

LONGITUDE 

-0.10 

-0.  14 

-0.18 

-0.23  -0.28 

-0.33  -0.39  -0.44 

RATE 

-1.30 

-2,  49 

-3.55 

-4.47  -5.23 

-5.82  -6.24  -6.46 

LATITUDE 

0.03 

0.09 

0.15 

0.20  0.26 

0.31  0.36  0.40 

163 

LONGITUDE 

-0.  12 

-0.  15 

-0.  19 

-0.23  -0.28 

-0.32  -0.37  -0.43 

RATE 

-1.38 

-2.67 

-3.84 

-4.89  -5.79 

-6.53  -7.09  -7.47 

LATITUDE 

0.05 

0.  12 

0.18 

0.24  0.30 

0.36  0.42  0.47 

164 

LONGITUDE 

-0.  12 

-0.14 

-0.16 

-0.19  -0.22 

-0.26  -0.30  -0.35 

RATE 

-1.49 

-2.92 

-4.28 

-5.53  -6.68 

-7.69  -8.55  -9.25 

LATITUDE 

0.08 

0.15 

0.21 

0.28  0.34 

0.41  0.47  0.53 

166 

LONGITUDE 

-0.20 

-0.20 

-0.21 

-0.22  -0.24 

-0.27  -0.29  -0.33 

SATE 

-1.54 

-3.08 

-4.57 

-6.01  -7.38 

-8.65  -9.81-10.83 

LATITUDE 

0.08 

0.  14 

0.21 

0.28  0.34 

0.40  0.47  0.53 

168 

LONGITUDE 

-0.  14 

-0.15 

-0.16 

-0. 17  -0. 19 

-0.22  -0.25  -0.28 

RATE 

-1.54 

-3.07 

-4.55 

-5.98  -7.32 

-8.57  -9.70-10.70 

; 

j 


«U 

ft 

M. ' 

1 


Table  4.2.  Ranges  of  track  latitude  produced  by  five  published  poles 
of  rotation.  This  table  indicates  some  of  the  uncertainty  in  back-tracked 
paleolatitudes  in  the  central  equatorial  Pacific. 

Rows  labelled  MAX  show  the  maximum  latitude  at  a given  time  (m.y.b.p.) 
among  the  five  tracks.  Rows  labelled  MIN  show  the  minimum  latitud.  among 
the  five  tracks.  Rows  labelled  RANGE  show  MAX  minus  MIN.  Positive  values 
indicate  northern  latitude,  negative  values  southern  latitude  for  MAX  and 
MIN  entries.  The  range  is  the  absolute  value  of  the  difference. 
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I A T I T U D E 


range 


AGE  (MlO 


5 

10 

15 

20 

25 

30 

35 

HOLE 

MAX 

2.96 

1.57 

0.19 

-1.17 

-2.52 

-3.83 

-5.10 

65 

MIN 

2.69 

0.92 

-0.96 

-2*92 

-4.95 

-7.03 

-9.15 

RANGE 

0.27 

0.66 

1.15 

l.?5 

2.43 

3.21 

4.06 

MAX 

1.02 

-0.34 

-1.66 

-p.95 

-4.20 

-5.40 

-6.54 

66 

MIN 

0.50 

-1.47 

-3.50 

-S.59 

-7.72 

-9.86 

-12.01 

RANGE 

0.51 

1.13 

1.84 

p.64 

3.52 

4.47 

5.47 

MAX 

16.26 

15.89 

15.62 

15.44 

15.35 

15.35 

15.45 

68 

MIN 

14.92 

13.23 

11.65 

10*19 

8.88 

7.71 

6.69 

RANGE 

1.34 

2.67 

3.97 

?. 24 

6.47 

7.65 

8.77 

MAX 

4.71 

3.47 

2.27 

1.13 

0.06 

-0.95 

-1.88 

69 

MIN 

3.93 

1 .81 

-0.33 

-2*47 

-4.61 

-6.73 

-8.80 

RANGE 

0.79 

1.66 

2.60 

3.61 

4.67 

5.78 

6.92 

max 

5.20 

4.12 

3.10 

p.  16 

1.30 

0.52 

-0.17 

70 

MJN 

4.19 

2.05 

-o.on 

-?.  18 

-4 .23 

-6.22 

-8.12 

range 

1.00 

2.07 

3.19 

4.34 

5.53 

6.73 

7.94 

MAX 

3.33 

2.25 

1.24 

0.30 

-0.56 

-1.34 

-2.02 

71 

MIN 

2.32 

0.18 

-1.95 

-4.05 

-6.10 

-8.08 

-9.98 

RANGE 

1.00 

2.07 

3.19 

4.35 

5.54 

6.74 

7.96 

Max 

-0.68 

-1.74 

-2.72 

-3.63 

-4,46 

-5.20 

-5.84 

72 

MIN 

-1.71 

-3.85 

-5.97 

-8.06 

-10.09 

-12.05 

-13.92 

RANGE 

1.03 

2.11 

3.25 

4.43 

5.63 

6.85 

8.08 

Max 

-3.01 

-4.04 

-5.00 

-5.88 

-6.68 

-7.38 

-7.99 

73 

min 

-4.06 

-8.31 

-10.38 

-12.40 

-14.33 

-16.17 

range 

1.05 

2.15 

3.31 

4.50 

5.72 

6.95 

6.18 

MAX 

-7.31 

-8.32 

-9.25 

-10.10 

-10.86 

-11.52 

-12.09 

74 

MIN 

-8.38 

-10.51 

-12.62 

-14.67 

-16.66 

-18.57 

-20.37 

RANGE 

1.07 

2.19 

3.37 

4.57 

5.81 

7.05 

8.26 

MAX 

-13.57 

-14.55 

-15. AS 

—16.26 

-16. 9e 

-17.60 

-18.12 

75 

MIN 

-14.66 

-16.78 

-18.87 

-20.91 

-22.88 

-24.75 

-26.52 

range 

1.09 

2.24 

3.43 

4.65 

5.90 

7.15 

8.40 

MAX 

-0.55 

-1.52 

-2.40 

-3.20 

-3.91 

-4.53 

-5.06 

77 

MIN 

-1.66 

-3.78 

-5,86 

-7.89 

-9.8.4 

-11.70 

-13.45 

range 

1.11 

2.26 

3,46 

4.68 

5.92 

7.16 

8.39 

MAX 

7.01 

6.16 

5.38 

4.69 

4.09 

3.58 

3.17 

78 

MIN 

5.84 

3.76 

1.75 

-0.20 

-2.05 

-3.80 

-5.44 

RANGE 

1.18 

2.39 

3.63 

4.89 

6.14 

7.38 

8.60 

40 


-6.32 
-11.30 
A.  98 

-7.61 

-14.14 

6.53 

15.65 

5.84 

9.81 

-2.74 

-10.82 

8.08 

-0.77 

-9.92 

9.15 

-2.62 

-11.78 

9.16 

-6.39 

-15.68 

9.29 

-8.50 

-17.89 

9.39 

-12.56 

-22.06 

9.50 

-18.53 
-28.  15 

9.62 

-5.98 

-15.08 

9.60 
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LATITUDE  range  ( C ONTO) 


AGE(mYBP) 


5 

10 

15 

20 

25 

30 

35 

60 

HOLE 

MAX 

1.72 

0.98 

0.33 

-0.23 

-0.69 

79 

MIN 

0.68 

-1.52 

-3.65 

-5.29 

-7.02 

RANGE 

1.26 

2.50 

3. 78 

5.06 

6.33 

Max 

-1.7P 

-2.5? 

-3.17 

-3.72 

-6.18 

80 

MIN 

-3.03 

-5.03 

-6.96 

-R . 79 

-10.51 

RANGE 

1 . 2 A 

2.51 

3.79 

5*07 

6.36 

MAX 

2.07 

1.65 

82 

MIN 

0.76 

-1.01 

RANGE 

1.33 

2.66 

MAX 

1 1 .69 

10.73 

10.05 

9.66 

8.96 

159 

MIN 

10.26 

e.26 

6.29 

6.63 

2.67 

RANGE 

1.23 

2.69 

3.76 

5.03 

6.29 

max 

10.67 

9.70 

e.79 

7.96 

7.22 

6.56 

5.99 

5.52 

160 

MIN 

9.5P 

7.65 

5.36 

3.33 

1.36 

-0.52 

-2.31 

-3.98 

RANGE 

1.10 

2.26 

3.63 

A.  63 

5.86 

7. OB 

8.30 

9.69 

max 

9.10 

8.03 

7.0? 

6.00 

5.22 

6.66 

3.75 

3.16 

161 

min 

8.09 

5.95 

3.8? 

1.72 

-0.32 

-2.30 

-6.20 

-6.00 

range 

1.01 

2.08 

3.20 

4.36 

5.56 

6.76 

7.95 

9.16 

MAX 

13.69 

12.55 

11.67 

1 0.A6 

9.52 

8.65 

7.87 

7.17 

16? 

MIN 

12.75 

10.60 

8.65 

6.33 

6.26 

2.21 

0.26 

-1.66 

RANGE 

0.96 

1 .95 

3.02 

4.13 

5.28 

6.65 

7.63 

8.81 

max 

9.99 

0.77 

7.51 

6.50 

5.66 

6.69 

3.59 

2.78 

163 

min 

9.15 

7.02 

6.88 

2.73 

0.60 

-1.69 

-3.54 

-5.52 

range 

0.83 

1.75 

2.73 

3.77 

6.66 

5.98 

7.13 

8.30 

max 

11.85 

10.5? 

9.23 

7.98 

6.78 

5.63 

4.55 

3.53 

1 66 

min 

11.26 

9.20 

7.1? 

4.99 

2.85 

0.70 

-1.63 

-3.56 

range 

0.61 

1.32 

2.1? 

2.99 

3.93 

6.93 

5.98 

7.07 

max 

2.37 

0.98 

-0.39 

-1.75 

-3.07 

-6.37 

-5.62 

-6.82 

166 

min 

2.06 

0.26 

-1.67 

-3.65 

-5.71 

-7.80 

-9.96 

-12.08 

range 

0.31 

0.76 

1.28 

1.91 

2.63 

3.66 

6.32 

5.26 

MAX 

9.27 

7.89 

6.52 

5-17 

3.85 

2.57 

1.36 

0.15 

168 

min 

8.92 

7.07 

5.13 

3.12 

1.06 

-1.07 

-3.21 

-5.36 

rangf 

0.35 

0.82 

1.39 

? . 05 

2.ei 

3.66 

6.56 

5.51 

K 
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Table  4.3.  Estimated  error  in  predicted  thickness  due  to  error  in 
basement  age.  The  rows  labelled  ERROR  are  the  estimated  thickness  errors 
in  meters  for  each  million-year  error  in  basement  age.  The  rows  labelled 
PERCENT  are  the  estimated  percent  thickness  errors.  The  column  labelled 
ZCCD  is  the  carbonate  compensation  depth,  zqq^-  The  errors  in  this  table 
were  calculated  from  eqs  (9)  and  (]0)  with  Sq  = 20  m/m.y.,  E = 3300  m, 
k = 40  m.y.,  and  C = 0. 

To  find  the  uncertainty  in  predicted  thickness  at  any  site,  use  the 
local  determined  from  Fig.  9,  or  use  the  local  average  sedimentation 

rate,  multiply  the  error  in  meters  from  this  table  by  the  ratio  of  the 
local  rate  and  20  m/m.y.;  this  correction  is  unnecessary  for  the  percent 
error. 
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Table  4.4  "Better  poles  of  rotation  for  four  pairs  of  inter- 
secting DSDP  tracks.  Column  gives  the  sedimentation  rates  (m/m.y.) 
for  each  track  of  the  hole  pair  at  the  point  of  intersection  of  the 
tracks.  Column  t^  gives  the  ages  (m.y.b.p.)  of  intersection.  S ^ and 
t^  are  from  Fig.  4 and  result  from  a pole  of  rotation,  72°N,  83°W, 
•8l°/m.y. 

We  found  the  "better"  times  of  crossing,  t^ , by  displacing  by  eye 
one  track  in  each  pair  relative  to  the  other  such  that  the  resulting  local 
sedimentation  pattern  agreed  more  with  the  general  pattern  of  Fig.  8 and 
such  that  the  S were  more  nearly  equal.  Given  these  t , we  found  the 

X x 

"better"  poles  by  searching  the  tables  in  the  addendum  to  this  appendix. 
The  possible  "better"  poles  that  we  found  are  indicated  by  a rate  and 
position,  the  position  symbolized  as  follows: 

h = high  latitude  (73-75°N) 

1 = low  latitude  (63-65°N) 
m = middle  latitude  (66-72°N) 

E = easterly  longitude  (40-50°W) 

W = westerly  longitude  (100-110°W) 

M = middle  longitude  (60-90°W) 


r 


wide  range  of  possible  latitudes  and  longitudes 


Addendum 


This  addendum  contains  tables  of  the  variation  of  the  time  of  track 
intersections  as  a function  of  pole  of  rotation  and  basement  age  for  four 
pairs  of  DSDP  holes:  holes  70  and  79,  71  and  79,  72  and  80,  and  77  and  80 
The  addendum  also  contains  the  computer  program  which  generated  the  tables 

Numbers  in  the  body  of  the  tables  are  the  ages  (m.y.b.p.)  of  inter- 
section for  the  given  hole  pair.  For  example,  with  a pole  rate  of  .6°/m.y 
pole  latitude  of  65°N,  and  pole  longitude  of  90°W  (-90),  the  hole  pair  70, 
79  will  intersect  at  three  possible  times,  depending  on  basement  age  error 

21,0 

24,3 

27,7. 

Ages  to  the  left  of  the  commas  refer  to  ages  along  the  track  of  the 
hole  listed  first  (70,  in  this  case);  ages  to  the  right  refer  to  the  hole 
listed  second  (79,  in  this  case). 

In  the  first  row  of  the  triplet,  the  intersections  are  for  tracks 
computed  with  a basement  age  10%  high  for  the  first  hole  and  10%  low  for 
the  second  hole.  In  the  second  row  of  the  triplet,  the  intersections  are 
for  tracks  computed  with  no  error,  using  basement  ages  shown  in  Table  I. 

In  the  third  row  of  the  triplet,  the  intersections  are  for  tracks  computed 
with  a basement  age  10%  low  for  the  first  hole  and  10%  high  for  the  second 
hole. 

When  tracks  do  not  intersect,  the  tables  indicate  this  as:  **,**. 
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VARIATION  OF  TRACK  INTERSECTIONS 
AS  A FUNCTION  OF  POLE  POSITION  AND  RATE 
AND  OF  ERROR  IN  BASEMENT  AGE 

HOLE  70  LAT  6.34  LONG  -140.36  DEPTH  5059 

HOLE  79  LAT  2.55  LONG  -121.57  DEPTH  4566 


RATE=  0.60 


LONG 

-110 

< 

o 

o 

i 

V^> 

O 

-80 

-70 

-60 

-50 

-40 

LAT 

♦*  #** 

35,  9 

21,  0 

**f  ** 

**,  ** 

** 

65 

39,14 

24,  3 

**,*» 

** 

43,20 

27,  7 

**,  *# 

**# 

♦ * # ** 

44,14 

27,  4 

**,  ** 

** 

** 

67 

**  f ** 

48,20 

30,  7 

**»  ** 

** 

**  f ** 

**.** 

**  # ** 

52,26 

35,13 

23,  3 

**,  ** 

*♦ 

**f  ** 

**  ,** 

56,20 

33,  8 

24,  2 

**,  ** 

** 

**,*<■ 

69 

37,13 

26,  4 

**,  *♦ 

**# 

♦ * 

«■*,** 

42,20 

31,  10 

23,  4 

*♦ 

**  r** 

**,*♦ 

41,13 

29,  6 

23,  1 

** 

** , ** 

**,*♦ 

71 

**#  ** 

46,  19 

33,  10 

25,  4 

20, 

0 

**  #** 

** 

53,26 

40,  18 

31,  11 

26, 

6 

24,  4 

24,  4 

♦*,** 

** 

52,19 

38,11 

30,  6 

26, 

4 

24,  2 

23,  2 

73 

**  # ♦ 

58,24 

43,  17 

34,  11 

29, 

7 

28,  6 

27,  6 

**#  ** 

**f  ** 

**  ,** 

51,25 

43,20 

37, 

15 

35,14 

37, 15 

++,** 

49,17 

39,  12 

34, 

9 

31,  7 

32,  7 

75 

♦ ♦ # ** 

**#  *♦ 

**  ,** 

54,23 

45,18 

40, 

14 

37,13 

38,13 

**  ,** 

**,  ** 

53,  26 

49, 

24 

48,  24 

48,24 

BASEHENT 
AGES  USED 

67,22 
61,  25 
55,27 


* 


* 

! > 


*• 
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INTERSECTIONS  POR  HOLES  70/79,  CONT*D 


RATE=  0.80 


LONG 

•1 10 

-100 

-90 

-80 

-70 

-60 

LAT 

**  § ** 

**f  *♦ 

•*#  ** 

**f 

** 

**# 

** 

**# 

65 

22,  1 

**# 

** 

**# 

♦ * 

**# 

25,  6 

**(  *• 

**r 

** 

**# 

** 

**# 

***** 

26,  3 

**,  ** 

**# 

** 

**# 

*♦ 

**# 

67 

29,  7 

♦ *,** 

** 

**# 

♦ ♦ 

**# 

** 

35,13 

*«■  f *# 

**• 

*♦ 

**# 

*♦ 

**# 

62,22 

***** 

**# 

** 

**• 

** 

**f 

69 

58,24 

*** 

** 

**# 

** 

*4# 

**f  *♦ 

53,27 

21,  2 

**# 

*♦ 

*** 

** 

**r 

**  #** 

**,♦* 

23,  1 

**# 

** 

**# 

** 

**# 

71 

27,  5 

**, 

** 

**# 

** 

**# 

**  r ** 

**t+* 

32,11 

**r 

** 

**# 

** 

*** 

33,  8 

**r 

** 

**# 

** 

**# 

73 

*♦,  ** 

40,15 

22, 

1 

**# 

♦ * 

**# 

**  9 ** 

**,  ** 

51,26 

27, 

7 

**# 

** 

**r 

62,22 

29, 

6 

**# 

♦ ♦ 

**r 

75 

♦ * ,*♦ 

• * ,** 

35, 

11 

23, 

2 

**# 

**  #*4 

♦ *,  «■* 

+*t++ 

45, 

22 

30, 

9 

22, 

3 

BASEKENT 
ACES  USED 

67,22 

61,25 

55,27 
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INTERSECTIONS  FOR  HOLES  70/79,  CONT'D 


LONG 

-no 

-100 

90 

RATE* 

-80 

1.00 

-70 

-60 

* 

50 

• 

© 

LAT 

* 

#** 

61f  22 

** 

t *♦ 

** 

** 

+ * 

* 

** 

65 

♦ 

#** 

**,  ** 

,** 

♦ * 

** 

** 

• 

** 

* 

, ** 

**f  ** 

** 

# ** 

** 

**# 

** 

♦ *p 

** 

♦ 

** 

* 

56,  20 

** 

** 

** 

**# 

♦ * 

* 

** 

67 

♦ 

56,24 

*♦ 

r ** 

** 

♦ * 

** 

* 

** 

**, 

♦ 

9 ** 

**#  ** 

** 

# ** 

*♦ 

** 

** 

* 

* * 

21,  0 

*♦ 

** 

♦ ♦ 

*♦ 

♦ 

** 

**, 

69 

♦ 

# ** 

**,  *♦ 

** 

** 

**  9 

** 

**  9 

** 

* 

** 

♦ 

#** 

** 

# ** 

*♦ 

** 

** 

* 

*♦ 

**# 

♦ 

f ♦* 

**#  ** 

** 

,** 

** 

** 

** 

* 

** 

**, 

71 

* 

,♦* 

♦ * # * * 

♦ * 

* ** 

** 

♦ * 

** 

♦ 

*♦ 

**, 

* 

** 

♦ * 

**  9 

** 

*♦ 

* 

** 

**, 

* 

**r  ♦* 

** 

# $t 

** 

**9 

♦ * 

**# 

** 

♦ 

** 

**, 

73 

♦ 

**#  ** 

** 

, ♦* 

** 

+ + 9 

** 

** 

* 

** 

**  9 

• 

,** 

27 

, 7 

** 

+ + 9 

** 

* * 

* 

** 

**( 

* 

*♦ 

# ** 

*♦ 

**  0 

** 

** 

* 

** 

7S 

♦ 

** 

r ** 

** 

**9 

** 

** 

♦ 

** 

♦ 

**#  ** 

** 

f ** 

19 

0 

** 

** 

* 

** 

**r 

BASEMENT 
AGES  USED 


67,  22 
61,25 
55,27 
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VARIATION  OF  TRACK  INTERSECTIONS 
AS  A FUNCTION  OF  POLE  POSITION  AND  RATE 
AND  OF  ERROR  IN  BASEMENT  AGE 

HOLE  71  LAT  4.47  LONG  -140.31  DEPTH  4419 

HOLE  79  LAT  2.55  LONG  -121.57  DEPTH  4566 


RATE=  0.60 


LONG 

-1  10 

-100 

-90 

0 

GO 

1 

-70 

-60 

-50 

-40 

LAT 

24 

,15 

21,  14 

19,13 

19,  13 

19,  13 

19,  13 

20,13 

21,14 

65 

24 

.17 

22,  16 

21,15 

20,15 

21,15 

21.15 

22,16 

24,17 

24 

.19 

23,  18 

22,  17 

22,  17 

23,  18 

24,  18 

25,  19 

27,20 

26 

,16 

23,  14 

21,14 

20.13 

20,13 

20,14 

21,14 

23,14 

67 

27 

,18 

24,16 

22,  16 

22,  16 

22,  16 

23,  16 

23,  16 

25,17 

27 

.20 

25,  19 

24,18 

24,18 

24,19 

25,  19 

26.20 

20,21 

30 

,17 

25,15 

23,  14 

22,  14 

22,  14 

22,  14 

22,14 

24,15 

69 

30 

.19 

26,18 

24,17 

23,  16 

24,16 

24,17 

25,17 

26,  18 

30 

,21 

27,20 

26,  19 

25,  19 

26.  19 

27,  20 

28,20 

30,21 

34 

.19 

28,  16 

25,  15 

23,15 

23,15 

23,  15 

24,15 

25,  15 

71 

34 

.21 

28,  18 

26,18 

25,  17 

25,  17 

26,  17 

26,  18 

28,  18 

34 

.23 

30,21 

28,20 

27,20 

28,20 

28,21 

29,21 

31,22 

42 

.21 

31,17 

27,  16 

26,  16 

25,  15 

25,  15 

26.  16 

27,16 

73 

42 

,29 

32,20 

29,  19 

28, 18 

28,18 

28,  18 

28,18 

30,19 

42 

.26 

33,  23 

31,22 

30,21 

30,21 

30,  22 

31,22 

33,23 

** 

,** 

36,19 

31,18 

29,17 

28,16 

28,16 

28.16 

29,17 

75 

** 

,** 

37,22 

33,20 

31,  19 

30,  19 

30,  19 

31,19 

32,20 

+* 

.•* 

39,25 

35,23 

33,23 

33,23 

33,23 

34,23 

36,24 

BASEMENT 
AGES  USED 

67,22 
61,  25 
55,27 
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INTERSECTIONS  FOR  HOLES  71/79,  CONT'D 


Rf.TE=  0.80 


LONG 

-1 10 

-100 

-90 

-80 

-70 

-60 

-50 

-40 

LAT 

18,12 

16,12 

16,11 

15,11 

16,11 

16, 

12 

17,12 

18,12 

65 

18,14 

17,13 

17,13 

17,  13 

17,  13 

18, 

14 

19,  14 

20,15 

18, 15 

18,15 

18,15 

18,15 

19,16 

20, 

16 

21,17 

23,18 

20,13 

18,  12 

17,  12 

16,  12 

16,  12 

17, 

12 

18,  12 

19,13 

67 

20,15 

18,  14 

18,  14 

18,14 

18,14 

19, 

14 

20,15 

21,15 

20,16 

19,16 

19,16 

19,  16 

20,  16 

21. 

17 

22,  17 

24,18 

22,14 

19,13 

18,13 

18,12 

17,12 

18, 

12 

19,13 

20,13 

69 

23,16 

20,15 

19,  14 

19,  14 

19,  14 

20, 

15 

21,15 

22,16 

22,  17 

21,17 

21,17 

21,17 

21,17 

22, 

17 

23,18 

25,19 

26 ,16 

21,14 

20.13 

19,  13 

19,  13 

19, 

13 

20,13 

21,14 

71 

26,17 

22,16 

21,15 

20,15 

20,15 

21, 

15 

22,16 

23,  16 

26 ,19 

23,  18 

23,  18 

22,  17 

23,  18 

24, 

18 

25,19 

26,19 

51,22 

24,15 

22,14 

21,14 

20,13 

21. 

14 

21,14 

22,14 

73 

31  ,20 

25,17 

23,  16 

22,  16 

22,  16 

23, 

16 

23,  16 

25,17 

30,21 

26,19 

24,19 

24,18 

24,19 

25, 

19 

26,20 

28,20 

**  ,** 

28,  16 

24,15 

23,  15 

22,  14 

23, 

14 

23,15 

24,15 

75 

#*,** 

29,19 

26,  17 

25,  17 

24,17 

25, 

17 

25,17 

26,18 

**  ,*♦ 

30,21 

27,20 

27,  20 

27,20 

28, 

20 

28,21 

30,21 

BASEMENT 
AGES  USED 

67,22 
61,  25 
55,27 


196 


INTERSECTIONS  FOR  HOLES  71/79,  CONT'D 


RATE3  1.00 


LONG 

-110 

-100 

-90 

© 

GO 

1 

-70 

-60 

0 

uO 

1 

0 

1 

LAT 

14,11 

14,10 

13,10 

13,10 

14,10 

14,11 

15,11 

16,  11 

65 

15,12 

14,  12 

14,  12 

14,  12 

15,  12 

16,  13 

17,  13 

18,14 

15,13 

15,13 

15,  14 

16,  14 

17,14 

18,15 

19,16 

20,  16 

16,12 

15,11 

14,11 

14,  11 

14,  11 

15,  11 

15,11 

16,12 

67 

16,13 

15,  12 

15, 12 

15,  12 

16,13 

16,13 

17,13 

18,14 

16 ,14 

16,14 

16,  14 

17,  14 

17,  15 

18,  15 

20,16 

21,17 

51,22 

16,  12 

15,11 

15,11 

15,11 

16,11 

16,12 

17,12 

69 

18,14 

17,13 

16,  13 

16,  13 

17,  13 

17,  13 

18,14 

19,14 

18,15 

17,15 

17,15 

18,15 

18,15 

19,  16 

21,17 

22,  17 

**  ,** 

18,12 

16,  12 

16,  12 

16,  12 

16,  12 

17,12 

18,12 

71 

49,24 

18,14 

17,13 

18.13 

18,14 

18,  14 

19,14 

20,  15 

20,16 

19,16 

19,15 

19,  16 

20,  16 

20,  16 

22,17 

23,18 

♦*  ,** 

20,13 

18,13 

18,12 

17,  12 

18,  12 

18,13 

19,  13 

73 

** , *♦ 

21,15 

19,14 

19,14 

19,14 

20,14 

20,15 

21.15 

23,18 

21,17 

20,  16 

21,  17 

21,  17 

22,  17 

23,  18 

24,19 

23,15 

21,14 

20,13 

19,13 

19,13 

20,13 

21,14 

75 

**  ,** 

24,  16 

22,  16 

21,  15 

21,  15 

21,  15 

22,16 

23,16 

** , ** 

25,  19 

23,18 

23,18 

23,18 

24,10 

25,19 

26,  19 

BASEMENT 
AGES  USED 

67,  22 
Cl, 25 
55,  27 
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VARIATION  OF  TRACK  INTERSECTIONS 
AS  A FUNCTION  OF  POLE  POSITION  AND  RATE 
AND  OF  ERROR  IN  BASEMENT  AGE 

HOLE  77  LAT  0.48  LONG  -133.23  DEPTH  4291 

HOLE  80  LAT  -0.96  LONG  -121.55  DEPTH  4399 


RATE=  0.60 


LONG 

-110 

-100 

-90 

-80 

-70 

-60 

-50 

-40 

LAT 

32 

19 

23,16 

20,14 

19,14 

18,13 

18,13 

18,13 

19,14 

65 

32 

21 

25,18 

22,  17 

21,  16 

21,  16 

21,  16 

21,16 

22,17 

31 

24 

27,22 

25,21 

24,20 

25.20 

25,21 

26,21 

27,22 

39 

20 

25,  17 

22,15 

21,  14 

20,  14 

19,  14 

19,  14 

20.14 

67 

37 

23 

27,19 

24,18 

22,17 

22,17 

22,  17 

22,17 

23,17 

35 

25 

29,23 

27,21 

26,  21 

26,  21 

26,  21 

27,22 

28,22 

*♦ 

*♦ 

28,18 

23,16 

22. 15 

21,14 

20,  14 

20,14 

21,15 

69 

** 

** 

30,21 

26,  19 

24,  18 

24,  18 

23,  18 

24,  18 

24,  18 

** 

** 

32,24 

29,23 

27,22 

27,22 

27,22 

28,22 

29.23 

** 

*♦ 

31,19 

26,17 

24,  16 

22,  15 

22,  15 

22,15 

22,15 

71 

*♦ 

33,22 

28,20 

26,19 

25,19 

25,18 

25,18 

26,  19 

♦ * 

** 

35,25 

31,24 

30,23 

29,  23 

29,  23 

29,23 

30,23 

** 

*♦ 

36,20 

28,18 

26,17 

24,16 

24,16 

23,16 

24,16 

73 

• * 

♦ * 

37,23 

31,21 

29,  20 

27,  19 

27,  19 

27,19 

28,20 

** 

** 

#*,  ** 

34,25 

32,24 

31,24 

31,24 

31,24 

32,24 

** 

** 

**f  ** 

32,  19 

29,  18 

27,  17 

26,  17 

25,  17 

26,17 

75 

** 

** 

35,22 

31,21 

30,21 

29,20 

29,20 

30,21 

** 

♦ ♦ 

** , ** 

34,25 

33,25 

33,  24 

33,24 

33,25 

BASEMENT 
AGES  USED 

46,21 
42,  23 
37,25 


198 


INTERSECTIONS  FOR  HOLES  77/30,  CONT'D 


RATE= 

0,80 

LONG 

-1 10 

-100 

-90 

» 

cc 

o 

-70 

-60 

-50 

-40 

LAT 

24 

.16 

18,13 

16,12 

15,12 

15,12 

15, 

12 

16,12 

16,12 

65 

24 

.18 

19,15 

18,  15 

17,  14 

17,  14 

IB. 

14 

18,15 

19,  15 

23 

. 19 

21,18 

20,18 

20,18 

20,18 

21, 

18 

22,19 

23,19 

** 

.** 

20,  14 

18,13 

16,  12 

16,  12 

16. 

12 

1 6,  12 

17,12 

67 

28 

.20 

21,16 

19,15 

19,15 

18,15 

19, 

15 

19,15 

20,  16 

26 

.21 

23,19 

22,  19 

21,  18 

21,  18 

22, 

19 

23,  19 

24,20 

** 

»** 

22,15 

19,  14 

18,13 

17,13 

17, 

13 

17,13 

18,13 

69 

** 

24,18 

21, 16 

20,  16 

19,  15 

20, 

16 

20,  16 

21,16 

35 

.25 

25,21 

24,20 

23,19 

23,19 

23, 

20 

24,20 

25,21 

♦ * 

.*♦ 

25,  17 

21,15 

19,  14 

18,  13 

18, 

13 

18,13 

19,13 

71 

** 

.♦* 

27,19 

23,17 

21,16 

21,16 

21, 

16 

21,16 

22,17 

** 

28,22 

26,21 

24.  20 

24,20 

25, 

20 

25,21 

26,21 

** 

29,18 

23,16 

21,14 

20,14 

20, 

14 

20,14 

20,  14 

73 

** 

31,21 

26,  19 

23,  18 

22,  17 

23, 

17 

23,17 

23,17 

** 

. ** 

33,24 

28,22 

27,21 

26,21 

26, 

21 

27,22 

28,22 

** 

39,20 

27,17 

23,  16 

22,  15 

21, 

15 

21.  15 

21,15 

75 

*♦ 

**,*♦ 

29,20 

26,19 

25,18 

24, 

18 

24.18 

25,  18 

** 

32,24 

29,  23 

29,  23 

28, 

22 

29,23 

29,23 

BASEMENT 
AGES  USED 

46,2 1 
42,  23 
37,25 
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INTE8SECII0NS 

o 

to 

HOLES 

77/80 

, CONT'D 

BATE= 

1.00 

LONG 

-1  10 

-100 

-90 

-80 

-70 

-60 

-50 

-40 

LAT 

♦ * 

15,12 

14,11 

13,11 

13,10 

14. 

11 

14,11 

14,11 

65 

37 

#23 

16,14 

15,  13 

15,  13 

15,  13 

15, 

13 

16,13 

17,14 

18 

, 16 

17,16 

17,15 

17.  16 

18,16 

19. 

17 

20,17 

21,18 

♦ * 

,** 

16,12 

15,11 

14,  11 

14,  11 

14, 

11 

14,11 

15,  1 1 

67 

*♦ 

t ** 

17,  14 

16,  14 

16,13 

16,13 

16, 

13 

17,  14 

18,  14 

20 

, 17 

19,17 

18,  16 

18,  16 

19,  17 

19, 

17 

21,18 

22,18 

*>* 

,** 

18,13 

16,12 

15,12 

15,11 

15, 

11 

15,12 

15,  12 

69 

*♦ 

19,15 

18,14 

17,  14 

17,  14 

17, 

14 

10,  14 

18,15 

35 

#25 

20,18 

20,17 

20,17 

20,17 

20. 

18 

21,18 

22,19 

** 

#** 

21,14 

18,  13 

16,  12 

16,  12 

16, 

12 

16,12 

16,12 

71 

*♦ 

,** 

22,17 

19,15 

19,15 

18,15 

18. 

15 

19,15 

19,15 

** 

#♦* 

23,  19 

21,18 

21,  18 

21,  18 

22, 

18 

23,19 

23,20 

** 

# ** 

25,16 

20,14 

18,13 

17,13 

17, 

13 

17,13 

17,13 

73 

** 

,** 

26,19 

22,17 

20,  16 

20,  15 

19. 

15 

20,16 

20,16 

** 

, ** 

27,21 

24,20 

23,  19 

23,19 

23, 

19 

24,20 

25,20 

** 

#♦* 

*♦,** 

23,  15 

20,  14 

19,  13 

19. 

13 

18,  13 

19,13 

75 

** 

# ** 

25,18 

23,17 

21,16 

21. 

16 

21,17 

22,  17 

** 

#** 

35,25 

27,22 

26,  21 

25,20 

25, 

20 

26,21 

26,21 

BASEMENT 
AGES  USED 

46,21 
42,  23 
37,25 
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VARIATION  OP  TRACK  INTERSECTIONS 
AS  A JUNCTION  OF  POLE  POSITION  AND  RATE 
AND  OF  ERROR  IN  BASEMENT  AGE 

HOLE  72  LAT  0.44  LONG  -138.87  DEPTH  4326 

HOLE  80  LAT  -0.96  LONG  -121.55  DEPTH  4399 


RATE®  0.60 

LONG  -110  -100  -90  -80  -70  -60  -50  -40 

LAT 


18,11 

15. 

10 

14, 

9 

14, 

9 

13, 

9 

14, 

9 

14, 

9 

15. 

10 

65 

18,12 

16. 

11 

15, 

11 

15, 

11 

15, 

10 

15, 

11 

16, 

11 

17. 

12 

19.14 

17. 

13 

16, 

12 

17, 

13 

17, 

13 

17. 

13 

18, 

13 

20, 

14 

19,12 

17, 

10 

15, 

10 

15, 

9 

14, 

9 

14, 

9 

15, 

9 

16, 

10 

67 

20,13 

17. 

12 

16, 

11 

16, 

11 

16, 

11 

16, 

11 

17, 

11 

18, 

12 

20,15 

18, 

14 

17. 

13 

18. 

13 

18, 

13 

18, 

13 

19, 

14 

21. 

15 

22,13 

18, 

11 

16, 

10 

16, 

10 

15, 

10 

15. 

10 

16, 

10 

17. 

10 

69 

22.  14 

19. 

13 

18, 

12 

17, 

12 

17, 

12 

17. 

12 

18. 

12 

19, 

13 

22,16 

20, 

15 

19, 

14 

19, 

14 

19, 

14 

19, 

14 

20, 

15 

22. 

16 

25,14 

20, 

12 

18, 

11 

17, 

11 

17, 

10 

17, 

10 

17, 

10 

18. 

11 

71 

25,10 

21. 

14 

19, 

13 

1.9, 

13 

18, 

12 

18. 

12 

19, 

13 

20, 

13 

25,17 

22, 

16 

21, 

15 

20, 

15 

20. 

15 

21, 

15 

22, 

15 

23, 

16 

29,16 

23. 

13 

20. 

12 

19, 

11 

18, 

11 

18, 

11 

18, 

11 

20, 

12 

73 

30,  10 

24, 

15 

21, 

14 

21, 

13 

20, 

13 

20, 

13 

21, 

13 

22, 

14 

30,20 

25, 

17 

23, 

16 

22, 

16 

22. 

16 

23. 

16 

23, 

16 

25. 

17 

36,18 

26, 

14 

23, 

13 

21, 

12 

20, 

12 

20, 

12 

20, 

12 

22. 

12 

75 

36,20 

28, 

17 

24. 

15 

23, 

15 

22, 

14 

22, 

14 

23, 

14 

24, 

15 

36,22 

29. 

19 

26, 

18 

25, 

17 

25, 

17 

25. 

17 

25, 

18 

27, 

19 

BASEHENT 
AGES  USED 


63,21 
56,  23 
52,25 


4. 

t 

A 
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INTERSECTIONS  FOR  HOLES  72/80,  COHT'D 


RATE*  0.80 

LONG  -110  -100  -90  -80  -70  -60  -50  -40 

LAT 


14,  9 

12, 

8 

11, 

8 

11. 

7 

11. 

8 

11. 

8 

11. 

6 

12. 

8 

65 

14 ,10 

13, 

9 

12. 

9 

12, 

9 

12. 

9 

12, 

9 

13, 

9 

14, 

10 

14,11 

13, 

11 

13, 

10 

13, 

10 

14, 

11 

14, 

11 

15, 

12 

16, 

12 

15,10 

13, 

9 

12, 

8 

12, 

8 

12. 

8 

12. 

8 

12, 

8 

13, 

8 

67 

15,11 

14, 

10 

13, 

9 

13, 

9 

13, 

9 

13, 

9 

14, 

10 

15, 

10 

15,12 

14, 

11 

14. 

11 

14, 

11 

14. 

11 

15. 

12 

16, 

12 

17. 

13 

16,10 

14, 

9 

13. 

9 

12, 

8 

12. 

8 

12, 

8 

13, 

0 

14. 

9 

69 

17,12 

15, 

11 

14, 

10 

13, 

10 

14, 

10 

14, 

10 

14, 

10 

15, 

1 1 

17,13 

16, 

12 

15, 

12 

15, 

11 

15, 

12 

16, 

12 

17, 

13 

18. 

13 

19,11 

15, 

10 

14, 

9 

13, 

9 

13, 

9 

13, 

9 

1«. 

9 

14. 

9 

71 

19,13 

16, 

11 

15, 

11 

15, 

10 

15. 

10 

15, 

10 

15. 

11 

16. 

11 

19,14 

17, 

13 

16, 

12 

16, 

12 

17, 

13 

17, 

13 

18, 

13 

19. 

14 

22,  13 

17, 

11 

16, 

10 

15, 

9 

15, 

9 

15, 

9 

15. 

9 

16, 

10 

73 

22,14 

18, 

12 

17, 

11 

16, 

11 

16, 

11 

16, 

11 

17. 

11 

18. 

12 

22,16 

19, 

14 

18, 

13 

18, 

13 

18, 

13 

18, 

14 

19, 

14 

20, 

15 

50,20 

20, 

12 

18, 

11 

16, 

10 

16, 

10 

16, 

10 

16, 

10 

17. 

10 

75 

27,16 

21, 

14 

19, 

13 

18, 

12 

18, 

12 

18, 

12 

18, 

12 

19, 

13 

26,18 

22, 

16 

20, 

15 

20, 

14 

19, 

14 

20, 

15 

21, 

15 

22, 

16 

BASEMENT 
AGES  USED 

63,21 
58,  23 
52,25 
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INTERSECTIONS  FOR  HOLES  72/80,  CONT*D 


RATE=  1.00 

LONG  -110  -100  -90  -80  -70  -60  -50  -40 

LET 


11 

. 8 

10, 

7 

9, 

7 

9. 

7 

9, 

7 

10. 

7 

10. 

7 

11. 

7 

65 

11 

, 8 

10, 

6 

10, 

8 

10. 

8 

10. 

8 

11, 

8 

11, 

8 

12, 

» 

11 

, 9 

11, 

9 

11. 

9 

11. 

9 

11, 

9 

12, 

10 

13. 

10 

14. 

11 

12 

, 8 

11. 

7 

10, 

7 

10, 

7 

10, 

7 

10, 

7 

10, 

7 

11. 

8 

67 

12 

, 9 

11, 

8 

11. 

8 

11, 

8 

11, 

8 

11, 

8 

12, 

9 

13, 

9 

13 

,10 

12, 

10 

12, 

10 

12, 

10 

12, 

10 

13. 

10 

14, 

11 

15. 

11 

50 

,20 

12, 

8 

11. 

7 

10, 

7 

10, 

7 

11, 

7 

11. 

8 

12, 

8 

69 

13 

,10 

12, 

9 

12, 

9 

11, 

9 

11. 

9 

12, 

9 

12. 

9 

13. 

10 

14 

,11 

13, 

10 

13, 

10 

13, 

10 

13, 

10 

14. 

11 

14, 

11 

15, 

12 

** 

>** 

13, 

8 

12. 

8 

11. 

8 

11, 

8 

11, 

8 

12. 

8 

12, 

8 

71 

15 

,11 

13, 

10 

13, 

9 

12, 

9 

12, 

9 

13, 

9 

13, 

10 

14, 

10 

15 

, 12 

14, 

11 

14. 

11 

14, 

11 

14, 

11 

15, 

11 

15, 

12 

16, 

12 

** 

14, 

9 

13, 

9 

12, 

8 

12, 

8 

12, 

8 

13, 

8 

13. 

9 

73 

** 

»** 

15, 

11 

14, 

10 

13, 

10 

13, 

10 

14, 

10 

14. 

10 

15. 

1 1 

18 

,13 

16, 

12 

15. 

12 

15, 

11 

15, 

12 

16, 

12 

16, 

12 

17, 

13 

** 

,** 

16, 

10 

15, 

9 

14, 

9 

13, 

9 

1«, 

9 

14, 

9 

14, 

9 

75 

** 

!** 

17, 

12 

16, 

11 

15, 

10 

15, 

10 

15, 

11 

16. 

11 

16, 

11 

46 

,25 

18, 

13 

17, 

13 

16, 

12 

16, 

12 

17. 

13 

18, 

13 

19, 

14 

BASEHENT 
AGES  USED 

63,21 
58,  23 
52,25 
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FILE  2 (KIND=.PRINTER) 
FILE  4 (KIND=READER) 


********************************************************************* 

THIS  PROGRAM  OSES  THE  METHOD  OF  FALSE  POSITION  TO  FIND  THE 
INTERSECTION  OF  A PAIR  OF  LATITUDE-DEPTH  TRACKS.  IT  DOES  THIS 
FOR  AN  INPUTTED  RANGE  OF  POLE  RATES,  LATITUDES  AND  LONGITUDES. 

IT  REPEATS  THIS  FOR  BASEMENT  AGES  WHICH  DIFFER  FROM  THE  INPUTTED 
AGES  BY  AN  ERROR  FACTOR,  AGERR.  IN  ONE  CASE  THE  BASEMENT  AGE  OF 
ONE  OF  THE  HOLES  IS  HIGH  AND  THE  OTHER  HOLE  IS  LON;  THE  SECOND 
CASE  IS  VICE  VERSA. 

THE  PRINT  OUT  LOOKS  LIKE  THIS: 

1 

17,25  +ERROR, -ERROR 

15,23  NO  ERROR.NO  ERROR 

12,20  -ERROR, ♦ERROR 

1 

THE  NUMBERS  REPRESENT  THE  TIMES  OF  INTERSECTION  (IN  MILLIONS  OF 
YEARS)  OF  THE  FIRST  AND  SECOND  HOLES  INPUTTED,  RESP. 

1 

NOTES: 

1)  WHEN  AN  INTERSECTION  DOES  NOT  EXIST  FOR  THE  GIVEN  TRACKS, 

HE  PRINT  OUT  **,•*.  1 

2)  WHEN  THE  NUMBER  OP  ITERATIONS  IN  THE  SOLUTION  FOR  THE 
INTERSECTION  EXCEEDS  "LIMIT",  HE  PRINT  OUT  99,99. 

3)  BECAUSE  OF  LINE  LENGTH  LIMITATIONS,  DO  NOT  ALLOW  THE 
NUMBER  OF  POSSIBILITIES  OF  LONGITUDE  FOR  THE  POLE  TO 
EXCEED  17. 

£******* ********** ******** ******** ******* ********* ****** *************** 

DIMENSION  XLAO  (2)  , 7.0  (2)  ,TB (2)  , A (2)  , R1  (2)  ,R2  (2)  , R3  (2)  ,XLO0  (2) 

1 , NHOLE  (2)  ,L  (2)  ,TE  (3.2)  ,AA  (3,2)  ,ITA  1(3,  20)  , I TA2 ( 3.  20) 

REAL  LO 

COMMON  PO,LO,TXX,PLA , PLO, XLAP, XLOPO, PI , Cl, C2 
C STATEMENT  FUNCTIONS 

T (J, Z) =TB  ( J)  - ( (Z-A  (J)  ) /B)  **2  XT  IS  AGE  BY  TREHU  FORMULA 
D(Z)  =R1(2)-R1(1)  + R 2 ( 2)  *COS  (R  3 (2)  *W*T  (2.Z)  )-R2  (1)*C0S  (R3  (1) 

1 ♦W*T(1,Z)) 

XD  IS  THE  DIFFERENCE  IN  THE  SINE  OF  LATITUDE  BETWEEN  TWO 
XTRACKS  AT  DEPTH  Z. 

C CONSTANTS 

PI=3. 1415926536 
DTR=PI/1 00. 

8=348.  XTOPOGRAPHY  CONSTANT 


n n 


V 
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AGERR=.1  XERROR  IN  AGES 


C INPUT  RANGE  AND  INCREHENTS  FOR  POLE  RATE,  LAT,  AND  LONG 
READ  (4,  108)  Hi,  W2,DW,PLAl,PLA2,DPLA,PL0l,PL02,DPL0,LIHIT 
IF  (LIMIT. EQ. 0)  L I MIT3 100 
108  FORMAT(11(F4.0,1X)) 

ILOl  = PL01  ♦SIGN  (.  0000  1.PL01) 

IL02=PL02  ♦SIGN  (.00001, PL02) 

ID LO*D PL 0+. 000001 
IN=IABS(IL02>IL0  1) /IDLO*1 
C CONVERT  DEGREES  TO  RADIANS 

W 1=DTF*W  1 
W2=DTR*W2 
DW=DTR*DW 
PLA1=DTR*PLA1 
PLA2=DTR*PLA2 
DPLA=PTR*DPLA 
PL0 1 = DTK*PL0  1 
PL02=D'i'R*PL02 
DPLO=DTR*DPLO 

C INPUT  INITIAL  LAT,  LONG,  DEPTH,  AND  AGE  BASBflENT 

1 DO  2 1*1,2 

READ  (4,100)  XL  AO  (I)  , XLOO  (I)  ,TB  (I)  ,Z0  (I)  ,NH0LE(I) 

100  F0RKAT(3F10.2,20X,F10.2, IX, 15  ) 

IF  (ZO  (I)  .EQ.  0.)  CALL  EXIT 

2 CONTINUE 


C OUTPUT  HEADING 

WRITE  (2,110) 

110  FORI!  AT  ( 1 H 1) 

WRITE  (2, 106) 

106  FORMAT (28X,  ‘VARIATION  OF  TRACK  INTERSECTIONS*/ 

1 24X,  ‘AS  A FUNCTION  OF  POLE  POSITION  AND  RATE'/ 

1 30X , ' AN D OF  ERROR  IN  BASEMENT  AGE’/) 

C CALCULATE  RANGE  OF  BASEMENT  AGES. 

DO  1000  ID=1 , 3 

TE  (ID,  1)  *TB(  1)  ♦ ( 1.-AGEF REFLOAT  (ID- 2)  ) 

TB (ID, 2) =TB (2) ♦ ( 1 . ♦ AGERR* FLO AT (I D-2) ) 

WE  USE  MORGAN'S  CURVE  FOR  TOPOGRAPHIC  SUBSIDENCE.  ARRAX  AA  IS  THE 
CONSTANT  DEPTH  CORRECTION. 
j DO  1000  1*1,2 

I AA (ID, I) =Z0 (I) -B*SQRT (TE (ID, I) ) 

1 1000  CONTINUE 


DO  10  1=1,2 


noon 
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WRITE {2, 101)  N HOLE  (I) ,ILA0  (I) , XLOO <I),Z0  <I) 

101  FORMAT(21X,*HOLE,,I4,3X»*LAT',F6.2,2X,‘LONG  ',F7.2,2X, 
1 ‘DEPTH  • , 14) 

C DEGREES  TO  RADIANS 

IU0  (I) =DTR*XLA0  (I) 

XLOO (I) = DTR*  XLOO  (I) 

C FIND  NUMBER  OF  DIGITS  IN  HOLE  NUMBER 

DO  5 N 3 1 , 1 0 

IF(NH01.E(I) /10**N  .EQ.  0)  GO  TO  B 
5 CONTINUE 
8 L ( I)  -N 
10  CONTINUE 


BEGIN  CALCULATIONS.  LOOP  THROUGH  POSSIBLE  ROATATIOH  RATES, 
POLE  LATITUDES,  BASEMENT  AGES,  AND  POLE  ROTATION  LONGITUDES. 


DO  6500  W=H1,H2‘. 00000001, DH 

C WRITE  COLUMN  HEADING 

M=0 

IF  (W.NE.W1)  GO  TO  1200 

WBITE(2,102)  W/DTR*.0001,IN, (I, I=IL01, IL02, IDLO) 

102  FORMAT (//3SX, ' RATE*' ,F5. 2// 

1 UX,‘  LONG*  ,3X,*  (IX,  14,  IX)/ 

1 10X , • LAT ' ) 

GO  TO  1250 
C NEK  PAGE 

12C0  WRITE  (2,110) 

WRITE  (2,  107)  (L  ( I)  , N HOLE  (I)  ,1=0.2) 

107  FORMAT  (10X,‘ INTERSECTIONS  FOR  HOLES • , 1 X, I*, •/' . I* , 8H,  CONT'D//) 
WRITE ( 2, 102)  W/DTR  + . 0001 ,IK.  (1, 1 = 1 L01 , IL02, I DLO) 

1250  DO  6000  PLA=PLA1 ,PLA2*. 00000001, DPLA 

K=0 
Mx(1*  1 

DO  5000  PLO=  PL01 ,PL02 '.0000001, DPLO 
KXK»1 


DO  2000  1=1,2 
L0  = XLO0  ( I) 

P0=XLA0  (I) 

CALL  GRIND  %GR IN  D DOES  MOST  OF  THE  WORK 

Rl  (I)=COS  (XLAP)*SIN(Pt,A) 

B2  (I)  =511.  ( XU  AP)  *COS  ( PLA) 

R3  (I) =XLOPO 


mi  iiii 


li A 
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2000  CONTINUE 

DO  SCOO  10*1,3  XAGE  E8BOH  LOOP 

DO  3000  1=1,2 
TB(I)=XE  (ID, I) 

A (I) =AA (ID, I) 

3000  CONTINUE 

C USE  HETUOD  OP  PALSE  POSITION  TO  FIND  CROSSING. 

C FIND  PI3ST  TWO  GUESSES:  THE  LIHITS  OP  THE  TRACKS. 

Z 1 — A f 5 A X 1 (A  (1 ) , A (2) ) XAMAXl  IS  SISTEH  INTRINSIC  TO  PIND  MAX 

22= A HI  N 1 (ZO  ( 1)  ,Z0(2))  XAHIN1  IS  SISTER  INTRINSIC  TO  TIND  HI N 
D1  = D(Z1) 

D2^D  (2  2) 

C IS  THERE  A CROSSING? 

IF (SIGN ( 1. , D 1) . EQ. SIGN ( 1 . ,D2) ) GO  TO  3700 

KOUNT*0 

3100  KOUNT=KOUNT*1 

IP  (ABS <D1) .LT..001)  GO  TO  3801 
IP(ABS(D2)  .LT..001)  GO  TO  3802 
IP  (KOUNT.GT. LI  HIT)  GO  TO  3803 
Z3=ZH(Z1-Z2)/(D2-D1)  *D1  XNEN  GUESS 

D3  = D ( Z 3 ) 

IF  (S IGH ( 1. , 0 3)  . EQ. SIGN  { 1 . ,02)  ) GO  TO  3200  XWHICH  SIDE? 

Z1  = Z3 
D 1 “D  3 

GO  TO  3100 
3200  Z2=2  3 
D2  = D3 

GO  TO  3100 

3700  IT1=1.B9 
IT2= 1. E9 
GO  TO  3850 
380  1 IT  1=T  ( 1,  Z 1)  ♦.  5 
IT2=T  (2 , Z1 ) ♦ .5 
GO  TO  3850 

3802  I T1  = T (1 ,Z2) ♦ . 5 
IT  2=T (2,2  2) ♦ . 5 
GO  TO  3850 

3803  IT  1 = 99 
IT2=99 

3850  ITA  1 (ID, K)  =IT 1 
ITA2  (ID,  K)  =IT2 


%NO  CROSSING 

^CROSSING 

^CROSSING 

XCOULDN*  T PIND  CROSSING  PAST  ENOUGH 


- . • . 


non 
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5000  CONTINUE 
OUTPUT 


IF  (H.LT.  9)  GO  TO  5100 
C NEW  PAGE 

WRITE  (2,110) 

WRITE (2,  107)  (L  (I)  , NHOLE  (I)  ,1*1,2) 

WRITE  (2 , 102)  W/DTR*.0001,IN,  (I,I*IL01, IL02, IDLO) 

5100  DO  550Q  ID*  1 , 3 

GOTO  (5110,5120,5130) , ID 

5110  WRITE(2, 109)  IN,  (ITA 1 (ID, IK)  ,ITA2(ID,IK) ,IK*1,IN) 

104  FORHAT  (20X,*  ( 1 X, 12  , • , • , 12) ) 

GO  TO  5500 

5120  IQ=PLA/DTR  ♦ SIGN  (.000001, PLA) 

WRITE! 2, 10  3)  IQ, IN,  (ITA 1 (ID, IK) ,ITA2  (ID , I K) , I K* 1 ,IN) 
103  FORHAT  (10 X,1 3,7 X,*  (1 X , 12 , * , • , 12) ) 

GO  TO  5500 

5130  URITE (2,105)  I N,  (ITA 1 (ID, IK) , ITA2 (ID, IK) , IK*  1, IN) 

105  FORHAT (20X,* (IX, 12,*,*, 12)/) 

5500  CONTINUE 
6000  CONTINUE 

WRITE  (2,117) 

117  FORHAT  (//  40X,'BASEHENT'/40X,'AGES  USED'/) 

DO  6300  ID* 1,3 

J=2 

IF  (TE  (ID,  2)  .GE.  100.)  J*3 

WRITE  (2, 11  8)  TE(ID,1)*.5,J,TE(ID,2>  ♦.S 

118  FORHAT (40X, IX, 13,' ,' ,1*) 

6300  CONTINUE 

6500  CONTINUE 

GO  TO  1 
END 


: a.  . 
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SUBROUTINE  GRIND 
REAL  LO, LPO, LX.LP1 

COMMON  PO, LO,TH, PS, LX, PH, LPO ,PI,C1,C2 

C THIS  SUBROUTINE  USES  THE  HALF  ANGLE  FORMULAS  FROM  SPHERICAL 
C TRIG  TO  CALCULATE  THE  LATITUDE  AND  LONGITUDE  OF  A POINT, 

C LO, PO,  IN  COORDINATES  OF  THE  POLE.  THE  POLE  IS  AT  LX, PS.  THE 

C LATITUDE  OF  THE  POINT  IN  PLATE  POLE  COORDINATES  IS  PH,  THE  LONG I- 

C TUDE  IS  LPO. 

C3*SIN(.5*(LX-L0  ) ) 

C5*SIN (. 5* (PS* PO  )) 

IF (CS. EQ.O.O)  CS=1.0E-9 
C63.C0S  (.5*  (PS*  PO  )) 

IF(C6. EQ.O.O)  CO-1 .OE-9 
C7*SI N (. 5*  (PS- PO  )) 

C8=C0S (.5* (PS-PO  )) 

C9*C6/C5 

C10*COTAN{.5* (LX-LO  )) 

W1*C10*C8/CS 
W2*C10*C7/C0 
C11*  ATAN(WI) 

IF (C 1 1.LT. 0. ) C11=C11»PI 
C12*  ATAN  (H2) 

IF (C12.LT. 0. ) C12*C12*PI 
C13*C0S  (Cl  1) 

C14»C0S(C12) 

IF  (CU. EQ.O.O)  Cl4*1.0E-9 
U=C9*C 13/C  14 
U*ATAN  (U) 

IF (U .LT. 0.  ) U»U*PI 

PH=*2.*U 

LP0*C1  UC12 

IF  (LPO. GT. PI)  LP0=2.*PI-LP0 

RETURN 

END 
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